Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP006881 1,666,671 T→A 18.3% L321* (TTA→TAA)  U712_08575 → Flagellar motor switch phosphatase FliY

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0068811,666,6710TA18.3% 927.7 / 188.8 366L321* (TTA→TAA) U712_08575Flagellar motor switch phosphatase FliY
Reads supporting (aligned to +/- strand):  ref base T (166/133);  new base A (36/31);  total (202/164)
Fisher's exact test for biased strand distribution p-value = 7.88e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

Only 100 of 375 total aligned reads displayed.

GCCAGTTCAGGTTTCACCGGTGGAATTTTCCGCGTTTGATCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTCGTTGTCATTGAAGAAAACTTCGGAGTGAGAGTCACTGACATTTTAAGCCAGG  >  CP006881/1666522‑1666820
                                                                                                                                                     |                                                                                                                                                     
  cAGTTCAGGTTTCACCGGTGGAATTTTCCGCGTTTGATCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTag                                                                                                                                                     >  1:2010176/1‑150 (MQ=255)
  cAGTTCAAGTTTCACCGGTGGAATTTTCCGCGTTTGATCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTag                                                                                                                                                     >  2:2056452/1‑150 (MQ=255)
    gTTCAGGTTTCACCGGTGGAATTTTCCGCGTTTGATCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTagag                                                                                                                                                   >  2:455831/1‑150 (MQ=255)
     ttCAGGTTTCACCGGTGGAATTTTCCGCGTTTGATCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCNGCTTTCGATNACTGTTGAGCTTGGCAGAACGAANAGAAGCGTAAAAGAGATTTTngagn                                                                                                                                                  <  1:4520381/150‑2 (MQ=255)
         ggTTTCNCCGGTGGAATTTTCCGCGTTTGATCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTc                                                                                                                                              <  1:3699093/150‑1 (MQ=255)
            ttCACCGGTGGAATTTTCCGCGTTTGATCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGc                                                                                                                                           >  2:2182062/1‑150 (MQ=255)
             tCACCGGTGGAATTTTCCGCGTTTGATCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGcc                                                                                                                                          >  2:809575/1‑150 (MQ=255)
               aCCGGTGGAATTTTCCGCGTTTGATCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCgg                                                                                                                                        <  1:209643/150‑1 (MQ=255)
               aCCGGTGGAATTTTCCGCGTTTGATCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCgg                                                                                                                                        <  2:1558439/150‑1 (MQ=255)
                 cGGTGGAATTTTCCGCGTTTGATCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGaa                                                                                                                                      >  2:5783816/1‑150 (MQ=255)
                 cGGTGGAATTTTCCGCGTTTGATCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGaa                                                                                                                                      <  2:3265399/150‑1 (MQ=255)
                 cGGTGGAATTTTCCGCGTTTGATCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTAGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGGGATTTTAGAGCTTTCTGCCGGaa                                                                                                                                      >  2:2412359/1‑150 (MQ=255)
                  ggTGGAATTTTCCGCGTTTGATCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTAAGAGCTTTCTGCCGGAAg                                                                                                                                     >  2:3146469/1‑150 (MQ=255)
                        aTTTTCCGCGTTTGATCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGtattat                                                                                                                               <  2:3854643/150‑1 (MQ=255)
                          tttCCGCGTTTTATACAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTAAGAGCTTTCTGCCGGAAGTATTATCg                                                                                                                             <  1:4339283/150‑1 (MQ=255)
                          tttCCGCGTTTGATCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCg                                                                                                                             <  2:24692/150‑1 (MQ=255)
                             ccGCGTTTGATCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGc                                                                                                                          >  1:1438693/1‑150 (MQ=255)
                                cgTTTGATCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTg                                                                                                                       <  2:297857/150‑1 (MQ=255)
                                cgTTTGATCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTg                                                                                                                       <  1:1602289/150‑1 (MQ=255)
                                 gTTTGATCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGa                                                                                                                      >  2:4202299/1‑150 (MQ=255)
                                 gTTTGATCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGa                                                                                                                      >  2:4924781/1‑150 (MQ=255)
                                    tGATCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATaa                                                                                                                   <  1:1842632/150‑1 (MQ=255)
                                    tGATCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATaa                                                                                                                   <  1:3628959/150‑1 (MQ=255)
                                      aTCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAAt                                                                                                                 <  2:4413072/150‑1 (MQ=255)
                                       tCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAAtt                                                                                                                <  1:1925241/150‑1 (MQ=255)
                                       tCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAAtt                                                                                                                <  2:765947/150‑1 (MQ=255)
                                       tCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAAtt                                                                                                                <  1:2419948/150‑1 (MQ=255)
                                        ccAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTAAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAGATTa                                                                                                               <  2:266381/150‑1 (MQ=255)
                                           aaCGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGcc                                                                                                            <  1:4171922/150‑1 (MQ=255)
                                           aaCGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGcc                                                                                                            >  2:5578011/1‑150 (MQ=255)
                                           aaCGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTAAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGcc                                                                                                            >  2:5037426/1‑150 (MQ=255)
                                             cGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCgg                                                                                                          >  2:3341562/1‑150 (MQ=255)
                                             cGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTAAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCgg                                                                                                          <  1:1301615/150‑1 (MQ=255)
                                              gAAGCTGTACAAGCCCCTATTCATAATCTGGGTATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGt                                                                                                         >  2:1828912/1‑150 (MQ=255)
                                              gAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGt                                                                                                         <  1:5526350/150‑1 (MQ=255)
                                              gAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGt                                                                                                         >  2:551328/1‑150 (MQ=255)
                                              gAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGt                                                                                                         <  2:4284666/150‑1 (MQ=255)
                                                    gTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGGATTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCt                                                                                                   >  1:4085112/1‑150 (MQ=255)
                                                    gTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCt                                                                                                   <  1:4911387/150‑1 (MQ=255)
                                                    gTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCt                                                                                                   <  1:4740875/150‑1 (MQ=255)
                                                    gTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTAAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCt                                                                                                   <  1:2655901/150‑1 (MQ=255)
                                                     tACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTg                                                                                                  >  2:5163832/1‑150 (MQ=255)
                                                     tACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTAAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTg                                                                                                  >  1:127793/1‑150 (MQ=255)
                                                      aCAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAACAGATTTAAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGt                                                                                                 >  1:4589740/1‑150 (MQ=255)
                                                       cAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTa                                                                                                <  1:671287/150‑1 (MQ=255)
                                                       cAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTa                                                                                                <  2:512693/150‑1 (MQ=255)
                                                        aaGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAg                                                                                               <  1:5407306/150‑1 (MQ=255)
                                                         aGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTAAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGa                                                                                              >  2:5308002/1‑150 (MQ=255)
                                                           ccccTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGata                                                                                            <  1:2379323/150‑1 (MQ=255)
                                                           ccccTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGata                                                                                            <  1:2615261/150‑1 (MQ=255)
                                                             ccTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATAtt                                                                                          <  2:3649733/150‑1 (MQ=255)
                                                                aTTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTa                                                                                       >  2:1349199/1‑150 (MQ=255)
                                                                   cATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGtt                                                                                    >  2:1053187/1‑150 (MQ=255)
                                                                   cATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTAAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGtt                                                                                    >  1:2805069/1‑150 (MQ=255)
                                                                    aTAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTa                                                                                   >  2:5377508/1‑150 (MQ=255)
                                                                     tAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTaa                                                                                  >  2:1455550/1‑150 (MQ=255)
                                                                      aaTCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATATATTAGCCGGTGAGCCTGTAGATATTTTAGTTAAt                                                                                 >  1:995057/1‑150 (MQ=255)
                                                                        tCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCa                                                                               <  1:4654459/150‑1 (MQ=255)
                                                                        tCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCa                                                                               >  2:1217687/1‑150 (MQ=255)
                                                                          tGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGc                                                                             >  1:5252568/1‑150 (MQ=255)
                                                                           ggATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCg                                                                            >  2:5214595/1‑150 (MQ=255)
                                                                           ggATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCg                                                                            >  1:5676567/1‑150 (MQ=255)
                                                                                tGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCg                                                                       <  1:5163832/150‑1 (MQ=255)
                                                                                tGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCg                                                                       <  1:2187668/150‑1 (MQ=255)
                                                                                 gCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGt                                                                      >  1:5054780/1‑150 (MQ=255)
                                                                                    tttGGATATCCCTCTTACGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTAAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGc                                                                   <  1:851847/150‑1 (MQ=255)
                                                                                        gATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCaaa                                                               >  1:296513/1‑150 (MQ=255)
                                                                                         atatCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAg                                                              <  1:4005107/150‑1 (MQ=255)
                                                                                          tatCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAgg                                                             >  2:403881/1‑150 (MQ=255)
                                                                                            tCCCGCTTTCGATTACTGTTGAGCTTGACAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTACTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGAGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCg                                                           <  2:3032038/150‑1 (MQ=255)
                                                                                             cccGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTAAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGa                                                          >  2:5645359/1‑150 (MQ=255)
                                                                                               cGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAgg                                                        <  2:3543716/150‑1 (MQ=255)
                                                                                               cGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCCGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAAGAGAgg                                                        >  1:2882069/1‑150 (MQ=255)
                                                                                                 cTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTc                                                      <  1:3209725/150‑1 (MQ=255)
                                                                                                 cTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTAAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTc                                                      <  1:1126865/150‑1 (MQ=255)
                                                                                                  tttCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTCg                                                     >  2:5273537/1‑150 (MQ=255)
                                                                                                    tCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTAGCCAAAGGCGAGGTCGtt                                                   >  2:958703/1‑150 (MQ=255)
                                                                                                      gATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTCGTTGt                                                 >  1:3864301/1‑150 (MQ=255)
                                                                                                         tACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTCGTTGTCAt                                              <  2:3727482/150‑1 (MQ=255)
                                                                                                         tACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTCGTTGTCAt                                              <  2:2010176/150‑1 (MQ=255)
                                                                                                          aCTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTCGTTGTCAtt                                             >  2:1759651/1‑150 (MQ=255)
                                                                                                           cTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTCGTTGTCATTg                                            >  1:1654507/1‑150 (MQ=255)
                                                                                                            tGTTGAGCTTGGCAGGACGAAGAGTGGCGTAAAAGGGATTTAAGTGCTTTCTGCCGGAAGCATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTCGTTGTCATTga                                           >  2:616855/1‑150 (MQ=255)
                                                                                                             gTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGAATCGTTGCCAAAGGCGAGGTCGTTGTCATTgaa                                          <  1:3038066/150‑1 (MQ=255)
                                                                                                              ttGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTCGTTGTCATTgaag                                         <  2:4639288/150‑1 (MQ=255)
                                                                                                                gAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTCGTTGTCATTgaagaa                                       >  1:1398568/1‑150 (MQ=255)
                                                                                                                 aGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTCGTTGTCATTGAAGaaa                                      >  1:579777/1‑150 (MQ=255)
                                                                                                                   cTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTCGTTGTCATTGAAGAAAAc                                    >  1:175251/1‑150 (MQ=255)
                                                                                                                         aGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTCGTTGTCATTGAAGAAAACTTCGGa                              >  2:3948588/1‑150 (MQ=255)
                                                                                                                          gAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTCGTTGTCATTGAAGAAAACTTCGGAg                             <  1:3341562/150‑1 (MQ=255)
                                                                                                                                 gagaAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTCGTTGTCATTGAAGAAAACTTCGGAGTGAGAGt                      <  1:1015393/150‑1 (MQ=255)
                                                                                                                                       cGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTCGTTGTCATTGAAGAAAACTTCGGAGTGAGAGTCACTGa                >  2:4613634/1‑150 (MQ=255)
                                                                                                                                       cGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTCGTTGTCATTGAAGAAAACTTCGGAGTGAGAGTCACTGa                >  2:4828778/1‑150 (MQ=255)
                                                                                                                                           aaaGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTCGTTGTCATTGAAGAAAACTTCGGAGTGAGAGTCACTGACAtt            >  2:845362/1‑150 (MQ=255)
                                                                                                                                           aaaGAGATTTAAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTCGTTGTCATTGAAGAAAACTTCGGAGTGAGAGTCACTGACAtt            <  2:3449654/150‑1 (MQ=255)
                                                                                                                                             agagATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTCGTTGTCATTGAAGAAAACTTCGGAGTGAGAGTCACTGACAtttt          >  2:185947/1‑150 (MQ=255)
                                                                                                                                              gagaTTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTCGTTGTCATTGAAGAAAACTTCGGAGTGAGAGTCACTGACATTTTa         <  2:262258/150‑1 (MQ=255)
                                                                                                                                               agaTTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTCGTTGTCATTGAAGAAAACTTCGGAGTGAGAGTCACTGACATTTTaa        <  1:4238914/150‑1 (MQ=255)
                                                                                                                                                 aTTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTCGTTGTCATTGAAGAAAACTTCGGAGTGAGAGTCACTGACATTTTAAGc      >  1:5322928/1‑150 (MQ=255)
                                                                                                                                                     tAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTCGTTGTCATTGAAGAAAACTTCGGAGTGAGAGTCACTGACATTTTAAGCCAgg  <  2:1642752/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCCAGTTCAGGTTTCACCGGTGGAATTTTCCGCGTTTGATCCAAACGAAGCTGTACAAGCCCCTATTCATAATCTGGATATGCTTTTGGATATCCCGCTTTCGATTACTGTTGAGCTTGGCAGAACGAAGAGAAGCGTAAAAGAGATTTTAGAGCTTTCTGCCGGAAGTATTATCGAGCTTGATAAATTAGCCGGTGAGCCTGTAGATATTTTAGTTAATCAGCGTATCGTTGCCAAAGGCGAGGTCGTTGTCATTGAAGAAAACTTCGGAGTGAGAGTCACTGACATTTTAAGCCAGG  >  CP006881/1666522‑1666820

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