Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 853,232 C→A 14.9% intergenic (+221/+20) ybiT → / ← ybiU putative transporter subunit of ABC superfamily: ATP‑binding component/DUF1479 family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP009273853,2320CA14.9% 866.7 / 101.1 279intergenic (+221/+20)ybiT/ybiUputative transporter subunit of ABC superfamily: ATP‑binding component/DUF1479 family protein
Reads supporting (aligned to +/- strand):  ref base C (131/106);  new base A (29/13);  total (160/119)
Fisher's exact test for biased strand distribution p-value = 1.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 284 total aligned reads displayed.

TACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCC  >  CP009273/853083‑853381
                                                                                                                                                     |                                                                                                                                                     
tACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATc                                                                                                                                                       <  2:3984136/150‑1 (MQ=255)
tACGCTTATCAGGCCTACGCAGTACCTGCAAGATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAACGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATc                                                                                                                                                       <  1:1434988/150‑1 (MQ=255)
 aCGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCa                                                                                                                                                      >  2:4402967/1‑150 (MQ=255)
   gCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAAt                                                                                                                                                    <  2:3091090/150‑1 (MQ=255)
    cTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAAtt                                                                                                                                                   <  1:659011/150‑1 (MQ=255)
     ttATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTc                                                                                                                                                  <  2:2611602/150‑1 (MQ=255)
        tCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATAAATTCACt                                                                                                                                               <  2:1905058/150‑1 (MQ=255)
           ggCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTcc                                                                                                                                            <  2:168416/150‑1 (MQ=255)
              cTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCAtt                                                                                                                                         <  1:4020844/150‑1 (MQ=255)
              cTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCAtt                                                                                                                                         >  1:1048703/1‑150 (MQ=255)
                 cGCAGTACCTGCAATCTCTTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATGAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAAt                                                                                                                                      <  1:1893957/150‑1 (MQ=255)
                 cGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAAt                                                                                                                                      <  1:3729992/150‑1 (MQ=255)
                   cAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATAAATTCACTTCCATTAATtc                                                                                                                                    >  2:4357723/1‑150 (MQ=255)
                      tACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCa                                                                                                                                 >  1:3208138/1‑150 (MQ=255)
                         cTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGACTCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATACTCAAAc                                                                                                                              >  1:171465/1‑150 (MQ=255)
                          tGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACa                                                                                                                             >  1:2656688/1‑150 (MQ=255)
                           gCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCCCATCAATTCACTTCCATTAATTCTCAAACAt                                                                                                                            >  1:3789589/1‑150 (MQ=255)
                             aaTATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATcc                                                                                                                          >  1:1005634/1‑150 (MQ=255)
                             aaTATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATcc                                                                                                                          >  1:2180772/1‑150 (MQ=255)
                             aaTATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATcc                                                                                                                          >  1:1573725/1‑150 (MQ=255)
                              atatatTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCa                                                                                                                         >  2:1808759/1‑150 (MQ=255)
                                atatTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATg                                                                                                                       <  2:146989/150‑1 (MQ=255)
                                atatTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATg                                                                                                                       >  1:3818121/1‑150 (MQ=255)
                                  atTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGcc                                                                                                                     >  2:545695/1‑150 (MQ=255)
                                      aaTCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATAAATTCACTTCCATTAATTCTCAAACATCCATGCCCAAc                                                                                                                 >  1:3255486/1‑150 (MQ=255)
                                        tCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAAcgc                                                                                                               >  1:2365743/1‑150 (MQ=255)
                                         cTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACCTCCATGCCCAAcgct                                                                                                              <  2:2144433/150‑1 (MQ=255)
                                         cTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAAcgct                                                                                                              >  2:3883832/1‑150 (MQ=255)
                                         cTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAAcgct                                                                                                              >  1:1298557/1‑150 (MQ=255)
                                            cAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTACCATTAATTCTCAAACAACCAATCCCAAccctcgc                                                                                                           >  1:2444647/1‑150 (MQ=38)
                                                aCTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTAc                                                                                                       >  1:1314632/1‑150 (MQ=255)
                                                 cTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTAcc                                                                                                      >  1:817973/1‑150 (MQ=255)
                                                    tGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATAAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTg                                                                                                   <  1:2190361/150‑1 (MQ=255)
                                                        ggCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATg                                                                                               <  1:4111940/150‑1 (MQ=255)
                                                          ccGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATAAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGtt                                                                                             <  1:4482877/150‑1 (MQ=255)
                                                             gATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGACTCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAGTGTTGAg                                                                                          >  2:2400763/1‑150 (MQ=255)
                                                               tAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGt                                                                                        <  1:1402564/150‑1 (MQ=255)
                                                                aaGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTc                                                                                       <  2:1029664/150‑1 (MQ=255)
                                                                 aGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCg                                                                                      <  1:4244928/150‑1 (MQ=255)
                                                                   gCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGc                                                                                    <  2:1240974/150‑1 (MQ=255)
                                                                     gTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCaa                                                                                  >  2:3306829/1‑150 (MQ=255)
                                                                     gTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCCCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCaa                                                                                  >  2:481435/1‑150 (MQ=255)
                                                                      ttAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAg                                                                                 >  1:1483604/1‑150 (MQ=255)
                                                                      ttAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATAAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAg                                                                                 >  1:434270/1‑150 (MQ=255)
                                                                        aaCGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCg                                                                               >  2:1720816/1‑150 (MQ=255)
                                                                          cgccgcATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATAAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTa                                                                             <  2:1764343/150‑1 (MQ=255)
                                                                             cgcATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATAAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAg                                                                          >  2:1846401/1‑150 (MQ=255)
                                                                                aTCAGTCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGt                                                                       <  2:2308590/150‑1 (MQ=255)
                                                                                   aGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCt                                                                    <  1:2835067/150‑1 (MQ=255)
                                                                                         tAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAg                                                              >  2:1597710/1‑150 (MQ=255)
                                                                                            cAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAgttt                                                           >  1:3617221/1‑150 (MQ=255)
                                                                                             aaaGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATAAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAgtttg                                                          >  2:3123657/1‑150 (MQ=255)
                                                                                             aaaGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATAAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAgtttg                                                          >  1:695692/1‑150 (MQ=255)
                                                                                             aaaGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATAAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAgtttg                                                          >  2:2717630/1‑150 (MQ=255)
                                                                                                 cgcATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTc                                                      >  1:2822963/1‑150 (MQ=255)
                                                                                                 cgcATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCAGGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTc                                                      <  2:3575240/150‑1 (MQ=255)
                                                                                                 cgcATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATAAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTc                                                      >  1:4184296/1‑150 (MQ=255)
                                                                                                  gcATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATAAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCa                                                     >  2:3110954/1‑150 (MQ=255)
                                                                                                   cATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCAt                                                    >  2:1919922/1‑150 (MQ=255)
                                                                                                   cATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATAAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCAt                                                    <  2:634203/150‑1 (MQ=255)
                                                                                                    aTTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATa                                                   <  1:3487275/150‑1 (MQ=255)
                                                                                                    aTTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATa                                                   >  1:507624/1‑150 (MQ=255)
                                                                                                    aTTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATAAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATa                                                   >  1:2328487/1‑150 (MQ=255)
                                                                                                       tGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTc                                                >  2:3762356/1‑150 (MQ=255)
                                                                                                         tCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTccg                                              <  2:2750018/150‑2 (MQ=255)
                                                                                                         tCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCt                                              >  2:823365/1‑150 (MQ=255)
                                                                                                         tCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCt                                              >  2:3262832/1‑150 (MQ=255)
                                                                                                          cAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTc                                             <  1:1959484/150‑1 (MQ=255)
                                                                                                           aGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTcg                                            <  2:2679387/150‑1 (MQ=255)
                                                                                                               aTCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTcgcgcg                                        >  1:130642/1‑150 (MQ=255)
                                                                                                                 cTGAGTAACGTCGAATCGATTTTCTTAATCCCACATAAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAAGGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCAAGTTTGTTTCATAGTCCTCGCGCggg                                      >  2:3967046/1‑150 (MQ=255)
                                                                                                                    aGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTAACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAg                                   >  1:1210135/1‑150 (MQ=255)
                                                                                                                      tAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACCTCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTc                                 <  2:4070424/150‑1 (MQ=255)
                                                                                                                          gTCGAATCGATTTTCTTAATCCCACATAAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGccc                             >  2:320362/1‑150 (MQ=255)
                                                                                                                           tCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCg                            >  2:2788010/1‑150 (MQ=255)
                                                                                                                           tCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCg                            <  1:3289831/150‑1 (MQ=255)
                                                                                                                            cGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCgg                           >  2:1650236/1‑150 (MQ=255)
                                                                                                                             gAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGc                          >  1:3765843/1‑150 (MQ=255)
                                                                                                                                 cGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATg                      >  1:591339/1‑150 (MQ=255)
                                                                                                                                 cGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATg                      >  1:1703084/1‑150 (MQ=255)
                                                                                                                                 cGATTTTCTTAATCCCACATAAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATg                      <  2:3698703/150‑1 (MQ=255)
                                                                                                                                  gATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATgc                     >  1:3632200/1‑150 (MQ=255)
                                                                                                                                  gATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATgc                     <  1:1607165/150‑1 (MQ=255)
                                                                                                                                   aTTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATgcg                    >  1:751309/1‑150 (MQ=255)
                                                                                                                                   aTTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATgcg                    >  2:116977/1‑150 (MQ=255)
                                                                                                                                      ttCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCt                 >  1:1869333/1‑150 (MQ=255)
                                                                                                                                      ttCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCt                 >  2:2652754/1‑150 (MQ=255)
                                                                                                                                         ttAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCtttt              <  1:2469340/150‑1 (MQ=255)
                                                                                                                                          tAATCCCACATAAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCttttt             <  2:434270/150‑1 (MQ=255)
                                                                                                                                           aaTCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTc            <  1:1359842/150‑1 (MQ=255)
                                                                                                                                           aaTCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTc            <  1:1022019/150‑1 (MQ=255)
                                                                                                                                               ccACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGt        >  2:2820524/1‑150 (MQ=255)
                                                                                                                                               ccACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGt        <  2:1300611/150‑1 (MQ=255)
                                                                                                                                               ccACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGt        <  2:2955326/150‑1 (MQ=255)
                                                                                                                                               ccACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGt        <  2:3329806/150‑1 (MQ=255)
                                                                                                                                               ccACATAAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGt        >  2:4001215/1‑150 (MQ=255)
                                                                                                                                                cacaTCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTg       >  1:991984/1‑150 (MQ=255)
                                                                                                                                                 acaTCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGc      >  2:4086311/1‑150 (MQ=255)
                                                                                                                                                 acaTCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGc      <  2:715712/150‑1 (MQ=255)
                                                                                                                                                     cAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGcc  <  2:1987706/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
TACGCTTATCAGGCCTACGCAGTACCTGCAATATATTGAATCTGCAGGACTTTGTAGGCCGGATAAGGCGTTAACGCCGCATCAGGCATTAACAAAGCGCATTTGTCAGCAATCTGAGTAACGTCGAATCGATTTTCTTAATCCCACATCAATTCACTTCCATTAATTCTCAAACATCCATGCCCAACGCTCGCTTACCGTGAATGTTGAGGTCGGCAAGCGTAAAGCGTCCTTCCCAGTTTGTTTCATAGTCCTCGCGCGGGAAGTCGCCCGGCGATGCGCCTTTTTCCAGTGCGGCC  >  CP009273/853083‑853381

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: