Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_014034 1,303,879 C→A 3.6% intergenic (+66/+91) urtE → / ← RCAP_rcc01231 urea ABC transporter ATP‑binding protein UrtE/UbiH/UbiF/VisC/COQ6 family ubiquinone biosynthesis hydroxylase

Read alignment evidence...
  seq id position change freq score reads annotation genes product
*NC_0140341,303,8790C→A 3.6% 3.0 ‑3.7 ‑0.1110intergenic (+66/+91)urtE/RCAP_rcc01231urea ABC transporter ATP‑binding protein UrtE/UbiH/UbiF/VisC/COQ6 family ubiquinone biosynthesis hydroxylase
Reads supporting (aligned to +/- strand):  new base (4/0):  ref base (41/65):  total (45/65)
Fisher's exact test strand distribution p-value = 2.58e-02
Kolmogorov-Smirnov test that lower quality scores support polymorphism than reference p-value = 7.39e-01

Only 100 of 121 total aligned reads displayed.

GCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCT.C..CGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCC....GCCAAGATGGGG..G..GC.CCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCA  >  NC_014034/1303782‑1303972
                                                                                                       |                                                                                                      
gcgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCt                                                                                                       <  1:3385630/101‑1 (MQ=37)
gcgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCt                                                                                                       <  1:1345128/101‑1 (MQ=37)
gcgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCt                                                                                                       >  1:324403/1‑101 (MQ=37)
gcgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCt                                                                                                       <  2:3374429/101‑1 (MQ=37)
 cgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCtt                                                                                                      <  2:4375600/101‑1 (MQ=37)
 cgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCtt                                                                                                      >  2:3809545/1‑101 (MQ=37)
 cgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCtt                                                                                                      <  2:3753589/101‑1 (MQ=37)
 cgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCtt                                                                                                      <  2:3151845/101‑1 (MQ=37)
 cgAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCtt                                                                                                      <  1:2985676/101‑1 (MQ=37)
   aaaCCCTGCGCGGCGTGGTCAGCGTTTAGGCAAGGGGGCGCTGCCCCC.C..CGGCTTCGCCTCACCCCCGGGCGATTTGGGGCAAGAT.AAACCCCATCCTTTc                                                                                                    >  1:2520620/1‑101 (MQ=37)
   aaaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTc                                                                                                    >  2:672659/1‑101 (MQ=37)
   aaaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTc                                                                                                    <  1:4546506/101‑1 (MQ=37)
   aaaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTc                                                                                                    <  2:133684/101‑1 (MQ=37)
   aaaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTc                                                                                                    >  2:2915590/1‑101 (MQ=37)
    aaCCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTcc                                                                                                   >  2:2933032/1‑100 (MQ=25)
      cccTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCAcc                                                                                                 >  1:3069401/1‑101 (MQ=37)
      cccTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCAcc                                                                                                 >  2:762776/1‑101 (MQ=37)
       ccTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCt                                                                                                >  2:1356017/1‑101 (MQ=37)
        tgcagcgcCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCtt                                                                                               >  2:2187067/5‑101 (MQ=25)
        cTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCtt                                                                                               <  1:2513536/101‑1 (MQ=37)
        cTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCtt                                                                                               <  1:3809830/101‑1 (MQ=37)
         tGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTg                                                                                              >  2:2185049/1‑101 (MQ=37)
          gcgcgcCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGc                                                                                             <  2:4118282/101‑1 (MQ=37)
           cgcggCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGcc                                                                                            <  2:2964949/101‑1 (MQ=25)
             cgcCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGc                                                                                          <  1:489274/101‑1 (MQ=37)
              gcCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCa                                                                                         <  2:3567447/101‑1 (MQ=37)
               ccGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCaa                                                                                        <  1:2874406/101‑1 (MQ=37)
               ccGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCaa                                                                                        >  1:1585085/1‑101 (MQ=37)
                cGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCCCCCCCCGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCaaa                                                                                       >  1:777301/1‑101 (MQ=37)
                 gTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGGTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAat                                                                                      >  2:2545248/1‑101 (MQ=37)
                 gTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAat                                                                                      <  2:2241231/101‑1 (MQ=37)
                 gTGGTCAGCGTTTAGGCAAGGGAGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAat                                                                                      <  2:4474152/101‑1 (MQ=37)
                  tGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAata                                                                                     <  2:4252110/101‑1 (MQ=37)
                  tGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAata                                                                                     <  1:2986212/101‑1 (MQ=37)
                   gggcAGGGTTTAGGGCAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAatat                                                                                    <  1:1531611/98‑1 (MQ=21)
                   ggTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAatat                                                                                    <  1:1129550/101‑1 (MQ=37)
                   ggTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAatat                                                                                    <  1:4116336/101‑1 (MQ=37)
                    gTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATc                                                                                   <  2:568293/101‑1 (MQ=37)
                    gTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATc                                                                                   <  1:2185199/101‑1 (MQ=37)
                    gTCAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATc                                                                                   >  2:3435157/1‑101 (MQ=25)
                      cAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATccc                                                                                 >  1:3110031/1‑101 (MQ=37)
                      cAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATccc                                                                                 <  1:664412/101‑1 (MQ=37)
                      cAGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATccc                                                                                 <  2:3555021/101‑1 (MQ=37)
                       aGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCg                                                                                <  1:3876810/101‑1 (MQ=37)
                       aGCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCg                                                                                >  1:3606527/1‑101 (MQ=25)
                        gCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCgg                                                                               <  1:442255/101‑1 (MQ=37)
                        gCGTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCgg                                                                               <  1:186243/101‑1 (MQ=37)
                          gTTTAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCgggg                                                                             <  1:447572/101‑1 (MQ=37)
                          ctgacggacggggggCTCTGC.CCC.C..CGGCTTCG.CTC.CCCCCGGGATATTTGCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCaggggg                                                                           >  2:2585515/11‑101 (MQ=2)
                          ctgacggacggggggCTCTGC.CCC.C..CGGCTTCG.CTC.CCCCCGGGATATTTGCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCaggggg                                                                           >  1:3167541/11‑101 (MQ=2)
                           tttAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCggggg                                                                            >  1:3705785/1‑101 (MQ=37)
                            ttAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCgggggg                                                                           <  2:4448374/101‑1 (MQ=37)
                             tAGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCCAATATCCCggggggc                                                                          >  2:567525/1‑100 (MQ=25)
                              aGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCgggggggg                                                                         >  1:541675/1‑99 (MQ=25)
                              aGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTg                                                                         <  1:2882259/101‑1 (MQ=37)
                              aGGCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTg                                                                         <  1:457214/101‑1 (MQ=37)
                               ggCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGAAAAAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGa                                                                        <  1:4653332/101‑1 (MQ=37)
                                gCAAGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGaa                                                                       <  1:1192725/101‑1 (MQ=37)
                                  aaGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAAtt                                                                     >  2:3190409/1‑101 (MQ=25)
                                  aaGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAAtt                                                                     <  2:4271463/101‑1 (MQ=37)
                                  aaGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAAtt                                                                     <  1:762851/101‑1 (MQ=37)
                                  aaGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCGATCACCTTGCCGCAAATATCCCGGGGGGTGAAtt                                                                     >  1:4005449/1‑101 (MQ=25)
                                   aGGGGGCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTg                                                                    <  1:3012741/101‑1 (MQ=37)
                                    gggggCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTgg                                                                   >  1:3012954/1‑101 (MQ=37)
                                      gggCTCTG..CCCCC..CGGCTTCG.CTC.CCCCCGGGATATTTGCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCCAGGGGG.GAATTGgccaaa                                                              >  1:3204317/1‑98 (MQ=2)
                                      gggCTCTG..CCCCC..CGGCTTCG.CTC.CCCCCGGGATATTTGCGGCAAGATGAAGGGGCATGATTTCATCTTGCCCGAAATATCCC.AGGGGTGAATTGgccaaa                                                              >  1:2983703/1‑98 (MQ=2)
                                      gggCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGgcc                                                                 >  2:4053354/1‑101 (MQ=37)
                                       ggCGCTG..CCC.CCTCGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGgcc....g                                                            <  2:3069157/101‑1 (MQ=37)
                                                                                 tGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGgggcggg                       <  1:2480518/101‑1 (MQ=14)
                                                                                 tGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGgggcggg                       >  2:169142/1‑101 (MQ=14)
                                                                                  ggggCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGgggcgggc                      <  2:1302747/101‑1 (MQ=14)
                                                                                   gggCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggcc                     >  2:4068038/1‑101 (MQ=14)
                                                                                       aaGAT.AACCCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcag                 <  2:324371/101‑1 (MQ=9)
                                                                                       aaGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcag                 <  2:3513518/101‑1 (MQ=14)
                                                                                       aaGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcag                 >  1:3519800/1‑101 (MQ=14)
                                                                                       aaGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGGGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcag                 >  1:3361931/1‑101 (MQ=9)
                                                                                        aGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcagg                >  1:4163937/1‑101 (MQ=14)
                                                                                        aGAT.AAACCCCATCCTTTCACCTTGCCGCAAATACCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcagg                <  2:2228843/101‑1 (MQ=9)
                                                                                         gAT.AAACCCCATCCTTTCCCCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcaggg               <  1:696988/101‑1 (MQ=9)
                                                                                         gAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCCGCGCCCCCTTTTCGGGGCgggccgccagg               >  2:2901542/1‑97 (MQ=9)
                                                                                          aT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGGGAATTGGCCGAAGGCCAGGAGGGGG..GCAGCGCCCCCTTTTAGGGGCgggccgcagggc              >  1:1387025/1‑101 (MQ=9)
                                                                                           t.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcagggcc             <  2:3334533/101‑1 (MQ=14)
                                                                                           t.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcagggcc             >  1:3835737/1‑101 (MQ=14)
                                                                                             aagccctgtCCTTTCGCCTTGCCGCAAAGGTTCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcagggccg            <  1:2545422/93‑1 (MQ=2)
                                                                                             aaaCCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcagggccg            >  2:1990618/1‑101 (MQ=14)
                                                                                             aaaCCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcagggccg            >  2:235268/1‑101 (MQ=14)
                                                                                              aaCCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcagggccgc           >  1:951809/1‑101 (MQ=14)
                                                                                              aaCCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcagggccgc           >  2:2594452/1‑101 (MQ=14)
                                                                                               aCCCCATCCTTTCTCCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCgggccgcagggccgca          <  1:2532924/101‑1 (MQ=9)
                                                                                                ccccATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCAc         <  2:211985/101‑1 (MQ=14)
                                                                                                ccccATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCAc         >  1:625181/1‑101 (MQ=14)
                                                                                                ccccATCCTTTCACCTTGCCGCAAATATCCCGGGGGGGGAATTGGCCGAAGGCCAAGAGGGGGTCA..GCGCGCCCTTTTCGGGGCGGGCCGCAGGGCCGCAc         >  2:2830202/1‑101 (MQ=9)
                                                                                                ccccATCCTTTCACCTTGCCGCAAATATCCCGGGGGGGGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGGACCCTTTTCGGGGCGGGCCGCAGGGCCGCAc         >  1:3428062/1‑101 (MQ=9)
                                                                                                 cccATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCAcc        >  2:4394458/1‑101 (MQ=14)
                                                                                                   cATCCTTTCACCTTGCCGCAAATCTCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCt      <  1:1934072/101‑1 (MQ=9)
                                                                                                    ctCCTTTCTCGTTGCGGTATATATCCTGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGT..GCAGCGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCCCCCTg     <  2:4553359/100‑1 (MQ=2)
                                                                                                      ccTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTc   >  1:244201/1‑101 (MQ=14)
                                                                                                      ccTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTc   >  2:2763581/1‑101 (MQ=14)
                                                                                                       cTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCa  <  2:3417652/101‑1 (MQ=14)
                                                                                                       cTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCCGAAGGCCAAGAGGGGG..GCAGCGCCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCa  >  2:3139949/1‑101 (MQ=14)
                                                                                                       |                                                                                                      
GCGAAACCCTGCGCGCCGTGGTCAGCGTTTAGGCAAGGGGGCGCTG..CCT.C..CGGCTTCGCCTCACCCCCGGGATATTTGGGGCAAGAT.AAACCCCATCCTTTCACCTTGCCGCAAATATCCCGGGGGGTGAATTGGCC....GCCAAGATGGGG..G..GC.CCCCCTTTTCGGGGCGGGCCGCAGGGCCGCACCCTGTCA  >  NC_014034/1303782‑1303972

Base quality scores:  ATCG  < 3 ≤  ATCG  < 16 ≤  ATCG  < 25 ≤  ATCG  < 33 ≤  ATCG  < 41 ≤  ATCG