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,7310GT98.5% 1063.4 / 4.8 327D182Y (GAT→TAT) mnmEtRNA U34 5‑methylaminomethyl‑2‑thiouridine modification GTPase
Reads supporting (aligned to +/- strand):  ref base G (3/2);  new base T (163/159);  total (166/161)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01
Rejected as polymorphism: Frequency below/above cutoff threshold.

Only 100 of 327 total aligned reads displayed.

ATCGACGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCGATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTAAATCG  >  CP009273/3880584‑3880877
                                                                                                                                                   |                                                                                                                                                  
aTCGACGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCTAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGTAGAGATCTAt                                                                                                                                                  >  1:1841427/1‑150 (MQ=255)
aTCGACGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTAt                                                                                                                                                  <  1:2320091/150‑1 (MQ=255)
  cGACGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTAttt                                                                                                                                                <  1:3733348/150‑1 (MQ=255)
   gACGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTc                                                                                                                                               <  1:3628311/150‑1 (MQ=255)
   gACGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTc                                                                                                                                               >  1:4390129/1‑150 (MQ=255)
    aCGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTcc                                                                                                                                              >  2:2425264/1‑150 (MQ=255)
       ccAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCctct                                                                                                                                           <  1:3680865/150‑1 (MQ=255)
       ccAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGAATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCctct                                                                                                                                           <  1:3059468/150‑1 (MQ=255)
        cAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCtctc                                                                                                                                          <  2:3921021/150‑1 (MQ=255)
         aGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTcc                                                                                                                                         >  1:1371480/1‑150 (MQ=255)
         aGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTcc                                                                                                                                         <  2:781883/150‑1 (MQ=255)
         aGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTcc                                                                                                                                         >  2:1028492/1‑150 (MQ=255)
         aGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTcc                                                                                                                                         <  2:1486828/150‑1 (MQ=255)
         aGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATGTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTcc                                                                                                                                         <  2:612204/150‑1 (MQ=255)
          gTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCg                                                                                                                                        <  2:4770947/150‑1 (MQ=255)
            tCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGAc                                                                                                                                      <  1:3492991/150‑1 (MQ=255)
                aaCAGGCGGCCCGGTCGTCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACGTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGaa                                                                                                                                  >  2:3033338/1‑150 (MQ=255)
                   aGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGaaaaa                                                                                                                               <  1:5215039/150‑1 (MQ=255)
                      cggcCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTg                                                                                                                            <  2:208312/150‑1 (MQ=255)
                      cggcCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCGATTTCCTCTCCGACGGAAAAATTg                                                                                                                            <  2:2775174/150‑1 (MQ=255)
                         cccGTTCGGCACTTAACTCGCTGCTAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCTATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAg                                                                                                                         <  2:4112867/150‑1 (MQ=255)
                           cGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCGATTTCCTCTCCGACGGAAAAATTGAAGcc                                                                                                                       >  1:2949553/1‑150 (MQ=255)
                            gTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGccc                                                                                                                      <  1:3012985/150‑1 (MQ=255)
                             ttCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCa                                                                                                                     <  1:540372/150‑1 (MQ=255)
                              tCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAg                                                                                                                    <  2:3100988/150‑1 (MQ=255)
                                    cTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAAt                                                                                                              <  2:4326599/150‑1 (MQ=255)
                                    cTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAAt                                                                                                              >  1:871758/1‑150 (MQ=255)
                                       aaCTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGAc                                                                                                           <  1:4306996/150‑1 (MQ=255)
                                        aCTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACg                                                                                                          <  1:1053765/150‑1 (MQ=255)
                                         cTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGt                                                                                                         >  1:1007069/1‑150 (MQ=255)
                                         cTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGt                                                                                                         >  1:734735/1‑150 (MQ=255)
                                             cTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTAtt                                                                                                     <  1:4107055/150‑1 (MQ=255)
                                                 aaGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCg                                                                                                 >  2:2827917/1‑150 (MQ=255)
                                                 aaGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCg                                                                                                 >  2:5248429/1‑150 (MQ=255)
                                                   ggCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGAt                                                                                               <  1:5009190/150‑1 (MQ=255)
                                                    gcgcATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATc                                                                                              <  1:2671581/150‑1 (MQ=255)
                                                      gcATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCtt                                                                                            <  2:1765770/150‑1 (MQ=255)
                                                       cATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTg                                                                                           >  1:3588963/1‑150 (MQ=255)
                                                        aTTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGa                                                                                          >  2:3772405/1‑150 (MQ=255)
                                                                cACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAgtgc                                                                                  >  1:3820385/1‑150 (MQ=255)
                                                                  cGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAgtgcgt                                                                                <  1:2161761/150‑1 (MQ=255)
                                                                   gggTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAgtgcgtg                                                                               <  2:184037/150‑1 (MQ=255)
                                                                   gggTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAgtgcgtg                                                                               >  2:3620982/1‑150 (MQ=255)
                                                                      ttAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTg                                                                            <  1:1627784/150‑1 (MQ=255)
                                                                          tcatcTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGc                                                                        <  2:2249333/150‑1 (MQ=255)
                                                                            atcTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCAc                                                                      <  2:5328931/150‑1 (MQ=255)
                                                                            atcTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCAc                                                                      <  1:786469/150‑1 (MQ=255)
                                                                                   aGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAggg                                                               >  2:2197458/1‑150 (MQ=255)
                                                                                   aGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAggg                                                               >  1:1880082/1‑150 (MQ=255)
                                                                                    gAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGt                                                              <  2:3841108/150‑1 (MQ=255)
                                                                                    gAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGt                                                              <  1:2399986/150‑1 (MQ=255)
                                                                                    gAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGt                                                              <  1:5076001/150‑1 (MQ=255)
                                                                                     aaGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTa                                                             <  1:1028492/150‑1 (MQ=255)
                                                                                         ccTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGttt                                                         >  1:1788599/1‑150 (MQ=255)
                                                                                          cTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTttg                                                        >  1:225271/1‑150 (MQ=255)
                                                                                          cTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCGATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTttg                                                        >  1:1615230/1‑150 (MQ=255)
                                                                                            cACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTttgtt                                                      <  1:3326946/150‑1 (MQ=255)
                                                                                            cACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTttgtt                                                      >  1:4705192/1‑150 (MQ=255)
                                                                                             aCCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTttgttg                                                     >  1:1642527/1‑150 (MQ=255)
                                                                                               ccACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTgcg                                                   >  1:1983201/1‑150 (MQ=255)
                                                                                                   ttGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGaa                                               >  1:4156847/1‑150 (MQ=255)
                                                                                                    tGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAg                                              <  1:2216263/150‑1 (MQ=255)
                                                                                                      cgcATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAggg                                            >  2:2720308/1‑150 (MQ=255)
                                                                                                       gcATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGa                                           >  1:1723368/1‑150 (MQ=255)
                                                                                                        cATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGAt                                          <  2:2595376/150‑1 (MQ=255)
                                                                                                         aTTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATg                                         >  2:1646927/1‑150 (MQ=255)
                                                                                                         aTTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGAAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATg                                         >  1:4900013/1‑150 (MQ=255)
                                                                                                               gTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAgtg                                   >  1:2483788/1‑150 (MQ=255)
                                                                                                                 cGAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAgtggt                                 <  2:1922746/150‑1 (MQ=255)
                                                                                                                  gAAGCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAgtggtg                                >  2:3369633/1‑150 (MQ=255)
                                                                                                                     gCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGAtt                             >  2:3946387/1‑150 (MQ=255)
                                                                                                                     gCGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGAtt                             >  2:4262860/1‑150 (MQ=255)
                                                                                                                      cGGCAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTg                            <  1:2503618/150‑1 (MQ=255)
                                                                                                                         cAATTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCg                         <  1:1500069/150‑1 (MQ=255)
                                                                                                                         cAATTGATTTCCCCGATGAAGAGATCTATAACCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGATATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGACAGGGTAGTTTGTTGCGCTAAGGGATGAAAGTGGTGATTGCCg                         >  1:4670203/1‑150 (MQ=255)
                                                                                                                          aaTTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCgg                        >  1:707586/1‑150 (MQ=255)
                                                                                                                          aaTTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCgg                        >  1:288234/1‑150 (MQ=255)
                                                                                                                          aaTTGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCgg                        >  2:5066650/1‑150 (MQ=255)
                                                                                                                             tGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACg                     >  2:3546559/1‑150 (MQ=255)
                                                                                                                             tGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACg                     >  2:1323086/1‑150 (MQ=255)
                                                                                                                             tGATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGAACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACg                     >  1:322989/1‑150 (MQ=255)
                                                                                                                             tGATGTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACg                     <  2:1154762/150‑1 (MQ=255)
                                                                                                                              gATTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGATGGGATGAAAGTGGTGATTGCCGGACGt                    >  2:4373357/1‑150 (MQ=255)
                                                                                                                               aTTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTc                   <  2:4897644/150‑1 (MQ=255)
                                                                                                                               aTTTCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTc                   >  1:4818257/1‑150 (MQ=255)
                                                                                                                                  tCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTa                <  2:1620321/150‑1 (MQ=255)
                                                                                                                                  tCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTa                <  2:2106857/150‑1 (MQ=255)
                                                                                                                                  tCCCCGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTa                <  2:325109/150‑1 (MQ=255)
                                                                                                                                    cccGATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAAc              >  2:4573542/1‑150 (MQ=255)
                                                                                                                                       gATGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGcc           >  1:4547263/1‑150 (MQ=255)
                                                                                                                                        aTGAAGAGATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGACGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCg          >  2:5166016/1‑150 (MQ=255)
                                                                                                                                        aTGAAGAGATCTATGTCCTCTCCGACGGAAAAATTGACGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCg          <  1:39343/150‑1 (MQ=255)
                                                                                                                                            agagATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTaa      >  1:3143391/1‑150 (MQ=255)
                                                                                                                                            agagATCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTaa      >  2:2549775/1‑150 (MQ=255)
                                                                                                                                              agaTCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTAAAt    >  2:286180/1‑150 (MQ=255)
                                                                                                                                               gaTCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTAAATc   >  2:452363/1‑150 (MQ=255)
                                                                                                                                               gaTCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTAAATc   >  1:1818065/1‑150 (MQ=255)
                                                                                                                                               gaTCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTAAATc   >  2:2951322/1‑150 (MQ=255)
                                                                                                                                               gaTCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTAAATc   >  1:4935324/1‑150 (MQ=255)
                                                                                                                                                aTCTATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTAAATCg  >  2:2508930/1‑150 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
ATCGACGCCAGTTCGGAACAGGCGGCCCGTTCGGCACTTAACTCGCTGCAAGGCGCATTCTCCGCACGGGTTAATCATCTGGTAGAAGCCCTCACCCACTTGCGCATTTACGTCGAAGCGGCAATTGATTTCCCCGATGAAGAGATCGATTTCCTCTCCGACGGAAAAATTGAAGCCCAGCTCAATGACGTTATTGCCGATCTTGATGCAGTGCGTGCTGAAGCACGTCAGGGTAGTTTGTTGCGCGAAGGGATGAAAGTGGTGATTGCCGGACGTCCTAACGCCGGTAAATCG  >  CP009273/3880584‑3880877

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: