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.6% 1099.9 / 18.9 353D182Y (GAT→TAT) mnmEtRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase
Reads supporting (aligned to +/- strand):  ref base G (8/5);  new base T (172/168);  total (180/173)
Fisher's exact test for biased strand distribution p-value = 5.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 354 total aligned reads displayed.

TTATCGACGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCGATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTAAATC  >  CP009273/3880582‑3880876
                                                                                                                                                     |                                                                                                                                                 
  aTCGACGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTAt                                                                                                                                                 <  2:3315983/150‑1 (MQ=255)
   tCGACGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTAtt                                                                                                                                                <  2:3277180/150‑1 (MQ=255)
       cGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCct                                                                                                                                            >  2:1002027/1‑150 (MQ=255)
        gCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCctc                                                                                                                                           <  2:2644132/150‑1 (MQ=255)
        gCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTAGTTCctc                                                                                                                                           >  1:4754641/1‑150 (MQ=255)
         ccAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCctct                                                                                                                                          >  2:883800/1‑150 (MQ=255)
          cAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCtctc                                                                                                                                         >  2:4610456/1‑150 (MQ=255)
           aGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTcc                                                                                                                                        >  2:2071554/1‑150 (MQ=255)
              tCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGAc                                                                                                                                     <  2:1268021/150‑1 (MQ=255)
              tCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGAc                                                                                                                                     <  1:4631971/150‑1 (MQ=255)
               cGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACg                                                                                                                                    >  1:424913/1‑150 (MQ=255)
               cGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACg                                                                                                                                    >  2:172563/1‑150 (MQ=255)
                 gAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGa                                                                                                                                  >  1:509070/1‑150 (MQ=255)
                  aaCAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGaa                                                                                                                                 >  1:2062470/1‑150 (MQ=255)
                  aaCAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGaa                                                                                                                                 >  1:1881127/1‑150 (MQ=255)
                    cAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGaaaa                                                                                                                               >  2:2951371/1‑150 (MQ=255)
                     aGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGaaaaa                                                                                                                              >  2:1031706/1‑150 (MQ=255)
                      ggcggcCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCAGTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAAt                                                                                                                             <  2:3611740/150‑1 (MQ=255)
                           cccGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAg                                                                                                                        >  2:4110097/1‑150 (MQ=255)
                           cccGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAg                                                                                                                        >  2:3775968/1‑150 (MQ=255)
                           cccGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAg                                                                                                                        >  1:1260636/1‑150 (MQ=255)
                            ccGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGc                                                                                                                       <  2:645570/150‑1 (MQ=255)
                               ttCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCa                                                                                                                    >  2:4425426/1‑150 (MQ=255)
                               ttCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCa                                                                                                                    >  1:352249/1‑150 (MQ=255)
                               ttCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCa                                                                                                                    <  2:33488/150‑1 (MQ=255)
                                tCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAg                                                                                                                   >  1:5163441/1‑150 (MQ=255)
                                  ggCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCt                                                                                                                 <  2:4262876/150‑1 (MQ=255)
                                    cACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCa                                                                                                               <  1:3106479/150‑1 (MQ=255)
                                        tAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGa                                                                                                           >  2:4733613/1‑150 (MQ=255)
                                          aCTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACg                                                                                                         <  1:1184854/150‑1 (MQ=255)
                                          aCTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACg                                                                                                         <  1:4085040/150‑1 (MQ=255)
                                             tgCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTa                                                                                                      >  1:4198897/2‑150 (MQ=255)
                                              gCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTAt                                                                                                     >  2:3859323/1‑150 (MQ=255)
                                               cTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTAtt                                                                                                    <  2:5041877/150‑1 (MQ=255)
                                               cTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTAtt                                                                                                    <  2:276604/150‑1 (MQ=255)
                                                  cAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGcc                                                                                                 <  2:4810050/150‑1 (MQ=255)
                                                             ctcCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGc                                                                                      >  1:3790507/1‑150 (MQ=255)
                                                              tcCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCa                                                                                     >  1:1980524/1‑150 (MQ=255)
                                                               ccGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTGGCGCATTGACGTCGAAGCGGCAATTGATTGCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAg                                                                                    <  2:2614709/150‑1 (MQ=255)
                                                                cGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCGATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAgt                                                                                   >  2:2500937/1‑150 (MQ=255)
                                                                  cACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAgtgc                                                                                 >  1:1309562/1‑150 (MQ=255)
                                                                  cACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAgtgc                                                                                 >  1:1826520/1‑150 (MQ=255)
                                                                   aCGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGAAGTTATTGCCGATCTTGATGCAgtgcg                                                                                >  2:57239/1‑150 (MQ=255)
                                                                       gTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCt                                                                            <  1:4895872/150‑1 (MQ=255)
                                                                        ttAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCTATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTg                                                                           <  1:720158/150‑1 (MQ=255)
                                                                          aaTCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGaa                                                                         <  2:4528849/150‑1 (MQ=255)
                                                                            tcatcTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGc                                                                       >  2:810275/1‑150 (MQ=255)
                                                                               tcTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACg                                                                    <  1:3236261/150‑1 (MQ=255)
                                                                                 tGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTc                                                                  <  2:1207852/150‑1 (MQ=255)
                                                                                  ggTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCa                                                                 <  1:3343967/150‑1 (MQ=255)
                                                                                  ggTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCa                                                                 <  1:4431993/150‑1 (MQ=255)
                                                                                    tAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAgg                                                               <  2:1432599/150‑1 (MQ=255)
                                                                                        aGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAg                                                           >  1:2080094/1‑150 (MQ=255)
                                                                                          cccTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGtt                                                         <  2:1355920/150‑1 (MQ=255)
                                                                                              cACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGGAGTttgtt                                                     >  2:4359952/1‑150 (MQ=255)
                                                                                               aCCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTttgttg                                                    >  2:495185/1‑150 (MQ=255)
                                                                                               aCCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTttgttg                                                    <  2:1780005/150‑1 (MQ=255)
                                                                                               aCCCACTTGCGCATGTACGTCGAAGCGGCAATTGATGTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTttgttg                                                    <  1:2198496/150‑1 (MQ=255)
                                                                                                cccACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTgc                                                   >  1:777119/1‑150 (MQ=255)
                                                                                                 ccACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTgcg                                                  <  1:1002027/150‑1 (MQ=255)
                                                                                                   aCTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGcgcg                                                >  1:2878904/1‑150 (MQ=255)
                                                                                                   aCTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGcgcg                                                <  2:1621407/150‑1 (MQ=255)
                                                                                                      tGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAg                                             >  2:1993453/1‑150 (MQ=255)
                                                                                                      tGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAg                                             <  2:1927990/150‑1 (MQ=255)
                                                                                                        cgcATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAggg                                           <  2:2744756/150‑1 (MQ=255)
                                                                                                         gcATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGa                                          <  2:509134/150‑1 (MQ=255)
                                                                                                         gcATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGa                                          <  2:1124068/150‑1 (MQ=255)
                                                                                                           aTTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATg                                        >  2:2201161/1‑150 (MQ=255)
                                                                                                           aTTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATg                                        <  1:4481356/150‑1 (MQ=255)
                                                                                                           aTTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATg                                        <  2:521268/150‑1 (MQ=255)
                                                                                                              tACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGGTGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGaaa                                     >  1:3065017/1‑150 (MQ=255)
                                                                                                              tACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGaaa                                     >  2:4478395/1‑150 (MQ=255)
                                                                                                              tACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTGTGTTGCGCGAAGGGATGaaa                                     >  2:4702091/1‑150 (MQ=255)
                                                                                                                 gTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAgtg                                  <  2:2537109/150‑1 (MQ=255)
                                                                                                                  tCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAgtgg                                 >  2:2031853/1‑150 (MQ=255)
                                                                                                                  tCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAgtgg                                 >  2:2405617/1‑150 (MQ=255)
                                                                                                                  tCGAAGCGGCAAGTGATTTCCCCGATGAAGAGATCTAGTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAgtgg                                 <  2:2277847/150‑1 (MQ=255)
                                                                                                                    gAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAgtggtg                               <  1:2951371/150‑1 (MQ=255)
                                                                                                                      aGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGAt                             <  2:1881185/150‑1 (MQ=255)
                                                                                                                       gCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAACAATTGAAGCCCAGCTCAAGGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGGTGCGCGAAGGGATGAAAGTGGTGAtt                            <  2:3219458/150‑1 (MQ=255)
                                                                                                                       gCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGAtt                            >  1:4550731/1‑150 (MQ=255)
                                                                                                                       gCGGCAATTGATTTCCCCGATGAAGAGATCGATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGAtt                            >  1:2444507/1‑150 (MQ=255)
                                                                                                                        cGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTg                           >  1:2561329/1‑150 (MQ=255)
                                                                                                                        cGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTg                           >  1:1747804/1‑150 (MQ=255)
                                                                                                                          gCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGcc                         <  2:899763/150‑1 (MQ=255)
                                                                                                                             aTTGATTTCCCCGATGAAGAGATCGATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGa                      >  1:2644739/1‑150 (MQ=255)
                                                                                                                                gATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGt                   >  2:3537578/1‑150 (MQ=255)
                                                                                                                                   ttCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCt                >  1:2296183/1‑150 (MQ=255)
                                                                                                                                    tCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTa               >  1:1085265/1‑150 (MQ=255)
                                                                                                                                     ccccGATGAAGAGATCGATTTCCTCTCCGACGGAAAAATAGAAGCCCAGCTCAATGACGTTATAGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGTGATGAAAGAGGTTATTGACGGACGTCCTaa              >  1:2328218/1‑150 (MQ=255)
                                                                                                                                      cccGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCCATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGACAGTGGTGATTGCCGGACGTCCTAAc             >  2:4823188/1‑150 (MQ=255)
                                                                                                                                      cccGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAAc             <  2:2365967/150‑1 (MQ=255)
                                                                                                                                         gATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGcc          <  1:4673610/150‑1 (MQ=255)
                                                                                                                                         gATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGcc          <  1:2564610/150‑1 (MQ=255)
                                                                                                                                           tGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCgg        >  1:642880/1‑150 (MQ=255)
                                                                                                                                           tGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCgg        >  1:3455943/1‑150 (MQ=255)
                                                                                                                                            gAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGt       <  1:146820/150‑1 (MQ=255)
                                                                                                                                             aaGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTa      <  1:3398151/150‑1 (MQ=255)
                                                                                                                                               gagaTCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTaaa    <  1:4011304/150‑1 (MQ=255)
                                                                                                                                                 gaTCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTAAATc  >  2:3153233/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
TTATCGACGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCGATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTAAATC  >  CP009273/3880582‑3880876

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: ‑