Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 4,520,215 T→A 14.8% N46Y (AAC→TAC)  sgcB ← putative enzyme IIB component of PTS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092734,520,2150TA14.8% 874.5 / 119.7 294N46Y (AAC→TAC) sgcBputative enzyme IIB component of PTS
Reads supporting (aligned to +/- strand):  ref base T (111/139);  new base A (28/16);  total (139/155)
Fisher's exact test for biased strand distribution p-value = 2.18e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 299 total aligned reads displayed.

CTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCCACAAGGATCTTTTTCATAATGGATGTGCCTC  >  CP009273/4520066‑4520364
                                                                                                                                                     |                                                                                                                                                     
ctcctTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGtt                                                                                                                                                       >  1:635044/1‑150 (MQ=255)
ctcctTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGtt                                                                                                                                                       <  1:51053/150‑1 (MQ=255)
  cctTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTaa                                                                                                                                                     <  1:1079247/150‑1 (MQ=255)
    ttttATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCTTAGGGATGCCGTAATCACTGTTGGTTCGCATTGAGGTTACAATAAGGTCCATACCGTTACAGTTAaga                                                                                                                                                   >  2:1028775/1‑150 (MQ=255)
    ttttATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAaga                                                                                                                                                   >  2:941714/1‑150 (MQ=255)
      ttATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAgg                                                                                                                                                 >  2:2562692/1‑150 (MQ=255)
       tATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAggg                                                                                                                                                <  2:137601/150‑1 (MQ=255)
       tATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAggg                                                                                                                                                >  2:3935112/1‑150 (MQ=255)
        aTTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTAAAGAGGGa                                                                                                                                               >  2:2252950/1‑150 (MQ=255)
         ttGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGAt                                                                                                                                              <  2:3762398/150‑1 (MQ=255)
          tGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGTAtc                                                                                                                                             <  2:1410280/150‑1 (MQ=255)
            aGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTAAAGAGGGAtctc                                                                                                                                           >  2:4501336/1‑150 (MQ=255)
                 aCAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCa                                                                                                                                      <  2:230117/150‑1 (MQ=255)
                  cAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAg                                                                                                                                     >  1:2203391/1‑150 (MQ=255)
                   agagCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGa                                                                                                                                    >  1:530611/1‑150 (MQ=255)
                   agagCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAAACACTGTTGGTTCGCATGGAGGTTACAAAAAGGTCCATACCGTTACAGTAAAGAGGGATCTCATTCAGa                                                                                                                                    >  2:1625772/1‑150 (MQ=255)
                    gagCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGAc                                                                                                                                   <  1:385868/150‑1 (MQ=255)
                     agCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTAAAGAGGGATCTCATTCAGACa                                                                                                                                  <  2:2513218/150‑1 (MQ=255)
                      gCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAg                                                                                                                                 <  1:746545/150‑1 (MQ=255)
                       ccTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGc                                                                                                                                <  1:3249790/150‑1 (MQ=255)
                        cTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCa                                                                                                                               >  1:1903501/1‑150 (MQ=255)
                        cTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTAAAGAGGGATCTCATTCAGACAGCa                                                                                                                               >  2:321105/1‑150 (MQ=255)
                         ttGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCAt                                                                                                                              <  2:3173711/150‑1 (MQ=255)
                               tgttgtTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCg                                                                                                                        >  2:3750400/1‑150 (MQ=255)
                                gttgttTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCgg                                                                                                                       >  1:508453/1‑150 (MQ=255)
                                    ttttAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGt                                                                                                                   <  1:1828458/150‑1 (MQ=255)
                                          cgcgTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGa                                                                                                             >  2:4317565/1‑150 (MQ=255)
                                             gtcgtcgtTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAAt                                                                                                          <  2:962932/150‑1 (MQ=255)
                                                  cgtTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCt                                                                                                     <  1:1738057/150‑1 (MQ=255)
                                                   gtTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTg                                                                                                    <  1:3279731/150‑1 (MQ=255)
                                                    ttGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGt                                                                                                   <  1:855757/150‑1 (MQ=255)
                                                     tGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGtt                                                                                                  <  1:789001/150‑1 (MQ=255)
                                                        tCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTtcag                                                                                               >  1:1178751/1‑150 (MQ=255)
                                                        tCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTtcag                                                                                               <  2:3912307/150‑1 (MQ=255)
                                                        tCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTtcag                                                                                               <  1:1887817/150‑1 (MQ=255)
                                                         cccTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTtcagt                                                                                              >  2:2974893/1‑150 (MQ=255)
                                                                aGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAAc                                                                                       >  2:3777206/1‑150 (MQ=255)
                                                                aGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAAc                                                                                       <  2:1773586/150‑1 (MQ=255)
                                                                   aGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAActct                                                                                    <  2:2073544/150‑1 (MQ=255)
                                                                   aGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAActct                                                                                    <  2:172265/150‑1 (MQ=255)
                                                                   aGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAActct                                                                                    <  2:901428/150‑1 (MQ=255)
                                                                   aGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGACCTGTTCAGTCAGGAActct                                                                                    >  1:4264434/1‑150 (MQ=255)
                                                                    gggCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTAAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCtt                                                                                   <  2:4220547/150‑1 (MQ=255)
                                                                      gcggcgCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGc                                                                                 <  2:4232322/150‑1 (MQ=255)
                                                                       cggcgCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCa                                                                                <  1:3445706/150‑1 (MQ=255)
                                                                        ggcgCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAg                                                                               >  2:208657/1‑150 (MQ=255)
                                                                        ggcgCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAg                                                                               <  2:2061734/150‑1 (MQ=255)
                                                                          cgcCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCt                                                                             >  2:1080918/1‑150 (MQ=255)
                                                                           gcCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCtt                                                                            >  1:1399076/1‑150 (MQ=255)
                                                                           gcCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCtt                                                                            <  2:53986/150‑1 (MQ=255)
                                                                              aTTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTtgtg                                                                         <  1:671272/150‑1 (MQ=255)
                                                                              aTTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTtgtg                                                                         <  1:671156/150‑1 (MQ=255)
                                                                              aTTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTAAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTtgtg                                                                         >  2:922228/1‑150 (MQ=255)
                                                                               ttAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTtgtgt                                                                        >  1:731901/1‑150 (MQ=255)
                                                                                    cGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAAt                                                                   >  1:4476530/1‑150 (MQ=255)
                                                                                        gggATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATg                                                               <  2:1429045/150‑1 (MQ=255)
                                                                                         ggATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATgg                                                              <  2:1277765/150‑1 (MQ=255)
                                                                                         ggATGCCGTAATCACTGTTGGTTCGCAGGGAGGTTACAATAAGGTCCATACCGTTACAGTAAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATgg                                                              <  2:305971/150‑1 (MQ=255)
                                                                                          gATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATTCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGTAATCATGGt                                                             <  2:2604684/150‑1 (MQ=255)
                                                                                          gATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGt                                                             >  2:3507135/1‑150 (MQ=255)
                                                                                          gATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGt                                                             >  2:607983/1‑150 (MQ=255)
                                                                                          gATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTAAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGt                                                             >  1:2486793/1‑150 (MQ=255)
                                                                                              ccGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGaa                                                         >  1:3028362/1‑150 (MQ=255)
                                                                                                gTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAACCATGGTTGAAGt                                                       <  1:2974893/150‑1 (MQ=255)
                                                                                                  aaTCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTAAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCg                                                     >  2:3673361/1‑150 (MQ=255)
                                                                                                    tCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTAAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGAc                                                   >  1:215908/1‑150 (MQ=255)
                                                                                                      aCTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTAAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACAt                                                 <  1:2861143/150‑1 (MQ=255)
                                                                                                       cTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTAAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCTACATg                                                >  2:1112882/1‑150 (MQ=255)
                                                                                                           tGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCgg                                            >  1:2032057/1‑150 (MQ=255)
                                                                                                            ggTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGt                                           >  1:1605347/1‑150 (MQ=255)
                                                                                                             gTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTa                                          >  1:4488581/1‑150 (MQ=255)
                                                                                                              ttCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTAc                                         >  1:3502611/1‑150 (MQ=255)
                                                                                                                 gCATGGAGGTTACAATAAGGTCCATACCGTTACAGTAAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGc                                      >  2:2432854/1‑150 (MQ=255)
                                                                                                                   aTGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCAt                                    <  2:4101956/150‑1 (MQ=255)
                                                                                                                   aTGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCAt                                    >  1:1441738/1‑150 (MQ=255)
                                                                                                                    tGGAGGTTACAATAAGGTCCATACCGTTACAGATAAGAGGGATCTCATTCAGACAGAATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATg                                   >  2:2766814/1‑150 (MQ=255)
                                                                                                                       aGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCca                                <  1:4359130/150‑1 (MQ=255)
                                                                                                                         gTTACAATAAGGTCCATACCGTTACAGTAAAGAGGGATCTCTTTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGATCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCcaca                              <  1:2260975/150‑1 (MQ=255)
                                                                                                                            aCAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCCACAAgg                           <  1:88835/150‑1 (MQ=255)
                                                                                                                            aCAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCCACAAgg                           <  1:2075085/150‑1 (MQ=255)
                                                                                                                            aCAATAAGGTCCATACCGTTACAGTAAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCCACCAgg                           >  2:1503297/1‑150 (MQ=255)
                                                                                                                             cAATAAGGCCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCCACAAGGa                          <  1:4078287/150‑1 (MQ=255)
                                                                                                                              aaTAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCCACAAGGAt                         <  2:810414/150‑1 (MQ=255)
                                                                                                                               aTAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCCACAAGGATc                        >  1:2554608/1‑150 (MQ=255)
                                                                                                                                tAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCCACAAGGATCt                       >  2:2630011/1‑150 (MQ=255)
                                                                                                                                tAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCCACAAGGATCt                       <  2:1972909/150‑1 (MQ=255)
                                                                                                                                 aaGGTCCATACCGTTACAGTAAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCCACAAGGATCtt                      >  2:2917796/1‑150 (MQ=255)
                                                                                                                                  aGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCCACAAGGATCttt                     >  2:3411898/1‑150 (MQ=255)
                                                                                                                                       cATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCCACAAGGATCTTTTTCAt                >  2:3458160/1‑150 (MQ=255)
                                                                                                                                       cATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCCACAAGGATCTTTTTCAt                >  1:3050704/1‑150 (MQ=255)
                                                                                                                                        aTACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCCACAAGGATCTTTTTCATa               >  1:3916587/1‑150 (MQ=255)
                                                                                                                                         tACCGTTACAGTAAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCCACAAGGATCTTTTTCATaa              <  2:3172351/150‑1 (MQ=255)
                                                                                                                                          aCCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCCACAAGGATCTTTTTCATAAt             <  2:1903445/150‑1 (MQ=255)
                                                                                                                                          aCCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCCACAAGGATCTTTTTCATAAt             <  1:1028775/150‑1 (MQ=255)
                                                                                                                                           ccGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCCACAAGGATCTTTTTCATAATg            <  1:3701057/150‑1 (MQ=255)
                                                                                                                                            cGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCCACAAGGATCTTTTTCATAATgg           <  2:3812083/150‑1 (MQ=255)
                                                                                                                                                aCAGTTAAGAGGGATCTCATTCAGACAGCAGTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCCACAAGGATCTTTTTCATAATGGAtgt       >  2:1936970/1‑150 (MQ=255)
                                                                                                                                                 cAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCCACAAGGATCTTTTTCATAATGGAtgtg      <  1:805333/150‑1 (MQ=255)
                                                                                                                                                  aGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCCACAAGGATCTTTTTCATAATGGATGTGc     >  2:4409730/1‑150 (MQ=255)
                                                                                                                                                  aGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCCACAAGGATCTTTTTCATAATGGATGTGc     >  1:3766669/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CTCCTTTTATTGAGTTAACAGAGCCTTGATTTGTTGTTTTAACGCGTCGTCGTTGATCCCTGTGAGCAGGGCGGCGCCATTAAGCGTAGGGATGCCGTAATCACTGTTGGTTCGCATGGAGGTTACAATAAGGTCCATACCGTTACAGTTAAGAGGGATCTCATTCAGACAGCATTGGGCGGTAGTGGCTGAAATGCCCTGTTCAGTCAGGAACTCTTGCAGCTTGTGTGCAATCATGGTTGAAGTCGACATGCCGGTACCGCATGCCACAAGGATCTTTTTCATAATGGATGTGCCTC  >  CP009273/4520066‑4520364

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: