Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 3,880,731 G→T 100% 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,7310GT96.7% 650.0 / 11.5 211D182Y (GAT→TAT) mnmEtRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase
Reads supporting (aligned to +/- strand):  ref base G (6/1);  new base T (110/94);  total (116/95)
Fisher's exact test for biased strand distribution p-value = 1.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.95e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

Only 100 of 214 total aligned reads displayed.

CGACGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCGATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTAAATCGA  >  CP009273/3880586‑3880878
                                                                                                                                                 |                                                                                                                                                   
cGACGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTAt                                                                                                                                                   >  1:1954446/1‑148 (MQ=255)
   cGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTc                                                                                                                                                >  2:1607425/1‑148 (MQ=255)
      cAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCctc                                                                                                                                             <  1:2388080/148‑1 (MQ=255)
       aGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCctct                                                                                                                                            <  2:2035617/148‑1 (MQ=255)
       aGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCctct                                                                                                                                            <  2:126378/148‑1 (MQ=255)
         ttCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTcc                                                                                                                                          >  1:1610897/1‑148 (MQ=255)
           cGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGa                                                                                                                                        <  1:629172/148‑1 (MQ=255)
           cGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGa                                                                                                                                        <  1:2747355/148‑1 (MQ=255)
             gAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACg                                                                                                                                      <  2:3184187/148‑1 (MQ=255)
              aaCAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACgg                                                                                                                                     >  2:1412481/1‑148 (MQ=255)
                   gcggcCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGaaaaa                                                                                                                                >  1:1047811/1‑148 (MQ=255)
                   gcggcCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGaaaaa                                                                                                                                >  1:1050087/1‑148 (MQ=255)
                       cccGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGa                                                                                                                            >  1:594529/1‑148 (MQ=255)
                       cccGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGa                                                                                                                            >  1:2618388/1‑148 (MQ=255)
                       cccGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGa                                                                                                                            <  2:123243/148‑1 (MQ=255)
                       cccGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGa                                                                                                                            <  1:1959773/148‑1 (MQ=255)
                       cccGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGa                                                                                                                            >  2:1839154/1‑148 (MQ=255)
                        ccGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGaa                                                                                                                           <  1:2024576/148‑1 (MQ=255)
                        ccGTTCGGCACTTAACTCGCTGCAAGGAGCATTCTCCGCACGGGTGAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGaa                                                                                                                           <  1:223607/148‑1 (MQ=255)
                            tCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGccc                                                                                                                       >  2:3428158/1‑148 (MQ=255)
                             cGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCa                                                                                                                      >  1:1398212/1‑148 (MQ=255)
                             cGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCa                                                                                                                      >  1:2907810/1‑148 (MQ=255)
                               gCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGc                                                                                                                    >  1:619287/1‑148 (MQ=255)
                               gCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGc                                                                                                                    >  1:2411534/1‑148 (MQ=255)
                               gCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGc                                                                                                                    >  2:2040787/1‑148 (MQ=255)
                                     aaCTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATg                                                                                                              <  1:673898/148‑1 (MQ=255)
                                     aaCTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATg                                                                                                              <  1:1889422/148‑1 (MQ=255)
                                      aCTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACGTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGa                                                                                                             <  2:2176237/148‑1 (MQ=255)
                                        tCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACg                                                                                                           <  2:1849358/148‑1 (MQ=255)
                                          gCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGtt                                                                                                         <  1:820640/148‑1 (MQ=255)
                                           cTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTa                                                                                                        >  1:1698381/1‑148 (MQ=255)
                                             gCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTAtt                                                                                                      >  1:1777975/1‑148 (MQ=255)
                                              cAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGTGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTg                                                                                                     >  2:1049173/1‑148 (MQ=255)
                                                  gcgcATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGa                                                                                                 >  2:2969695/1‑148 (MQ=255)
                                                  gcgcATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGa                                                                                                 >  1:3511980/1‑148 (MQ=255)
                                                       ttgtcCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTg                                                                                            >  2:3349978/4‑148 (MQ=255)
                                                       ttCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTg                                                                                            >  2:3349702/1‑148 (MQ=255)
                                                       ttCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTg                                                                                            <  1:1013639/148‑1 (MQ=255)
                                                        tctcCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGa                                                                                           <  2:1374836/148‑1 (MQ=255)
                                                           ccGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGc                                                                                        <  2:2450091/148‑1 (MQ=255)
                                                                    ttAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAgtgcgtgc                                                                               <  2:1616307/148‑1 (MQ=255)
                                                                       atcatcTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGa                                                                            >  1:1330309/1‑148 (MQ=255)
                                                                          atcTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGc                                                                         <  2:1246361/148‑1 (MQ=255)
                                                                           tcTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCa                                                                        >  1:374857/1‑148 (MQ=255)
                                                                                 aGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAg                                                                  <  1:1176759/148‑1 (MQ=255)
                                                                                   aaGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAggg                                                                <  1:1296587/148‑1 (MQ=255)
                                                                                   aaGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAggg                                                                >  1:1599879/1‑148 (MQ=255)
                                                                                   aaGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAggg                                                                <  1:3237253/148‑1 (MQ=255)
                                                                                      cccTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGtcg                                                             >  1:43511/1‑146 (MQ=255)
                                                                                      cccTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAg                                                             >  2:1156746/1‑148 (MQ=255)
                                                                                       ccTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGt                                                            <  1:2493675/148‑1 (MQ=255)
                                                                                       ccTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGt                                                            <  2:2189748/148‑1 (MQ=255)
                                                                                       ccTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGt                                                            >  2:831735/1‑148 (MQ=255)
                                                                                        cTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGtt                                                           <  2:619287/148‑1 (MQ=255)
                                                                                          cACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGTGTAGTttg                                                         >  1:307811/1‑148 (MQ=255)
                                                                                          cACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTttg                                                         <  2:2931321/148‑1 (MQ=255)
                                                                                          cACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTttg                                                         <  1:600626/148‑1 (MQ=255)
                                                                                           aCCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTttgt                                                        >  2:3524652/1‑148 (MQ=255)
                                                                                             ccACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTttgttg                                                      >  2:2373701/1‑148 (MQ=255)
                                                                                             ccACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTttgttg                                                      >  2:583672/1‑148 (MQ=255)
                                                                                               aCTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTgcg                                                    <  1:1118706/148‑1 (MQ=255)
                                                                                                  tGCGCAGTTACGTCGCAGCGGCAATTGATGTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGa                                                 <  2:2899497/148‑1 (MQ=255)
                                                                                                    cgcATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAg                                               >  2:2799720/1‑148 (MQ=255)
                                                                                                       aTTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGa                                            >  2:1457346/1‑148 (MQ=255)
                                                                                                          tACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGACGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGa                                         >  2:319177/1‑148 (MQ=255)
                                                                                                           aCGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGaa                                        <  2:2908184/148‑1 (MQ=255)
                                                                                                           aCGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGaa                                        <  1:3356234/148‑1 (MQ=255)
                                                                                                            cGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGaaa                                       >  1:441896/1‑148 (MQ=255)
                                                                                                            cGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGaaa                                       <  2:710838/148‑1 (MQ=255)
                                                                                                            cGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGaaa                                       >  1:2267502/1‑148 (MQ=255)
                                                                                                            cGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGaaa                                       >  1:3000686/1‑148 (MQ=255)
                                                                                                            cGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGaaa                                       <  1:503045/148‑1 (MQ=255)
                                                                                                             gTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAg                                      >  2:2611318/1‑148 (MQ=255)
                                                                                                              tCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAgt                                     <  1:2125571/148‑1 (MQ=255)
                                                                                                              tCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAgt                                     <  2:1047811/148‑1 (MQ=255)
                                                                                                                  aGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAgtggtg                                 >  1:3203392/1‑148 (MQ=255)
                                                                                                                   gCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGa                                <  1:1515341/148‑1 (MQ=255)
                                                                                                                    cGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGAt                               <  1:1014022/148‑1 (MQ=255)
                                                                                                                    cGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTGTGTTGCGCGAAGGGATGAAAGTGGTGAt                               >  2:2851891/1‑148 (MQ=255)
                                                                                                                       cAATTGATTTCCCTGATGAAGAGATCGATTTCCTCTCCGACGGCAAAATTGAAGCCCAGCTTAATGACGTGATTGCCGATCTCGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAGGTGGTGATTGc                            >  1:2461961/1‑148 (MQ=255)
                                                                                                                       cAATTGATTTCCCTGATGAAGAGATCGATTTCCTCTCCGACGGCAAAATTGAAGCCCAGCTTAATGACGTGATTGCCGATCTCGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAGGTGGTGATTGc                            >  1:3056260/1‑148 (MQ=255)
                                                                                                                         aTTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCg                          >  2:2106920/1‑148 (MQ=255)
                                                                                                                         aTTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGATTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCg                          >  2:787228/1‑148 (MQ=255)
                                                                                                                             agttCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACg                      <  2:981071/146‑1 (MQ=255)
                                                                                                                             aTTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACg                      >  2:2219259/1‑148 (MQ=255)
                                                                                                                               ttCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTc                    <  2:1104869/148‑1 (MQ=255)
                                                                                                                               ttCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTc                    >  1:2049920/1‑148 (MQ=255)
                                                                                                                               ttCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTc                    >  1:2045753/1‑148 (MQ=255)
                                                                                                                                  cccGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTa                 <  2:534728/148‑1 (MQ=255)
                                                                                                                                   ccGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTaa                <  2:10558/148‑1 (MQ=255)
                                                                                                                                          agagATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGt         >  2:1575008/1‑148 (MQ=255)
                                                                                                                                           gagaTCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGCTGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTa        >  2:1685270/1‑148 (MQ=255)
                                                                                                                                           gagaTCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTa        >  2:1616997/1‑148 (MQ=255)
                                                                                                                                            agaTCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTaa       <  2:2444092/148‑1 (MQ=255)
                                                                                                                                            agaTCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTaa       >  1:253130/1‑148 (MQ=255)
                                                                                                                                            agaTCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTAGTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTaa       <  2:2049036/148‑1 (MQ=255)
                                                                                                                                              aTCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTAAAt     <  1:998392/148‑1 (MQ=255)
                                                                                                                                              aTCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTAAAt     <  1:1407007/148‑1 (MQ=255)
                                                                                                                                                 gATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTAAATCGa  <  2:940812/148‑1 (MQ=255)
                                                                                                                                                 gATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGCAAATCGa  <  2:2722881/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: