Predicted mutation
evidence seq id position mutation freq annotation gene description
RA CP009273 1,238,721 C→A 5.0% A29E (GCG→GAG)  emtA → lytic murein endotransglycosylase E

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*CP0092731,238,7210CA 5.0% 1002.9 / 19.7 258A29E (GCG→GAG) emtAlytic murein endotransglycosylase E
Reads supporting (aligned to +/- strand):  ref base C (128/117);  new base A (8/5);  total (136/122)
Fisher's exact test for biased strand distribution p-value = 5.79e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.00e-01

Only 100 of 267 total aligned reads displayed.

ATCATGCCGCTTAGGTGTGCCGTTGTCACCTCAACGGCGATTCCAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTG  >  CP009273/1238572‑1238870
                                                                                                                                                     |                                                                                                                                                     
atCATGCCGCTTAGGTGTGCCGTTGTCACCTCAACGGCGATTCCAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAAcgc                                                                                                                                                       <  2:2892125/150‑1 (MQ=255)
 tCATGCCGCTTAGGTGTGCCGTTGTCACCTCAACGGCGATTCCAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAAcgcg                                                                                                                                                      >  1:3700079/1‑150 (MQ=255)
 tCATGCCGCTTAGGTGTGCCGTTGTCACCTCAACGGCGATTCCAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAAcgcg                                                                                                                                                      >  1:2962121/1‑150 (MQ=255)
    tGCCGCTTAGGTGTGCCGTTGTCACCTCAACGGCGATTCCAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGaaa                                                                                                                                                   >  2:2983238/1‑150 (MQ=255)
     gCCGCTTAGGTGTGCCGTTGTCACCTCAACGGCGATTCCAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAg                                                                                                                                                  >  1:23374/1‑150 (MQ=255)
      ccGCTTAGGTGTGCCGTTGTCACCTCAACGGCGATTCCAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGAGAAAGt                                                                                                                                                 <  2:4181688/150‑1 (MQ=255)
         cTTAGGTGTGCCGTTGTCACCTCAACGGCGATTCCAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTcc                                                                                                                                              >  2:548237/1‑150 (MQ=255)
          ttAGGTGTGCCGTTGTCACCTCAACGGCGATTCCAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCg                                                                                                                                             >  1:1776004/1‑150 (MQ=255)
           tAGGTGTGCCGTTGTCACCTCAACGGCGATTCCAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCgg                                                                                                                                            >  1:919285/1‑150 (MQ=255)
            aGGTGTGCCGTTGTCACCTCAACGGCGATTCCAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGt                                                                                                                                           >  2:3503556/1‑150 (MQ=255)
              gtgtGCCGTTGTCACCTCACCGGCGATTCCAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGc                                                                                                                                         <  1:1636544/150‑1 (MQ=255)
                 tgCCGTTGTCACCTCAACGGCGATTCCAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAAc                                                                                                                                      <  2:4414982/150‑1 (MQ=255)
                  gCCGTTGTCACCTCAACGGCGATTCCAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACg                                                                                                                                     >  2:3090678/1‑150 (MQ=255)
                       tGTCACCTCAACGGCGATTCCAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGa                                                                                                                                <  1:311173/150‑1 (MQ=255)
                                aaCGGCGATTCCAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGa                                                                                                                       >  2:3917933/1‑150 (MQ=255)
                                aaCGGCGATTCCAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGa                                                                                                                       >  1:1228899/1‑150 (MQ=255)
                                   ggCGATTCCAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGc                                                                                                                    >  1:2862470/1‑150 (MQ=255)
                                   ggCGATTCCAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGc                                                                                                                    <  2:4113274/150‑1 (MQ=255)
                                    gCGATTCCAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGcc                                                                                                                   <  1:2024094/150‑1 (MQ=255)
                                    gCGATTCCAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGcc                                                                                                                   <  2:2165649/150‑1 (MQ=255)
                                      gATTCCAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCaa                                                                                                                 >  2:3174942/1‑150 (MQ=255)
                                      gATTCCAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCaa                                                                                                                 >  1:1417118/1‑150 (MQ=255)
                                          ccAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAg                                                                                                             <  1:4567971/150‑1 (MQ=255)
                                           cAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGc                                                                                                            <  2:2349595/150‑1 (MQ=255)
                                               cTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGa                                                                                                        <  2:1096930/150‑1 (MQ=255)
                                                  taATGATAGAATAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAg                                                                                                     <  1:4086909/150‑1 (MQ=255)
                                                  taAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAg                                                                                                     >  2:3277405/1‑150 (MQ=255)
                                                   aaGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGc                                                                                                    <  2:732829/150‑1 (MQ=255)
                                                    aGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGcc                                                                                                   >  2:141279/1‑150 (MQ=255)
                                                     ggATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCg                                                                                                  >  2:1281025/1‑150 (MQ=255)
                                                       aTAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGt                                                                                                <  2:579826/150‑1 (MQ=255)
                                                          gaagaaGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCa                                                                                             <  1:1056706/150‑1 (MQ=255)
                                                              aagTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCt                                                                                         >  2:1187401/1‑150 (MQ=255)
                                                               agTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTg                                                                                        >  1:849679/1‑150 (MQ=255)
                                                                 tGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTggg                                                                                      <  1:3793981/150‑1 (MQ=255)
                                                                   aaaTTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGc                                                                                    <  1:4387016/150‑1 (MQ=255)
                                                                    aaTTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCg                                                                                   <  1:360784/150‑1 (MQ=255)
                                                                    aaTTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCg                                                                                   >  1:565057/1‑150 (MQ=255)
                                                                     aTTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGt                                                                                  >  1:971608/1‑150 (MQ=255)
                                                                      ttGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTc                                                                                 <  1:3290124/150‑1 (MQ=255)
                                                                      ttGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTc                                                                                 <  2:1844136/150‑1 (MQ=255)
                                                                          gaTGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATc                                                                             <  2:450742/150‑1 (MQ=255)
                                                                           aTGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATcc                                                                            >  2:1496453/1‑150 (MQ=255)
                                                                              gTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCcaca                                                                         <  1:1190245/150‑1 (MQ=255)
                                                                                ttGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAAt                                                                       <  2:4470382/150‑1 (MQ=255)
                                                                                   ccTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGAGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGa                                                                    >  2:2094918/1‑150 (MQ=255)
                                                                                      ttttGATTGTGTTATGAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCa                                                                 <  2:2330224/150‑1 (MQ=255)
                                                                                         tGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACgg                                                              <  1:432309/150‑1 (MQ=255)
                                                                                         tGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACgg                                                              >  1:257391/1‑150 (MQ=255)
                                                                                            ttGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAATCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGa                                                           >  1:1580351/1‑150 (MQ=255)
                                                                                            ttGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGTGTCGATCCACAATTGATCACGGCGa                                                           >  2:3629585/1‑150 (MQ=255)
                                                                                             tgtgTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGAt                                                          <  1:2680630/150‑1 (MQ=255)
                                                                                             tgtgTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGAGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGAt                                                          >  2:2821561/1‑150 (MQ=255)
                                                                                              gtgtTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGAtt                                                         <  2:4345694/150‑1 (MQ=255)
                                                                                               tgtTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTa                                                        <  1:1428936/150‑1 (MQ=255)
                                                                                               tgtTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTa                                                        >  1:3575774/1‑150 (MQ=255)
                                                                                                gtTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTAt                                                       <  1:2094310/150‑1 (MQ=255)
                                                                                                gtTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTAt                                                       <  1:2462479/150‑1 (MQ=255)
                                                                                                   attaGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGc                                                    <  1:499884/150‑1 (MQ=255)
                                                                                                   attaGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGc                                                    >  1:1786965/1‑150 (MQ=255)
                                                                                                   attaGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGc                                                    >  1:4158145/1‑150 (MQ=255)
                                                                                                   attaGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGc                                                    >  2:627344/1‑150 (MQ=255)
                                                                                                   attaGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGAGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGc                                                    >  2:3727433/1‑150 (MQ=255)
                                                                                                    ttaGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCt                                                   >  1:881921/1‑150 (MQ=255)
                                                                                                       gCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATc                                                <  2:942897/150‑1 (MQ=255)
                                                                                                            ttgttCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATc                                           >  2:2403679/1‑150 (MQ=255)
                                                                                                             tgttCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCg                                          >  1:2047303/1‑150 (MQ=255)
                                                                                                                 catcaAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCTGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGggtg                                      >  2:3133251/1‑150 (MQ=255)
                                                                                                                 catcaAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGggtg                                      >  2:3344732/1‑150 (MQ=255)
                                                                                                                  atcaAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGggtgg                                     >  2:1997583/1‑150 (MQ=255)
                                                                                                                    caAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTa                                   >  2:3825129/1‑150 (MQ=255)
                                                                                                                     aaaGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTaa                                  <  2:84959/150‑1 (MQ=255)
                                                                                                                     aaaGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTaa                                  >  1:1763572/1‑150 (MQ=255)
                                                                                                                     aaaGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTaa                                  <  2:1046939/150‑1 (MQ=255)
                                                                                                                      aaGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAAt                                 >  2:4019271/1‑150 (MQ=255)
                                                                                                                      aaGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAAt                                 <  1:187630/150‑1 (MQ=255)
                                                                                                                      aaGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAAt                                 >  1:909471/1‑150 (MQ=255)
                                                                                                                      aaGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAAt                                 >  2:1709340/1‑150 (MQ=255)
                                                                                                                         cATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCg                              <  1:4593486/150‑1 (MQ=255)
                                                                                                                          aTGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGa                             >  2:3372273/1‑150 (MQ=255)
                                                                                                                          aTGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGa                             >  1:1362719/1‑150 (MQ=255)
                                                                                                                           tGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGaa                            <  1:2183940/150‑1 (MQ=255)
                                                                                                                              cTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAAcgc                         <  1:2929238/150‑1 (MQ=255)
                                                                                                                                   cGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGa                    >  1:2072947/1‑150 (MQ=255)
                                                                                                                                    gAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAg                   >  1:3951092/1‑150 (MQ=255)
                                                                                                                                     aaCCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGt                  <  2:1686514/150‑1 (MQ=255)
                                                                                                                                     aaCCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGt                  >  1:2555595/1‑150 (MQ=255)
                                                                                                                                     aaCCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGt                  <  2:2919579/150‑1 (MQ=255)
                                                                                                                                      aCCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTa                 <  1:2283765/150‑1 (MQ=255)
                                                                                                                                        ccgccgTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTaaa               <  1:2604383/150‑1 (MQ=255)
                                                                                                                                        ccgccgTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTaaa               >  2:1367734/1‑150 (MQ=255)
                                                                                                                                         cgccgTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAAt              >  2:673998/1‑150 (MQ=255)
                                                                                                                                          gccgTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCTGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATc             <  1:3148324/150‑1 (MQ=255)
                                                                                                                                          gccgTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATc             <  1:3950539/150‑1 (MQ=255)
                                                                                                                                            cgTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGa           <  2:377193/150‑1 (MQ=255)
                                                                                                                                                 aaCGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGcc      >  1:1567278/1‑150 (MQ=255)
                                                                                                                                                  tcgcgAAAGTTCCTGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCa     <  2:2495935/149‑1 (MQ=255)
                                                                                                                                                  aCGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCa     <  2:4574459/150‑1 (MQ=255)
                                                                                                                                                  aCGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCa     >  1:663251/1‑150 (MQ=255)
                                                                                                                                                     cgAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTg  <  1:2236965/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
ATCATGCCGCTTAGGTGTGCCGTTGTCACCTCAACGGCGATTCCAGGCTATAAGGATAGAAGAAGTGAAATTGAGATGGTTTGCCTTTTTGATTGTGTTATTAGCGGGTTGTTCATCAAAGCATGACTATACGAACCCGCCGTGGAACGCGAAAGTTCCGGTGCAACGTGCGATGCAGTGGATGCCAATAAGCCAGAAAGCCGGTGCAGCCTGGGGCGTCGATCCACAATTGATCACGGCGATTATCGCTATCGAATCGGGTGGTAATCCGAACGCGGTGAGTAAATCGAATGCCATTG  >  CP009273/1238572‑1238870

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: