Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,880,731 G→T 92.7% D182Y (GAT→TAT)  mnmE → tRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092733,880,7310GT92.7% 595.8 / 35.4 215D182Y (GAT→TAT) mnmEtRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase
Reads supporting (aligned to +/- strand):  ref base G (9/7);  new base T (92/106);  total (102/113)
Fisher's exact test for biased strand distribution p-value = 6.04e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 216 total aligned reads displayed.

CGACGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCGATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTAAATCGA  >  CP009273/3880586‑3880878
                                                                                                                                                 |                                                                                                                                                   
cGACGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTAt                                                                                                                                                   <  2:1346009/148‑1 (MQ=255)
cGACGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCGAt                                                                                                                                                   <  2:1015763/148‑1 (MQ=255)
  aCGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTAttt                                                                                                                                                 <  1:461126/148‑1 (MQ=255)
    gCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCGATTTcc                                                                                                                                               >  2:1764439/1‑148 (MQ=255)
     ccAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCct                                                                                                                                              >  1:2631216/1‑148 (MQ=255)
     ccAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCct                                                                                                                                              >  1:2631186/1‑148 (MQ=255)
         ttCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTcc                                                                                                                                          >  1:2943099/1‑148 (MQ=255)
             gAACAGGCGGCCCGTTCGGAACTTAACTCGCTGCAAGGCGCATTCTCCGCTCTGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTAAGTAGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACg                                                                                                                                      <  2:1959143/148‑1 (MQ=255)
              aaCAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACgg                                                                                                                                     <  2:2191605/148‑1 (MQ=255)
               aCAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGa                                                                                                                                    >  1:1864620/1‑148 (MQ=255)
                   gcggcCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGaaaaa                                                                                                                                >  2:465772/1‑148 (MQ=255)
                   gcggcCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGaaaaa                                                                                                                                <  1:2799350/148‑1 (MQ=255)
                   gcggcCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCGATTTCCTCTCCGACGGaaaaa                                                                                                                                >  1:2416139/1‑148 (MQ=255)
                    cggcCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAAt                                                                                                                               <  2:2504059/148‑1 (MQ=255)
                          gTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGc                                                                                                                         <  2:502261/148‑1 (MQ=255)
                           ttCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGcc                                                                                                                        >  2:2693945/1‑148 (MQ=255)
                                cACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCt                                                                                                                   >  1:1140883/1‑148 (MQ=255)
                                cACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCt                                                                                                                   >  1:2687077/1‑148 (MQ=255)
                                   ttAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGAGTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCaa                                                                                                                <  2:793781/148‑1 (MQ=255)
                                    tAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAAt                                                                                                               <  2:1845419/148‑1 (MQ=255)
                                     aaCTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATGTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATg                                                                                                              <  1:204111/148‑1 (MQ=255)
                                      aCTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGa                                                                                                             <  2:1326188/148‑1 (MQ=255)
                                      aCTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGa                                                                                                             <  2:1326178/148‑1 (MQ=255)
                                      aCTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGa                                                                                                             >  1:1015375/1‑148 (MQ=255)
                                        tCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACg                                                                                                           <  1:611857/148‑1 (MQ=255)
                                        tCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACg                                                                                                           <  1:1319053/148‑1 (MQ=255)
                                          gCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGtt                                                                                                         <  1:1848545/148‑1 (MQ=255)
                                           cTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTa                                                                                                        >  2:861663/1‑148 (MQ=255)
                                              cAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTg                                                                                                     <  1:3004068/148‑1 (MQ=255)
                                              cAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTg                                                                                                     <  1:1357135/148‑1 (MQ=255)
                                                 ggCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCg                                                                                                  >  2:1938623/1‑148 (MQ=255)
                                                 ggCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCg                                                                                                  >  1:1811608/1‑148 (MQ=255)
                                                 ggCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCg                                                                                                  >  1:1431355/1‑148 (MQ=255)
                                                  gcgcATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGa                                                                                                 >  2:1534798/1‑148 (MQ=255)
                                                  gcgcATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGa                                                                                                 >  2:2020850/1‑148 (MQ=255)
                                                  gcgcATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGa                                                                                                 <  1:349897/148‑1 (MQ=255)
                                                   cgcATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGAt                                                                                                <  1:1596153/148‑1 (MQ=255)
                                                   cgcATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGAt                                                                                                <  2:1652069/148‑1 (MQ=255)
                                                    gcATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATc                                                                                               <  1:2885944/148‑1 (MQ=255)
                                                        tctcCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGa                                                                                           <  2:1089995/148‑1 (MQ=255)
                                                        tctcCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGa                                                                                           <  2:448005/148‑1 (MQ=255)
                                                        tctcCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCGTGa                                                                                           >  2:2798233/1‑148 (MQ=255)
                                                           ccGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGc                                                                                        >  1:512028/1‑148 (MQ=255)
                                                           ccGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCGATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGc                                                                                        >  1:2055812/1‑148 (MQ=255)
                                                            cGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCCATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCa                                                                                       >  1:285583/1‑148 (MQ=255)
                                                              cACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAgt                                                                                     >  1:1271816/1‑148 (MQ=255)
                                                                 gggTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAgtgcg                                                                                  <  2:3162698/148‑1 (MQ=255)
                                                                   gTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAgtgcgtg                                                                                <  2:2312046/148‑1 (MQ=255)
                                                                   gTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAgtgcgtg                                                                                <  2:1301762/148‑1 (MQ=255)
                                                                      aaTCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTg                                                                             >  1:903036/1‑148 (MQ=255)
                                                                       atcatcTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGa                                                                            <  1:708399/148‑1 (MQ=255)
                                                                           tcTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCa                                                                        >  1:1012696/1‑148 (MQ=255)
                                                                               gTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTc                                                                    >  1:3194933/1‑148 (MQ=255)
                                                                                  gAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAgg                                                                 >  1:1724106/1‑148 (MQ=255)
                                                                                     gCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTa                                                              >  1:517623/1‑148 (MQ=255)
                                                                                        cTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGtt                                                           >  1:1126282/1‑148 (MQ=255)
                                                                                         tCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGttt                                                          >  1:571992/1‑148 (MQ=255)
                                                                                          cACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTttg                                                         <  1:419687/148‑1 (MQ=255)
                                                                                           aCCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCt                                                                              <  1:2631513/126‑1 (MQ=255)
                                                                                           aCCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCt                                                                              >  2:2631513/1‑126 (MQ=255)
                                                                                               aCTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTgcg                                                    <  1:2586872/148‑1 (MQ=255)
                                                                                               aCTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTgcg                                                    >  2:2084649/1‑148 (MQ=255)
                                                                                                    cgcATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAg                                               >  1:2456620/1‑148 (MQ=255)
                                                                                                    cgcATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAg                                               <  1:1831687/148‑1 (MQ=255)
                                                                                                    cgcATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAg                                               <  1:229143/148‑1 (MQ=255)
                                                                                                      cATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAggg                                             <  2:3277445/148‑1 (MQ=255)
                                                                                                           aCGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGACGCACGTCAGGGTAGTGGGTTGCGCGAAGGGATGaa                                        >  2:3267896/1‑148 (MQ=255)
                                                                                                            cGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGaaa                                       <  1:465772/148‑1 (MQ=255)
                                                                                                            cGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGaaa                                       <  1:1546659/148‑1 (MQ=255)
                                                                                                               cGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAgtg                                    <  2:2942920/148‑1 (MQ=255)
                                                                                                                gAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAgtgg                                   <  1:2002822/148‑1 (MQ=255)
                                                                                                                 aaGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAgtggt                                  >  2:1726206/1‑148 (MQ=255)
                                                                                                                 aaGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAgtggt                                  >  2:350044/1‑148 (MQ=255)
                                                                                                                   gCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGa                                <  1:3067651/148‑1 (MQ=255)
                                                                                                                    cGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGGAGTTTGTTGCGCGAAGGGATGAAAGTGGTGAt                               >  2:753228/1‑148 (MQ=255)
                                                                                                                      gCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTg                             >  2:2612508/1‑148 (MQ=255)
                                                                                                                      gCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTg                             >  2:2615720/1‑148 (MQ=255)
                                                                                                                        aaTTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGcc                           <  1:2308788/148‑1 (MQ=255)
                                                                                                                        aaTTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGcc                           >  1:314114/1‑148 (MQ=255)
                                                                                                                           tGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGa                        >  2:651543/1‑148 (MQ=255)
                                                                                                                              tttCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGt                     <  2:1743864/148‑1 (MQ=255)
                                                                                                                                    cGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAAc               >  1:1194444/1‑148 (MQ=255)
                                                                                                                                    cGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAAc               >  1:60633/1‑148 (MQ=255)
                                                                                                                                    cGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAAc               >  1:2180956/1‑148 (MQ=255)
                                                                                                                                      aTGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGc             <  2:3114976/148‑1 (MQ=255)
                                                                                                                                      aTGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGc             <  2:1901277/148‑1 (MQ=255)
                                                                                                                                      aTGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGc             <  2:1475925/148‑1 (MQ=255)
                                                                                                                                         aaGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCgg          <  2:3044551/148‑1 (MQ=255)
                                                                                                                                         aaGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCgg          <  2:3044677/148‑1 (MQ=255)
                                                                                                                                          agagATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGt         >  1:745068/1‑148 (MQ=255)
                                                                                                                                          agagATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGt         >  1:745082/1‑148 (MQ=255)
                                                                                                                                           gagaTCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTa        >  2:1174396/1‑148 (MQ=255)
                                                                                                                                            agaTCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTaa       >  2:694216/1‑148 (MQ=255)
                                                                                                                                            agaTCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTaa       >  2:1447123/1‑148 (MQ=255)
                                                                                                                                            agaTCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTaa       >  2:605807/1‑148 (MQ=255)
                                                                                                                                            agaTCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTaa       <  2:1140883/148‑1 (MQ=255)
                                                                                                                                            agaTCGATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTaa       >  2:2997426/1‑148 (MQ=255)
                                                                                                                                              aTCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTAAAt     >  2:1245646/1‑148 (MQ=255)
                                                                                                                                              aTCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTAAAt     >  1:876912/1‑148 (MQ=255)
                                                                                                                                                 gATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTAAATCGa  <  1:1764439/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CGACGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCGATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTAAATCGA  >  CP009273/3880586‑3880878

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: