Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,828,022 C→T 5.9% intergenic (+168/+31) mtgA → / ← rpoH monofunctional biosynthetic peptidoglycan transglycosylase/RNA polymerase factor sigma‑32

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,828,0220CT 5.9% 469.3 / 12.3 153intergenic (+168/+31)mtgA/rpoHmonofunctional biosynthetic peptidoglycan transglycosylase/RNA polymerase factor sigma‑32
Reads supporting (aligned to +/- strand):  ref base C (74/70);  new base T (3/6);  total (77/76)
Fisher's exact test for biased strand distribution p-value = 3.27e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.61e-01

Only 100 of 155 total aligned reads displayed.

CGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTTCTCTTCCGCCA  >  NC_002947/5827873‑5828171
                                                                                                                                                     |                                                                                                                                                     
   gggTAAACCCGCTCCTACGGGAAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAg                                                                                                                                                    >  1:1493336/1‑150 (MQ=255)
    ggTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGc                                                                                                                                                   >  2:1015968/1‑150 (MQ=255)
    ggTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCACCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGAGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGc                                                                                                                                                   <  2:764866/150‑1 (MQ=255)
     gTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGCAGCGGGCTTACCCGCGAAGAGGCCGGCACAGCCAAACAAAAAAGCCGCTCACCCCCTAGGGCCcgcg                                                                                                                                                  >  1:2188120/1‑150 (MQ=255)
          cccGCTCGTACGGCCAACCGCATACCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGAAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGTCGGCAAAGCCAAACATAAAAGCCGCTGGCCCTATGGGGTGAGCGGCttt                                                                                                                                             <  2:896741/150‑1 (MQ=255)
          cccGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCttt                                                                                                                                             >  1:1439981/1‑150 (MQ=255)
          cccGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCttt                                                                                                                                             >  2:2377660/1‑150 (MQ=255)
           ccGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCtttt                                                                                                                                            >  2:1655911/1‑150 (MQ=255)
             gCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTGGCCCTATGGGGTGAGCGGCtttttt                                                                                                                                          >  2:1668253/1‑148 (MQ=39)
              cTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCtt                                                                                                                                         <  1:1655911/150‑1 (MQ=255)
              cTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCtt                                                                                                                                         <  1:2377660/150‑1 (MQ=255)
               tgcTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCttt                                                                                                                                        <  2:1127452/148‑1 (MQ=255)
                ccTACGGGCAACCGAATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTa                                                                                                                                       <  1:288035/150‑1 (MQ=255)
                 cTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCTGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTGGCCCTATGGGGTGAGCGGCttttttgttt                                                                                                                                      <  1:291552/150‑7 (MQ=37)
                 cTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAGAAAAGCCGCTGGCCCTATGGGGTGAGCGGCttttttgttt                                                                                                                                      <  1:1517162/150‑7 (MQ=37)
                 cTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTGGCCCTATGGGGTGAGCGGCttttttgttt                                                                                                                                      <  1:1668253/150‑7 (MQ=38)
                   aCGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACaa                                                                                                                                    >  2:11468/1‑150 (MQ=255)
                      ggCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTGCAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAAccc                                                                                                                                 <  1:214422/150‑1 (MQ=255)
                       gCAACCGCATCCCCGGCCTCCTCGCGGGCCAAACACCGGCTGGCCATAACCTCAGCCCCTGTAGGAGCGGGAAACCTCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCa                                                                                                                                <  1:1739473/150‑1 (MQ=255)
                               aTCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAAt                                                                                                                        <  2:101295/150‑1 (MQ=255)
                                tCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAAAAGCATCACCTCCACGCGCGCCATCACCAGGGAACCCACCTGCAATGTCTAAGAAAAAAGCCGCTCACCCCAGAGGTCCAGCGGCTTTTCTTTACAACCCAGACTCAATc                                                                                                                       <  2:2079817/150‑1 (MQ=255)
                                       cctcctCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGc                                                                                                                <  1:624525/150‑1 (MQ=255)
                                        ctcctcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTc                                                                                                                                           <  2:945388/122‑1 (MQ=255)
                                        ctcctcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTc                                                                                                                                           >  1:945388/1‑122 (MQ=255)
                                         tcctcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGa                                                                                                              >  1:340725/1‑150 (MQ=255)
                                         tcctcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGa                                                                                                              <  2:1750805/150‑1 (MQ=255)
                                         tcctcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGa                                                                                                              >  1:2058368/1‑150 (MQ=255)
                                            tcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCaa                                                                                                                         >  1:163958/1‑136 (MQ=255)
                                            tcGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCa                                                                                                           <  2:758360/150‑1 (MQ=255)
                                               cgGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTg                                                                                                        <  1:1020773/150‑1 (MQ=255)
                                               cgGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTg                                                                                                        <  1:1015968/150‑1 (MQ=255)
                                                        ccGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTa                                                                                                                                       >  2:1452693/1‑110 (MQ=255)
                                                        ccGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCt                                                                                               >  1:1507857/1‑150 (MQ=255)
                                                           cgcCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGt                                                                                            <  2:1552695/150‑1 (MQ=255)
                                                           cgcCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGt                                                                                            <  1:1001304/150‑1 (MQ=255)
                                                             ccTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTc                                                                                          >  1:1240782/1‑150 (MQ=255)
                                                             ccTACAGGCAACCGCAGCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTc                                                                                          <  2:2341995/150‑1 (MQ=255)
                                                                aCAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATc                                                                                       <  2:981057/150‑1 (MQ=255)
                                                                aCAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATc                                                                                       >  1:369912/1‑150 (MQ=255)
                                                                 cAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGACGAGCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCg                                                                                      >  2:2131660/1‑150 (MQ=255)
                                                                  aGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAACCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCGAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGc                                                                                     >  2:1737653/1‑150 (MQ=255)
                                                                      aaCCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCAttc                                                                                 <  1:95102/150‑1 (MQ=255)
                                                                      aaCCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAAACCAGACTCAATCAGGCAGAGACCAGTGCCTTGACCTTGTTCACCGCCttc                                                                                 >  1:1704380/1‑150 (MQ=255)
                                                                         cgaaTCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAGGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCAttcttc                                                                              >  1:398706/4‑150 (MQ=255)
                                                                          gCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCTGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCt                                                                                               >  2:119048/1‑132 (MQ=255)
                                                                          gCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATtcttct                                                                             <  2:2380858/150‑1 (MQ=255)
                                                                          gCATCCCCTGAAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCACTCAGGCAGCGATCAGTGCCTTGACCTGGTTCATCGCATtcttct                                                                             <  2:216278/150‑1 (MQ=255)
                                                                           cATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCAttcttc                                                                              <  1:2037523/148‑1 (MQ=255)
                                                                               cccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAAc                                                                        <  2:641149/150‑1 (MQ=255)
                                                                               cccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAAc                                                                        <  1:11468/150‑1 (MQ=255)
                                                                                ccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTtctc                                                                            >  1:69047/1‑145 (MQ=255)
                                                                                ccTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACt                                                                       >  1:319262/1‑150 (MQ=255)
                                                                                 cTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTg                                                                      >  2:2276985/1‑150 (MQ=255)
                                                                                 cTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTg                                                                      >  2:107488/1‑150 (MQ=255)
                                                                                   gTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTGGCCCTATGGGGTGAGCGGCtttt                                                                                                                                            >  1:604601/1‑78 (MQ=18)
                                                                                   gTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGc                                                                    <  1:1898190/150‑1 (MQ=255)
                                                                                       gAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAAt                                                                >  1:1067328/1‑150 (MQ=255)
                                                                                            ggCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCt                                                           >  1:1661503/1‑150 (MQ=255)
                                                                                               ttACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCt                                                                                                     <  1:669679/105‑1 (MQ=255)
                                                                                                   ccGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGaa                                                                 >  2:1631884/1‑137 (MQ=255)
                                                                                                   ccGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGaa                                                                 <  1:1631884/137‑1 (MQ=255)
                                                                                                        aaGAGGCCGGCAAAGCCAAAAAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGc                                               <  1:2170254/150‑1 (MQ=255)
                                                                                                          gagGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTg                                             <  1:1176032/150‑1 (MQ=255)
                                                                                                            ggCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTa                                           <  1:1228511/150‑1 (MQ=255)
                                                                                                               cGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCCGGGAACGATTAGTGTTTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACtt                                        >  1:2097442/1‑150 (MQ=255)
                                                                                                               cGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTAc                                          >  2:712170/1‑148 (MQ=255)
                                                                                                               cGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACtt                                        <  1:2276985/150‑1 (MQ=255)
                                                                                                                ggCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTg                                       >  2:1296866/1‑150 (MQ=255)
                                                                                                                 gCAAAGCCAAACAAAAAAGCCTCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGACTACGCCCAGCCGCTACGCTGTACTTGt                                      <  2:1240782/150‑1 (MQ=255)
                                                                                                                  cAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTc                                     <  2:369912/150‑1 (MQ=255)
                                                                                                                   aaaGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCt                                    >  1:729936/1‑150 (MQ=255)
                                                                                                                      gCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATc                                                                                                                       <  1:433262/64‑1 (MQ=255)
                                                                                                                      gCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATc                                                                                                                       >  2:433262/1‑64 (MQ=255)
                                                                                                                      gCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGa                                                                                                  <  2:1617152/85‑1 (MQ=255)
                                                                                                                      gCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGa                                                                                                  >  1:1617152/1‑85 (MQ=255)
                                                                                                                       ccAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGt                                                                                            >  2:2206258/1‑90 (MQ=255)
                                                                                                                       ccAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGt                                                                                            <  1:2206258/90‑1 (MQ=255)
                                                                                                                       ccAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCa                                <  2:2173398/150‑1 (MQ=255)
                                                                                                                        caaacaaaAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACtt                                        >  2:2266073/1‑141 (MQ=255)
                                                                                                                         aaacaaaAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGc                                                                    <  1:574048/112‑1 (MQ=255)
                                                                                                                         aaacaaaAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGc                                                                    >  2:574048/1‑112 (MQ=255)
                                                                                                                              aaaaaGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATg                         <  2:1067328/150‑1 (MQ=255)
                                                                                                                               aaaaGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGc                        >  1:144774/1‑150 (MQ=255)
                                                                                                                                aaaGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTCCCACCCAGACTCAATCAGCCAGCGACAAGTGCCTTAACCTTGTTCATCGCATTCTTCTCCAACTGCCGAATAAGCTCCGCCGCTACGCTCTACTTGTATGCCAGTTCATGcc                       >  2:680052/1‑149 (MQ=255)
                                                                                                                                aaaGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCa                       >  1:1927530/1‑150 (MQ=255)
                                                                                                                                 aaGCCGCTGGCCCTATGGGGTGAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAg                      >  2:739088/1‑150 (MQ=255)
                                                                                                                                   gCCGCTAACCCCATGGGGCCATCGGTTTTGCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCg                    >  2:1956215/1‑150 (MQ=255)
                                                                                                                                    ccGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTAATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGt                   <  2:2058368/150‑1 (MQ=255)
                                                                                                                                      gCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATa                                                  >  2:2235425/1‑117 (MQ=255)
                                                                                                                                          aCCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGt                                                                                                         >  1:2237929/1‑58 (MQ=255)
                                                                                                                                          aCCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGc                                                            >  1:728453/1‑103 (MQ=255)
                                                                                                                                          aCCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGc                                                            <  2:728453/103‑1 (MQ=255)
                                                                                                                                          aCCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCtt             <  2:1555475/150‑1 (MQ=255)
                                                                                                                                          aCCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCtt             <  2:1439981/150‑1 (MQ=255)
                                                                                                                                            cccATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTtct           >  2:836281/1‑150 (MQ=255)
                                                                                                                                             ccATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGGTACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTtc            >  1:622961/1‑148 (MQ=255)
                                                                                                                                             ccATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGGTACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTtc            <  2:622961/148‑1 (MQ=255)
                                                                                                                                               aTAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTTCTCtt        <  2:2127879/150‑1 (MQ=255)
                                                                                                                                                   ggCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGaa                                                                 <  1:987184/89‑1 (MQ=255)
                                                                                                                                                     ccAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTTCTCTTCCgcca  <  2:319262/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCGGGTAAACCCGCTCCTACGGGCAACCGCATCCCCGGCCTCCTCGCGGGCAAACCCGCGCCTACAGGCAACCGCATCCCCTGTAGGAGCGGGCTTACCCGCGAAGAGGCCGGCAAAGCCAAACAAAAAAGCCGCTCACCCCATAGGGCCAGCGGCTTTTCTTTACAACCCAGACTCAATCAGGCAGCGATCAGTGCCTTGACCTTGTTCATCGCATTCTTCTCCAACTGGCGAATACGCTCAGCCGATACGCTGTACTTGTCTGCCAGTTCATGCAGCGTGGCCTTCTCTTCCGCCA  >  NC_002947/5827873‑5828171

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑