Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,939,832 C→A 14.1% E236* (GAA→TAA)  mltA ← membrane‑bound lytic murein transglycosylase A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092732,939,8320CA14.1% 714.6 / 63.5 227E236* (GAA→TAA) mltAmembrane‑bound lytic murein transglycosylase A
Reads supporting (aligned to +/- strand):  ref base C (92/103);  new base A (13/19);  total (105/122)
Fisher's exact test for biased strand distribution p-value = 5.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 231 total aligned reads displayed.

ATAATGGAACGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATA  >  CP009273/2939684‑2939980
                                                                                                                                                    |                                                                                                                                                    
ataatGGAACGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTc                                                                                                                                                      <  1:2927452/149‑1 (MQ=255)
      gAACGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCa                                                                                                                                                >  1:3447948/1‑149 (MQ=255)
      gAACGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCa                                                                                                                                                >  1:1854838/1‑149 (MQ=255)
           aTCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGAc                                                                                                                                           >  2:3532043/1‑149 (MQ=255)
             cAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGa                                                                                                                                         <  1:3680900/149‑1 (MQ=255)
             cAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGa                                                                                                                                         >  2:1575940/1‑149 (MQ=255)
              agagGCATCTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGaa                                                                                                                                        <  2:726207/149‑1 (MQ=255)
              agagGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGaa                                                                                                                                        <  2:51516/149‑1 (MQ=255)
               gagGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAAt                                                                                                                                       >  1:764490/1‑149 (MQ=255)
               gagGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGAAGAAAGACGGGTTCTGTTCCAGCAGCTCGCGGACCGCGGCTTCACTGTGGGTTTAGCCCCAGTGACGAAt                                                                                                                                       >  2:3854947/1‑149 (MQ=255)
                agGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATc                                                                                                                                      <  2:2969671/149‑1 (MQ=255)
                   cAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGcc                                                                                                                                   <  2:2186826/149‑1 (MQ=255)
                   cAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGcc                                                                                                                                   <  2:821450/149‑1 (MQ=255)
                   cAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGcc                                                                                                                                   <  1:3476212/149‑1 (MQ=255)
                   cAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGcc                                                                                                                                   <  1:2455297/149‑1 (MQ=255)
                   cAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGcc                                                                                                                                   <  1:239614/149‑1 (MQ=255)
                    aaCTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCt                                                                                                                                  >  2:1726435/1‑149 (MQ=255)
                      cTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGc                                                                                                                                >  2:2809528/1‑149 (MQ=255)
                          cgcgcgACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGATCAAAAGATTGCGGTTTTAAGAAGACGATAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCTTCg                                                                                                                            <  2:2402934/149‑1 (MQ=255)
                               gACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACata                                                                                                                       <  2:2628402/149‑1 (MQ=255)
                               gACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATCGACata                                                                                                                       <  1:2516058/149‑1 (MQ=255)
                                accaaccaGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACatat                                                                                                                      >  1:1878326/1‑149 (MQ=255)
                                 ccaaccaGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATAtc                                                                                                                     >  1:3452157/1‑149 (MQ=255)
                                  caaccaGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGATGCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGACTCGCCTGCATCGACATAtct                                                                                                                    <  2:1836321/149‑1 (MQ=255)
                                      caGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATcttctt                                                                                                                <  1:2415903/149‑1 (MQ=255)
                                      caGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCttt                                                                                                               <  2:1515938/149‑1 (MQ=255)
                                      caGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCttt                                                                                                               <  2:3564017/149‑1 (MQ=255)
                                        gCGGCACCGCAATTGCCCCTTTCACCGGAGCAAAAGATTGCTGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCtttt                                                                                                              <  1:1428707/149‑1 (MQ=255)
                                           gCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTc                                                                                                           >  2:1547002/1‑149 (MQ=255)
                                           gCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTc                                                                                                           <  1:1075684/149‑1 (MQ=255)
                                           gCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTc                                                                                                           >  2:122585/1‑149 (MQ=255)
                                           gCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTc                                                                                                           >  1:1906884/1‑149 (MQ=255)
                                            cACCGCACGTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCa                                                                                                          <  2:3445096/149‑1 (MQ=255)
                                             aCCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCAc                                                                                                         >  1:1726352/1‑149 (MQ=255)
                                             aCCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCAc                                                                                                         <  2:2701235/149‑1 (MQ=255)
                                               cGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACtt                                                                                                       <  2:3685779/149‑1 (MQ=255)
                                               cGCACTTGCCCCTTTCACCGGAGCAAAAGAGTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACtt                                                                                                       <  2:1006754/149‑1 (MQ=255)
                                                  aCTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAATGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGc                                                                                                    >  2:887251/1‑149 (MQ=255)
                                                      gCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACg                                                                                                >  2:1091342/1‑149 (MQ=255)
                                                         ccTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTc                                                                                             <  1:1547002/149‑1 (MQ=255)
                                                         ccTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTc                                                                                             <  1:3110592/149‑1 (MQ=255)
                                                         ccTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTc                                                                                             <  1:3110939/149‑1 (MQ=255)
                                                         ccTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTc                                                                                             >  1:2910007/1‑149 (MQ=255)
                                                           tttCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGa                                                                                           >  1:2186767/1‑149 (MQ=255)
                                                           tttCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGa                                                                                           <  1:2819247/149‑1 (MQ=255)
                                                             tCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATc                                                                                         <  2:2203252/149‑1 (MQ=255)
                                                                 cGGAGCAAAAGCTTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCa                                                                                     <  2:2715215/149‑1 (MQ=255)
                                                                    aGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCt                                                                                  <  2:3650206/149‑1 (MQ=255)
                                                                     gCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCtt                                                                                 >  2:1657058/1‑149 (MQ=255)
                                                                      cAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTa                                                                                <  2:928124/149‑1 (MQ=255)
                                                                      cAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTa                                                                                <  2:2986943/149‑1 (MQ=255)
                                                                        aaaGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTAcc                                                                              >  2:3361211/1‑149 (MQ=255)
                                                                        aaaGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACAGATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTAcc                                                                              >  2:1844728/1‑149 (MQ=255)
                                                                          aGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCaa                                                                            <  1:306564/149‑1 (MQ=255)
                                                                           gATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAAt                                                                           <  1:406343/149‑1 (MQ=255)
                                                                           gATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAAt                                                                           >  2:2577072/1‑149 (MQ=255)
                                                                           gATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAAt                                                                           <  1:1303866/149‑1 (MQ=255)
                                                                                 ggTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCg                                                                     >  1:2541786/1‑149 (MQ=255)
                                                                                      aaaGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAgg                                                                >  2:430643/1‑149 (MQ=255)
                                                                                      aaaGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAgg                                                                >  1:131471/1‑149 (MQ=255)
                                                                                       aagaagACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGc                                                               <  1:3583405/149‑1 (MQ=255)
                                                                                       aagaagACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGc                                                               <  2:2050863/149‑1 (MQ=255)
                                                                                        agaagaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCa                                                              >  2:3660466/1‑149 (MQ=255)
                                                                                        agaagaCGAAAGACGGGTGCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCa                                                              >  2:2105047/1‑149 (MQ=255)
                                                                                          aagaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATg                                                            >  1:1573197/1‑149 (MQ=255)
                                                                                          aagaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATg                                                            >  1:2962331/1‑149 (MQ=255)
                                                                                           agaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATTCTGCGATAGGCATGa                                                           >  2:3666874/1‑149 (MQ=255)
                                                                                           agaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGa                                                           <  1:544707/149‑1 (MQ=255)
                                                                                            gaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGAc                                                          >  2:3539437/1‑149 (MQ=255)
                                                                                            gaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGAc                                                          >  2:2073666/1‑149 (MQ=255)
                                                                                               gAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGt                                                       <  1:3539437/149‑1 (MQ=255)
                                                                                                 aaGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGttt                                                     >  2:654472/1‑149 (MQ=255)
                                                                                                  aGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGtttt                                                    >  1:3938031/1‑149 (MQ=255)
                                                                                                     cGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTcc                                                 >  2:3509960/1‑149 (MQ=255)
                                                                                                      gggTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTccc                                                <  2:2349019/149‑1 (MQ=255)
                                                                                                       ggTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCt                                               >  1:3701845/1‑149 (MQ=255)
                                                                                                       ggTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGACTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCt                                               <  2:2910007/149‑1 (MQ=255)
                                                                                                         ttCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGc                                             <  2:74264/149‑1 (MQ=255)
                                                                                                          tCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCAt                                           <  1:3340218/149‑1 (MQ=255)
                                                                                                           cTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCAt                                           >  1:2586424/1‑149 (MQ=255)
                                                                                                           cTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCAt                                           >  2:2457946/1‑149 (MQ=255)
                                                                                                           cTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATa                                          >  2:1941643/1‑149 (MQ=255)
                                                                                                            tGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATa                                          <  1:3020856/149‑1 (MQ=255)
                                                                                                              ttCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGAATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGc                                        >  2:2149078/1‑149 (MQ=255)
                                                                                                                 cagcagCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGa                                     >  1:458261/1‑149 (MQ=255)
                                                                                                                    cagcTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGaaaa                                  <  1:982887/149‑1 (MQ=255)
                                                                                                                    cagcTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCGTCTTTTTTCACTTCGCCACGCACGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGaaaa                                  <  1:4035961/149‑1 (MQ=255)
                                                                                                                       cTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGt                               <  1:2544330/149‑1 (MQ=255)
                                                                                                                          gcgAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAATAGTTaa                            >  2:3808725/1‑149 (MQ=255)
                                                                                                                            gAACCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGc                          >  1:3198851/1‑149 (MQ=255)
                                                                                                                             aaCCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCg                         <  1:1698177/149‑1 (MQ=255)
                                                                                                                             aaCCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCg                         >  2:3903293/1‑149 (MQ=255)
                                                                                                                               ccTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGa                       <  1:627793/149‑1 (MQ=255)
                                                                                                                                 tCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACt                     >  2:1328662/1‑149 (MQ=255)
                                                                                                                                      ttCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCAt                >  1:577502/1‑149 (MQ=255)
                                                                                                                                         aCTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCAcc            >  2:1401506/1‑149 (MQ=255)
                                                                                                                                              gtgtTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAg        <  1:1726435/149‑1 (MQ=255)
                                                                                                                                                   tCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGat   >  2:418064/1‑149 (MQ=255)
                                                                                                                                                   tCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGat   >  1:1201516/1‑149 (MQ=255)
                                                                                                                                                    cGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGata  <  2:2962331/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
ATAATGGAACGATCAGAGGCAACTGACGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATA  >  CP009273/2939684‑2939980

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: