Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,046,600 G→T 5.7% pseudogene (210/1053 nt) yeeL ← pseudogene, glycosyltransferase homology

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,046,6000GT 5.7% 806.5 / 18.7 210pseudogene (210/1053 nt)yeeLpseudogene, glycosyltransferase homology
Reads supporting (aligned to +/- strand):  ref base G (107/91);  new base T (4/8);  total (111/99)
Fisher's exact test for biased strand distribution p-value = 2.34e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 222 total aligned reads displayed.

TTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGATGACAACAAAGTCTGTTTTTTCAATATTGAAATT  >  CP009273/2046451‑2046749
                                                                                                                                                     |                                                                                                                                                     
ttACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAg                                                                                                                                                       >  1:2814978/1‑150 (MQ=255)
ttACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAg                                                                                                                                                       >  1:2081498/1‑150 (MQ=255)
  aCGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGAc                                                                                                                                                     >  1:3102862/1‑150 (MQ=255)
      aTCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCt                                                                                                                                                 >  1:444679/1‑150 (MQ=255)
          ccTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAaat                                                                                                                                             >  2:1737515/1‑150 (MQ=255)
          ccTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAaat                                                                                                                                             >  1:3083860/1‑150 (MQ=255)
           cTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAATACTTCTAaatt                                                                                                                                            <  2:2969662/150‑1 (MQ=255)
           cTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAaatt                                                                                                                                            <  2:3169867/150‑1 (MQ=255)
           cTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCAGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAaatt                                                                                                                                            >  1:3331637/1‑150 (MQ=255)
             cAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAATACTTCTAaattaa                                                                                                                                          <  1:3615962/150‑1 (MQ=255)
               aTGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAaattaatt                                                                                                                                        <  1:2519021/150‑1 (MQ=255)
                tGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTg                                                                                                                                       <  2:842930/150‑1 (MQ=255)
                 gCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGt                                                                                                                                      <  2:1241096/150‑1 (MQ=255)
                  cAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTc                                                                                                                                     >  1:3975428/1‑150 (MQ=255)
                    aaGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGa                                                                                                                                   >  2:1401470/1‑150 (MQ=255)
                       aTCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGt                                                                                                                                <  2:1376992/150‑1 (MQ=255)
                          cGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTAc                                                                                                                             >  1:1893942/1‑150 (MQ=255)
                                aCCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAATACTTCTAAATTAATTGTCGACGTTACAATTaa                                                                                                                       <  2:3437418/150‑1 (MQ=255)
                                  cGTATTGCTGTGTTATGTTCTTTGCTAAATNTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTaaaa                                                                                                                     >  1:4279609/1‑150 (MQ=255)
                                  cGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTaaaa                                                                                                                     >  1:3173883/1‑150 (MQ=255)
                                     aTTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGAt                                                                                                                  >  2:1585527/1‑150 (MQ=255)
                                       tGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGc                                                                                                                <  1:1950801/150‑1 (MQ=255)
                                         cTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCtt                                                                                                              >  1:2807309/1‑150 (MQ=255)
                                         cTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCtt                                                                                                              <  2:2337635/150‑1 (MQ=255)
                                          tgtgTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTa                                                                                                             <  2:4519199/150‑1 (MQ=255)
                                            tgtTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAAt                                                                                                           <  2:1267281/150‑1 (MQ=255)
                                             gtTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAAtt                                                                                                          <  2:3720932/150‑1 (MQ=255)
                                             gtTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAAtt                                                                                                          <  2:3720647/150‑1 (MQ=255)
                                                aTGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTa                                                                                                       <  1:4174986/150‑1 (MQ=255)
                                                 tGTTCTTTGCTAAATGTGTTATTTCTTTTAAGGTTGCGTGATCCTGCATAGTCATGGTGACTAGTGTTAATTTAGTCTGTGCAAGGTTAAGCGCATTAATGACTTCTCAATTACTTGTCGACGTTACGATTACAAGATGCTTAGTTTTAt                                                                                                      >  2:3537868/1‑150 (MQ=255)
                                                 tGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTAt                                                                                                      >  2:3537904/1‑150 (MQ=255)
                                                  gTTCTTTGCTAAATCTTTTATTACTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATg                                                                                                     <  2:1555202/150‑1 (MQ=255)
                                                    tCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCa                                                                                                   >  2:582895/1‑150 (MQ=255)
                                                      tttGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAAt                                                                                                 <  2:2012726/150‑1 (MQ=255)
                                                       ttGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAAtt                                                                                                >  2:4242206/1‑150 (MQ=255)
                                                         gCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCa                                                                                              >  2:3927458/1‑150 (MQ=255)
                                                          cTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCaa                                                                                             >  1:4356384/1‑150 (MQ=255)
                                                            aaaTCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAgc                                                                                           >  1:1716510/1‑150 (MQ=255)
                                                            aaaTCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAgc                                                                                           >  1:1805150/1‑150 (MQ=255)
                                                            aaaTCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAgc                                                                                           <  1:3550670/150‑1 (MQ=255)
                                                            aaaTCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAgc                                                                                           <  1:2588347/150‑1 (MQ=255)
                                                                 ttttATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCg                                                                                      <  1:2177118/150‑1 (MQ=255)
                                                                  tttATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGa                                                                                     >  2:1449559/1‑150 (MQ=255)
                                                                   ttATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGaa                                                                                    <  2:4183051/150‑1 (MQ=255)
                                                                     aTTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAATACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATa                                                                                  >  2:2260980/1‑150 (MQ=255)
                                                                     aTTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATa                                                                                  >  1:271622/1‑150 (MQ=255)
                                                                          ttttAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAgg                                                                             <  2:3007538/150‑1 (MQ=255)
                                                                           tttAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGa                                                                            >  1:1087060/1‑150 (MQ=255)
                                                                           tttAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGa                                                                            >  2:1412567/1‑150 (MQ=255)
                                                                           tttAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGa                                                                            >  1:1087074/1‑150 (MQ=255)
                                                                               aTGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGa                                                                        >  1:249656/1‑150 (MQ=255)
                                                                                tGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGAt                                                                       >  2:1811712/1‑150 (MQ=255)
                                                                                 gTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAATACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATg                                                                      >  2:1849196/1‑150 (MQ=255)
                                                                                 gTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATg                                                                      >  2:4191050/1‑150 (MQ=255)
                                                                                  ttGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATgg                                                                     <  1:1034820/150‑1 (MQ=255)
                                                                                   tGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGc                                                                    >  2:264865/1‑150 (MQ=255)
                                                                                     cGTGATCCTGCATAGTCTTTGTGTCTAATGTTAATTTATTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGAAGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAATCGCCCTAATAACAGGAAAGATGGCCa                                                                  <  2:2705333/150‑1 (MQ=255)
                                                                                       tGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATa                                                                >  1:4231214/1‑150 (MQ=255)
                                                                                       tGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATa                                                                >  2:2538876/1‑150 (MQ=255)
                                                                                        gATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAg                                                               <  2:1432706/150‑1 (MQ=255)
                                                                                        gATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAg                                                               >  2:582388/1‑150 (MQ=255)
                                                                                        gATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAg                                                               >  2:2884844/1‑150 (MQ=255)
                                                                                         aTCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGc                                                              >  2:469725/1‑150 (MQ=255)
                                                                                         aTCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGc                                                              >  2:3356622/1‑150 (MQ=255)
                                                                                          tCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCa                                                             >  1:1504408/1‑150 (MQ=255)
                                                                                           ccTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCAt                                                            >  1:4060142/1‑150 (MQ=255)
                                                                                            cTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATc                                                           <  1:1161935/150‑1 (MQ=255)
                                                                                               cATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGcc                                                        <  1:3840268/150‑1 (MQ=255)
                                                                                                aTAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCa                                                       >  1:3105746/1‑150 (MQ=255)
                                                                                                      aTTGTGACTAATGTTAATTTAGTCTGTTCAAGTTAAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGa                                                 >  2:3204129/1‑150 (MQ=255)
                                                                                                         gtgACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCg                                              <  2:2157679/150‑1 (MQ=255)
                                                                                                              tAATGTTAATTTAGTCTGTTCAAGTTTATGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAat                                         >  1:879202/1‑150 (MQ=255)
                                                                                                              tAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAat                                         <  2:4260218/150‑1 (MQ=255)
                                                                                                                 tGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGGTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATg                                      <  1:2858887/150‑1 (MQ=255)
                                                                                                                    tAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGAt                                   <  1:3693043/150‑1 (MQ=255)
                                                                                                                    tAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGAt                                   >  2:3774176/1‑150 (MQ=255)
                                                                                                                        ttAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGTATGaca                               <  2:508900/150‑1 (MQ=255)
                                                                                                                        ttAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGATGaca                               >  1:3045204/1‑150 (MQ=255)
                                                                                                                         tAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGATGacaa                              >  2:4065984/1‑150 (MQ=255)
                                                                                                                           gTCTGTTCAAGTTTAAGCGCATTAAAGAATTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGATGacaaca                            <  1:3949913/150‑1 (MQ=255)
                                                                                                                            tCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGATGAcaacaa                           >  1:2408335/1‑150 (MQ=255)
                                                                                                                            tCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGATGAcaacaa                           >  2:1174427/1‑150 (MQ=255)
                                                                                                                                 tCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGATGACAACAAAGTCt                      >  2:4370951/1‑150 (MQ=255)
                                                                                                                                  cAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGATGACAACAAAGTCTg                     <  2:739408/150‑1 (MQ=255)
                                                                                                                                  cAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGATGACAACAAAGTCTg                     >  1:601098/1‑150 (MQ=255)
                                                                                                                                   aaGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGATGACAACAAAGTCTGt                    <  1:1360104/150‑1 (MQ=255)
                                                                                                                                      tttAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGATGACAACAAAGTCTGtttt                 <  2:109338/150‑1 (MQ=255)
                                                                                                                                       ttAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGATGACAACAAAGTCTGttttt                <  1:580647/150‑1 (MQ=255)
                                                                                                                                       ttAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGATGACAACAAAGTCTGttttt                <  2:4268906/150‑1 (MQ=255)
                                                                                                                                         aaGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGATGACAACAAAGTCTGTTTTTTc              >  2:2741143/1‑150 (MQ=255)
                                                                                                                                          aGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGATGACAACAAAGTCTGTTTTTTCa             <  1:692905/150‑1 (MQ=255)
                                                                                                                                            cgcATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGATGACAACAAAGTCTGTTTTTTCAat           >  2:2647173/1‑150 (MQ=255)
                                                                                                                                              cATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGATGACAACAAAGTCTGTTTTTTCAatat         <  2:1392674/150‑1 (MQ=255)
                                                                                                                                              cATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGATGACAACAAAGTCTGTTTTTTCAatat         <  2:3083977/150‑1 (MQ=255)
                                                                                                                                              cATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGATGACAACAAAGTCTGTTTTTTCAatat         <  2:3500388/150‑1 (MQ=255)
                                                                                                                                                ttAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGATGACAACAAAGTCTGTTTTTTCAATATTg       <  2:37632/150‑1 (MQ=255)
                                                                                                                                                 tAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGATGACAACAAAGTCTGTTTTTTCAATATTGa      <  1:3920717/150‑1 (MQ=255)
                                                                                                                                                   aaGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGATGACAACAAAGTCTGTTTTTTCAATATTGaaa    >  2:392110/1‑150 (MQ=255)
                                                                                                                                                   aaGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGATGACAACAAAGTCTGTTTTTTCAATATTGaaa    <  1:2002482/150‑1 (MQ=255)
                                                                                                                                                     gACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGATGACAACAAAGTCTGTTTTTTCAATATTGAAAtt  <  1:1263092/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TTACGCATCCCCTCAATGCAAAGATCCGGCGTACCGTATTGCTGTGTTATGTTCTTTGCTAAATCTTTTATTTCTTTTAATGTTGCGTGATCCTGCATAGTCATTGTGACTAATGTTAATTTAGTCTGTTCAAGTTTAAGCGCATTAAAGACTTCTAAATTAATTGTCGACGTTACAATTAAAAGATGCTTAATTTTATGCAATTCAAGCGCCCGAATAACAGGAAAGATGGCCATAGCATCGCCAATCTGATCGGGAATATGGATGACAACAAAGTCTGTTTTTTCAATATTGAAATT  >  CP009273/2046451‑2046749

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: