Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,809,962 A→G 5.6% intergenic (+119/+36) mrcA → / ← maeB penicillin‑insensitive transglycosylase/penicillin‑sensitive transpeptidase/malic enzyme B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,809,9620AG 5.6% 511.0 / 12.8 162intergenic (+119/+36)mrcA/maeBpenicillin‑insensitive transglycosylase/penicillin‑sensitive transpeptidase/malic enzyme B
Reads supporting (aligned to +/- strand):  ref base A (76/75);  new base G (4/5);  total (81/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.22e-01

Only 100 of 163 total aligned reads displayed.

GAAGCGACGCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAACGACGGTGAT  >  NC_002947/5809817‑5810111
                                                                                                                                                 |                                                                                                                                                     
  aGCGACGCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAg                                                                                                                                                 >  1:1772137/1‑150 (MQ=255)
      aCGCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCgg                                                                                                                                             <  2:2980607/150‑1 (MQ=255)
      aCGCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCgg                                                                                                                                             <  1:1267105/150‑1 (MQ=255)
       cGCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCGACATAAAAAAGCCCCGGCTCGCATGAGCCggg                                                                                                                                            <  1:1975136/150‑1 (MQ=39)
       cGCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCGCATGAGCCggg                                                                                                                                            <  2:254641/150‑1 (MQ=255)
       cGCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCggg                                                                                                                                            <  2:1878279/150‑1 (MQ=255)
       cGCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCggg                                                                                                                                            <  1:210212/150‑1 (MQ=255)
       cGCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCggg                                                                                                                                            <  1:1710989/150‑1 (MQ=255)
       cGCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCCGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCggg                                                                                                                                            <  1:2928052/150‑1 (MQ=255)
        gCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCGCATGAGCCgggg                                                                                                                                           <  1:495674/150‑1 (MQ=255)
        gCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCgggg                                                                                                                                           <  1:2916136/150‑1 (MQ=255)
           gATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCtt                                                                                                                                        >  2:3293378/1‑150 (MQ=255)
              ggATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCttttt                                                                                                                                     >  2:1528864/1‑150 (MQ=255)
              ggATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCttttt                                                                                                                                     <  2:155121/150‑1 (MQ=255)
                     ttCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTAtttag                                                                                                                              >  2:2881134/1‑146 (MQ=255)
                      tCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGc                                                                                                                             <  1:3561028/150‑1 (MQ=255)
                      tCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGc                                                                                                                             >  1:795581/1‑150 (MQ=255)
                      tCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTGTTATGTCGc                                                                                                                             >  1:795558/1‑150 (MQ=255)
                       cTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCt                                                                                                                            <  1:899872/150‑1 (MQ=255)
                         aaGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTtg                                                                                                                          >  2:1938928/1‑150 (MQ=255)
                             cAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTtgctgc                                                                                                                      <  2:3576185/150‑1 (MQ=255)
                                gccgcCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCaaa                                                                                                                   >  2:3527287/3‑149 (MQ=255)
                                        aCTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGtt                                                                                                           <  1:571911/150‑1 (MQ=255)
                                          tCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGa                                                                                                         <  2:3755557/150‑1 (MQ=255)
                                          tCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGa                                                                                                         >  2:2539070/1‑150 (MQ=255)
                                          tCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGa                                                                                                         >  2:4436814/1‑150 (MQ=255)
                                          tCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGa                                                                                                         >  1:3285086/1‑150 (MQ=255)
                                           cAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGaa                                                                                                        <  2:1758217/150‑1 (MQ=255)
                                           cAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGaa                                                                                                        >  1:63363/1‑150 (MQ=255)
                                            aGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACACGTGCCCTGACCAACACCTTGCAGCCCTGACACATGCCCACGCCCCATAAAAAAGCCCCGGCTCCTGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAac                                                                                                       <  2:2160555/150‑1 (MQ=255)
                                            aGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAACCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAac                                                                                                       <  1:3097562/150‑1 (MQ=255)
                                               ttttGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAATCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACt                                                                                                    >  1:1857350/1‑150 (MQ=255)
                                                         gcgcTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGc                                                                                          <  2:3000455/150‑1 (MQ=255)
                                                         gcgcTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGc                                                                                          >  2:2627526/1‑150 (MQ=255)
                                                           gctcgctcTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTg                                                                                        <  1:1012530/150‑1 (MQ=255)
                                                           gcTCGATCTGATAGGCGCCACATGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCACGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTg                                                                                        <  1:3527287/150‑1 (MQ=255)
                                                           gcTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTg                                                                                        >  1:2998898/1‑150 (MQ=255)
                                                                 tCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGc                                                                                  <  1:1551122/150‑1 (MQ=255)
                                                                  cTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCg                                                                                 <  2:291105/150‑1 (MQ=255)
                                                                   tGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCgg                                                                                <  2:764347/150‑1 (MQ=255)
                                                                   tGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCgg                                                                                >  1:2279456/1‑150 (MQ=255)
                                                                   tGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCgg                                                                                >  2:2704397/1‑150 (MQ=255)
                                                                    gATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCCGCggg                                                                               >  2:1180896/1‑150 (MQ=255)
                                                                      tAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTa                                                                             <  1:2278238/150‑1 (MQ=255)
                                                                            ccACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCtt                                                                       <  2:3892201/150‑1 (MQ=255)
                                                                             cacaAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAACAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTc                                                                      <  2:3973183/150‑1 (MQ=255)
                                                                             cacaAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTc                                                                      <  1:387485/150‑1 (MQ=255)
                                                                             cacaAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCAGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTc                                                                      <  1:3739463/150‑1 (MQ=255)
                                                                             cacaAGTGCCCTGACCAACACCTTGCAGCCCTGACACAGTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCATTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTc                                                                      <  1:387191/150‑1 (MQ=255)
                                                                              acaAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCTCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCg                                                                     <  1:2026847/150‑1 (MQ=255)
                                                                              acaAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCg                                                                     <  1:4010711/150‑1 (MQ=255)
                                                                              acaAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCg                                                                     <  2:2712124/150‑1 (MQ=255)
                                                                                aaGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGa                                                                   <  1:2866698/150‑1 (MQ=255)
                                                                                aaGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGa                                                                   <  1:1238124/150‑1 (MQ=255)
                                                                                aaGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGa                                                                   <  1:4351207/150‑1 (MQ=255)
                                                                                      ccTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGc                                                             >  2:583425/1‑150 (MQ=255)
                                                                                            cAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCGGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCCGCGGGTAATGCTTCGGATAAGGCAGGGTa                                                       >  2:3057180/1‑150 (MQ=255)
                                                                                            cAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTCATGCTTCGGATAAGGCAGGGTa                                                       >  2:2776128/1‑150 (MQ=255)
                                                                                            cAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACGTCATCCACGCTGGTCAGCGGGTAATGCGTCGGATAAGGCCGGGTa                                                       >  1:1411/1‑150 (MQ=255)
                                                                                            cAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTGTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCGGGTCCGCGGGTCATGCTTCGGATAAGGCAGGGtc                                                       >  2:2072771/1‑149 (MQ=255)
                                                                                            cAACACCTTGCAGACCTGACACATTCCCAATCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGTGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTa                                                       <  2:2740134/150‑1 (MQ=255)
                                                                                                   ttGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGc                                                <  1:4258606/150‑1 (MQ=255)
                                                                                                     gCAGCCCTGACACATTCCCAAGCCCCATAACAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCg                                              <  1:210780/150‑1 (MQ=255)
                                                                                                       aGCCCTGACACATTCCCAAGCCCCATAAAAATGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGa                                            >  1:3467697/1‑150 (MQ=255)
                                                                                                       aGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGa                                            >  1:1027039/1‑150 (MQ=255)
                                                                                                       aGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCATTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGa                                            <  2:2559787/150‑1 (MQ=255)
                                                                                                          ccTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCCTCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTc                                         <  2:4230590/150‑1 (MQ=255)
                                                                                                            tGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGa                                       >  1:2980634/1‑150 (MQ=255)
                                                                                                            tGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGa                                       >  1:3327402/1‑150 (MQ=255)
                                                                                                              acacATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGACCACTTCATCCACGCTGGTCAGCGGGTCATGCTTCGGATCAGGCCGGGTCGCCACGCCGGATTCGagg                                     >  1:3332322/1‑148 (MQ=255)
                                                                                                               cacaTTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATgg                                    <  1:575499/150‑1 (MQ=255)
                                                                                                                 caTTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCa                                  <  2:1914894/150‑1 (MQ=255)
                                                                                                                 caTTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCa                                  <  1:1911616/150‑1 (MQ=255)
                                                                                                                  aTTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAg                                 >  2:4489972/1‑150 (MQ=255)
                                                                                                                  aTTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAg                                 >  1:4176592/1‑150 (MQ=255)
                                                                                                                    tCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGcc                               >  1:909764/1‑150 (MQ=255)
                                                                                                                           ccccATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCca                        >  1:1453463/1‑150 (MQ=255)
                                                                                                                           ccccATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCca                        >  1:3729481/1‑150 (MQ=255)
                                                                                                                            cccATCACACAGCCCCGGCTCATGCGAGCCGGGGCTTATTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCcac                       <  1:4402617/150‑1 (MQ=255)
                                                                                                                             ccATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCcaca                      <  2:2682846/150‑1 (MQ=255)
                                                                                                                               aTAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGc                    >  1:2402072/1‑150 (MQ=255)
                                                                                                                                tAAAAAAGCCCCGGCTCGCATGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGGAGCCACGCCGGATTCGATGGCAGCCTTGGCCACGGCa                   >  2:16859/1‑150 (MQ=38)
                                                                                                                                tAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCa                   >  1:3059507/1‑150 (MQ=255)
                                                                                                                                tAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCa                   >  1:789903/1‑150 (MQ=255)
                                                                                                                                  aaaaaGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATc                 >  2:2627592/1‑150 (MQ=255)
                                                                                                                                  aaaaaGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATc                 >  1:2889337/1‑150 (MQ=255)
                                                                                                                                   aaaaGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCg                <  2:4430237/150‑1 (MQ=255)
                                                                                                                                    aaaGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCgg               >  1:2218499/1‑150 (MQ=255)
                                                                                                                                        ccccGGCTCGCATGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAac           >  2:25881/1‑150 (MQ=255)
                                                                                                                                        ccccGGCTCGCATGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAac           >  2:2566947/1‑150 (MQ=255)
                                                                                                                                        ccccGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAac           >  1:1060082/1‑150 (MQ=255)
                                                                                                                                        ccccGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAac           >  1:2718583/1‑150 (MQ=255)
                                                                                                                                        ccccGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAac           >  2:3172106/1‑150 (MQ=255)
                                                                                                                                         cccGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAacg          >  2:3623687/1‑150 (MQ=255)
                                                                                                                                          ccGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAacga         >  1:1925/1‑150 (MQ=255)
                                                                                                                                          ccGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAacga         >  1:3790716/1‑150 (MQ=255)
                                                                                                                                          ccGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAacga         >  1:711893/1‑150 (MQ=255)
                                                                                                                                             gCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCAGCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAACGACgg      <  2:1772137/150‑1 (MQ=255)
                                                                                                                                                cATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAACGACGGTGa   <  1:2826134/150‑1 (MQ=255)
                                                                                                                                                cATGCGAGCCGGGGATTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGTATAAGGCTGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAACGACGGTGa   <  1:602555/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GAAGCGACGCCGATGGATCTGTTCTAAGGCAGGGCGCCTGACTCAGGTTTTGCGTGGGCGCTCGATCTGATAGGCGCCACAAGTGCCCTGACCAACACCTTGCAGCCCTGACACATTCCCAAGCCCCATAAAAAAGCCCCGGCTCATGCGAGCCGGGGCTTTTTTATGTCGCTTGCTGCAATCAGCCGTTGAACACTTCATCCACGCTGGTCAGCGGGTAATGCTTCGGATAAGGCAGGGTAGCCACGCCGGATTCGATGGCAGCCTTGGCCACAGCATCGGAAACGACGGTGAT  >  NC_002947/5809817‑5810111

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 25 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: