Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 2,939,859 Δ1 bp 13.7% coding (679/1098 nt) 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,8590A.13.7% 731.5 / 128.5 249coding (679/1098 nt)mltAmembrane‑bound lytic murein transglycosylase A
Reads supporting (aligned to +/- strand):  ref base A (88/126);  new base . (16/18);  total (104/145)
Fisher's exact test for biased strand distribution p-value = 5.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.80e-01

Only 100 of 252 total aligned reads displayed.

CGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATACCCACTACCCTGCACATCCATAATGA  >  CP009273/2939710‑2940008
                                                                                                                                                     |                                                                                                                                                     
 gcgcgACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGAc                                                                                                                                                      <  2:2857211/150‑1 (MQ=255)
  cgcgACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACa                                                                                                                                                     <  1:2774324/150‑1 (MQ=255)
  cgcgACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACa                                                                                                                                                     <  1:4542682/150‑1 (MQ=255)
  cgcgACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACa                                                                                                                                                     <  2:1299510/150‑1 (MQ=255)
   gcgACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACat                                                                                                                                                    <  1:1641646/150‑1 (MQ=255)
    cgACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACata                                                                                                                                                   <  2:2664608/150‑1 (MQ=255)
    cgACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACata                                                                                                                                                   <  2:3154093/150‑1 (MQ=255)
         aaccaGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATAtcttct                                                                                                                                             <  1:3926098/150‑1 (MQ=255)
              gCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCttttt                                                                                                                                         <  2:2181142/150‑1 (MQ=255)
               cGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTc                                                                                                                                       <  2:91444/150‑1 (MQ=255)
               cGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTc                                                                                                                                       <  2:1444463/150‑1 (MQ=255)
                ggCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTc                                                                                                                                       <  2:3206509/150‑1 (MQ=255)
                ggCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACAGATCTTCTTTTTTc                                                                                                                                       >  2:572656/1‑150 (MQ=255)
                ggCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCa                                                                                                                                      <  2:141066/150‑1 (MQ=255)
                ggCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTc                                                                                                                                       >  2:124021/1‑150 (MQ=255)
                     cGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTc                                                                                                                                  <  2:417315/150‑1 (MQ=255)
                     cGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTc                                                                                                                                  >  2:4714807/1‑150 (MQ=255)
                     cGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACTTCg                                                                                                                                 <  1:4521623/150‑1 (MQ=255)
                       cACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGc                                                                                                                                >  2:4109998/1‑150 (MQ=255)
                        aCTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGcc                                                                                                                               >  1:4679895/1‑150 (MQ=255)
                         cTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCa                                                                                                                              <  2:1144135/150‑1 (MQ=255)
                         cTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCa                                                                                                                              >  1:4621655/1‑150 (MQ=255)
                         cTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCa                                                                                                                              >  1:2688769/1‑150 (MQ=255)
                           tGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACg                                                                                                                            >  1:2298993/1‑150 (MQ=255)
                           tGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTGCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACg                                                                                                                            <  1:2272776/150‑1 (MQ=255)
                             ccccTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGt                                                                                                                          <  2:816784/150‑1 (MQ=255)
                               ccTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCg                                                                                                                        <  2:1019978/150‑1 (MQ=255)
                               ccTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCg                                                                                                                        <  1:355200/150‑1 (MQ=255)
                               ccTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCg                                                                                                                        <  1:4385012/150‑1 (MQ=255)
                                   tCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCa                                                                                                                    <  1:1927286/150‑1 (MQ=255)
                                      ccGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCa                                                                                                                 <  2:3435661/150‑1 (MQ=255)
                                          aGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCtt                                                                                                             >  2:3727644/1‑150 (MQ=255)
                                              aaaGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATTGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCa                                                                                                         >  1:235970/1‑150 (MQ=255)
                                                aGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAAt                                                                                                       <  1:1160507/150‑1 (MQ=255)
                                                    tGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAAtgctg                                                                                                   >  2:4892183/1‑150 (MQ=255)
                                                     gCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAAtgctgc                                                                                                  <  2:2047814/150‑1 (MQ=255)
                                                       ggTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGa                                                                                                >  2:1990426/1‑150 (MQ=255)
                                                       ggTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGa                                                                                                <  2:4701405/150‑1 (MQ=255)
                                                       ggTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGAt                                                                                               <  2:945176/150‑1 (MQ=255)
                                                          ttAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAg                                                                                             <  1:3509782/150‑1 (MQ=255)
                                                            aaaGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGc                                                                                           >  1:1012622/1‑150 (MQ=255)
                                                             aagaagACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCAt                                                                                         <  1:3879763/150‑1 (MQ=255)
                                                               gaagaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATg                                                                                        <  2:2140900/150‑1 (MQ=255)
                                                               gaagaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATg                                                                                        <  2:634941/150‑1 (MQ=255)
                                                                aagaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGa                                                                                       >  1:1558201/1‑150 (MQ=255)
                                                                aagaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTAGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGa                                                                                       <  2:4445529/150‑1 (MQ=255)
                                                                  gaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGAcn                                                                                     <  1:3344839/150‑2 (MQ=255)
                                                                  gaCGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCg                                                                                    <  1:610569/150‑1 (MQ=255)
                                                                    cGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGt                                                                                   >  2:811798/1‑150 (MQ=255)
                                                                     gAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGtt                                                                                  >  1:1533504/1‑150 (MQ=255)
                                                                     gAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGtt                                                                                  >  1:2408176/1‑150 (MQ=255)
                                                                          aCGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTccc                                                                            >  2:852077/1‑150 (MQ=255)
                                                                          aCGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTccc                                                                            >  2:4743325/1‑150 (MQ=255)
                                                                            gggTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCt                                                                           >  2:4002091/1‑150 (MQ=255)
                                                                            gggTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCt                                                                           <  1:2985561/150‑1 (MQ=255)
                                                                             ggTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTg                                                                          >  2:152759/1‑150 (MQ=255)
                                                                             ggTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGc                                                                         >  1:3204791/1‑150 (MQ=255)
                                                                               ttCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCa                                                                        >  1:2806060/1‑150 (MQ=255)
                                                                                  tGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAg                                                                     >  1:209918/1‑150 (MQ=255)
                                                                                     tCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTg                                                                  >  2:782660/1‑150 (MQ=255)
                                                                                      ccAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGa                                                                 <  1:1701232/150‑1 (MQ=255)
                                                                                        agcagcTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGaaa                                                               <  1:3021959/150‑1 (MQ=255)
                                                                                           agcTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAg                                                            >  1:4784902/1‑150 (MQ=255)
                                                                                               cgcgAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTaa                                                        <  1:3604728/150‑1 (MQ=255)
                                                                                                gcgAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTaan                                                       <  1:3389628/150‑2 (MQ=255)
                                                                                                gcgAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAg                                                       <  2:559663/150‑1 (MQ=255)
                                                                                                   aaCCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATTGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCgg                                                    <  2:235970/150‑1 (MQ=255)
                                                                                                   aaCCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCgg                                                    <  2:381053/150‑1 (MQ=255)
                                                                                                     ccTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACt                                                 <  1:67880/150‑1 (MQ=255)
                                                                                                        cGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGc                                               >  2:4760183/1‑150 (MQ=255)
                                                                                                         ggCTTCACGGTGTGTTTCGCCCCAGTGACGAATCGCCGGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGCTCAGCACCGTACCAAGGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACGGcc                                              <  2:3041125/150‑1 (MQ=255)
                                                                                                          gCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCAt                                            >  2:2365595/1‑150 (MQ=255)
                                                                                                            ttCACTGTGTGTTTCGCCCCAGTGACGACTCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCa                                          <  2:2970185/150‑1 (MQ=255)
                                                                                                              cACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACAGATCTTCTTTTTTCAATTCGCCACGGTCGATCAGCAACTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAATCGGACTGCCATCAc                                         >  1:1994306/1‑150 (MQ=255)
                                                                                                                     tgtTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTc                                  <  2:2287876/150‑1 (MQ=255)
                                                                                                                      gtTTAGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCg                                 <  1:124021/150‑1 (MQ=255)
                                                                                                                       tttCGCCCCAGTGACGAATCGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGat                               <  1:1103290/150‑1 (MQ=255)
                                                                                                                        ttCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACGGCCATCACCAAAGTCGat                               <  2:3384056/150‑1 (MQ=255)
                                                                                                                        ttCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACGGCCATCACCAAAGTCGat                               <  1:2886210/150‑1 (MQ=255)
                                                                                                                           gCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACGGCCATCACCAAAGTCGAtata                            >  1:503090/1‑150 (MQ=255)
                                                                                                                            ccacAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATAc                           <  1:1898159/147‑1 (MQ=255)
                                                                                                                             cccAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATAcc                          >  2:4540785/1‑150 (MQ=255)
                                                                                                                               cAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCCTAGCAGAAAAAGTTAAGCGGACGTCCATCACCACAGTCGAGATAcact                        >  2:5073660/1‑147 (MQ=255)
                                                                                                                                  tGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTCAGCGGACGGCCATCACCAAAGTCGATATACCCACTa                     >  2:2869990/1‑150 (MQ=255)
                                                                                                                                      gAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATACCCACTACCCt                 <  1:7624/150‑1 (MQ=255)
                                                                                                                                      gAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACGGCCATCACCAAAGTCGATATACCCACTACCCt                 >  1:4212794/1‑150 (MQ=255)
                                                                                                                                         tCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACGGCCATCACCAAAGTCGATATACCCACTACCCTGca              <  1:1774302/150‑1 (MQ=255)
                                                                                                                                          cGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATACCCACTACCCTGcaca            >  2:1727861/1‑150 (MQ=255)
                                                                                                                                          cGCCTGCATCGCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATACCCACTACCCTGcaca            >  2:1203811/1‑150 (MQ=255)
                                                                                                                                           gCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACGGCCATCACCAAAGTCGATATACCCACTACCCTGcaca            >  2:2287381/1‑150 (MQ=255)
                                                                                                                                           gCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACGGCCATCACCAAAGTCGATATACCCACTACCCTGcaca            <  2:3915834/150‑1 (MQ=255)
                                                                                                                                             cTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATACCCACTACCCTGCACATc          <  1:3695972/150‑1 (MQ=255)
                                                                                                                                             cTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACGGCCATCACCAAAGTCGATATACCCACTACCCTGCACATc          <  1:93490/150‑1 (MQ=255)
                                                                                                                                               gCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACGGCCATCACCAAAGTCGATATACCCACTACCCTGCACATCCa        <  2:1966488/150‑1 (MQ=255)
                                                                                                                                                 aTCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAATCGGACGGCCATCACCAAAGTCGATATACCCACAACCCTGCACATCCATa      >  1:2865780/1‑150 (MQ=255)
                                                                                                                                                 aTCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATACCCACTACCCTGCACATCCATa      >  2:1804693/1‑150 (MQ=255)
                                                                                                                                                 aTCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACGGCCATCACCAAAGTCGATATACCCACTACCCTGCACATCCATa      >  1:2918675/1‑150 (MQ=255)
                                                                                                                                                   cGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATACCCACTACCCTGCACATCCATAAt    <  1:2058160/150‑1 (MQ=255)
                                                                                                                                                   cGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACGGCCATCACCAAAGTCGATATACCCACTACCCTGCACATCCATAAt    >  1:2542957/1‑150 (MQ=255)
                                                                                                                                                     aCATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATACCCACTACCCTGCACATCCATAATGa  <  1:3293382/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCGCGACCAACCAGCGGCACCGCACTTGCCCCTTTCACCGGAGCAAAAGATTGCGGTTTAAAGAAGACGAAAGACGGGTTCTGTTCCAGCAGCTCGCGAACCTCGGCTTCACTGTGTGTTTCGCCCCAGTGACGAATCGCCTGCATCGACATATCTTCTTTTTTCACTTCGCCACGGTCGATCAGCACCTTACCAATGCTGCGATAGGCATGACCGTTTTTCCCTGCATAGCTGAAAAAGTTAAGCGGACTGCCATCACCAAAGTCGATATACCCACTACCCTGCACATCCATAATGA  >  CP009273/2939710‑2940008

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: