Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,863,154 A→C 11.2% intergenic (‑138/+99) ligA ← / ← PP_4275 DNA ligase/cell division protein ZipA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,863,1540AC11.2% 581.0 / 48.8 221intergenic (‑138/+99)ligA/PP_4275DNA ligase/cell division protein ZipA
Reads supporting (aligned to +/- strand):  ref base A (130/59);  new base C (5/21);  total (141/80)
Fisher's exact test for biased strand distribution p-value = 1.72e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.90e-01

Only 100 of 222 total aligned reads displayed.

TCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATACGCTGGCG  >  NC_002947/4863006‑4863303
                                                                                                                                                      |                                                                                                                                                     
tCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTAc                                                                                                                                                      <  2:8434926/150‑1 (MQ=255)
tCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTAc                                                                                                                                                      <  1:3798055/150‑1 (MQ=255)
 cggcggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCACAAGGGCGGGCCATTCCCTAGGGCGGTGGTGTAcc                                                                                                                                                     >  1:526622/1‑150 (MQ=255)
  ggcggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCt                                                                                                                                                    <  1:4189615/150‑1 (MQ=255)
          tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCACCCACAAGGGCGTCACATGCCCAAGGGCGGTGATGTACCAGTAGGAGc                                                                                                                                            >  1:7422507/1‑150 (MQ=255)
          tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAAACCCACCCACAAAGGCGTCACAATCCCAAGGGCGGGGATGTAACAGTATGAGc                                                                                                                                            >  1:4649316/1‑150 (MQ=255)
            ggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTTAACCCCCGCCCACAAAGGCGGCACAAGCCCAAGGGCGGGGTTGTAACAGTAGGAGCgg                                                                                                                                          >  1:8420756/1‑150 (MQ=255)
            ggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCgg                                                                                                                                          >  2:43006/1‑150 (MQ=255)
            ggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCCCGCCCACCACGGCGTGACAAGCCCTACGGCGGGGATGTAACAGTATGAGCgg                                                                                                                                          >  1:8421440/1‑150 (MQ=255)
             gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCggg                                                                                                                                         <  2:4539692/150‑1 (MQ=255)
             gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCggg                                                                                                                                         <  2:7165788/150‑1 (MQ=255)
             gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCCTGAGGGCGTGCTGGCCGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCgg                                                                                                                                          <  1:2046957/150‑1 (MQ=37)
              aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGt                                                                                                                                        <  1:6096719/150‑1 (MQ=14)
              aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGt                                                                                                                                        <  1:7090696/150‑1 (MQ=255)
                ttGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGttt                                                                                                                                      <  2:312793/150‑1 (MQ=14)
                  gTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTAc                                                                                                                                    <  2:1049204/150‑1 (MQ=14)
                   tcttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTAcc                                                                                                                                   >  2:7743231/1‑150 (MQ=255)
                   tcttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTAcc                                                                                                                                   >  2:5039004/1‑150 (MQ=255)
                   tcttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGCGCGGGTTTAcc                                                                                                                                   >  2:7085723/1‑150 (MQ=255)
                    cttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTATCTGTAGGAGCGGGTTTAccc                                                                                                                                  >  2:2728755/1‑150 (MQ=255)
                    cttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTAccc                                                                                                                                  >  1:3526456/1‑150 (MQ=255)
                     ttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCcg                                                                                                                                 <  1:9027265/150‑1 (MQ=14)
                     ttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCcg                                                                                                                                 >  2:6822529/1‑150 (MQ=255)
                     ttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCcg                                                                                                                                 >  2:1882831/1‑150 (MQ=255)
                     ttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCcg                                                                                                                                 >  1:3744103/1‑150 (MQ=255)
                      tctcTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATTCCCTAGGGCGGGGGTGTACCGGTAGGTGCGGGGTTTCCcgc                                                                                                                                >  1:6572425/1‑150 (MQ=38)
                      tctcTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCcgc                                                                                                                                <  2:3378502/150‑1 (MQ=255)
                        tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGa                                                                                                                              >  2:6566248/1‑150 (MQ=255)
                        tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCAAATGCCCTAGGGCGGTGATTTACCTGTAGGAGCGGGGTTAACCGCGa                                                                                                                              >  1:2075651/1‑150 (MQ=255)
                         ctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGaa                                                                                                                             <  2:2758813/150‑1 (MQ=11)
                         ctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                             >  2:2828032/1‑150 (MQ=255)
                          ttGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGAag                                                                                                                            >  2:7644688/1‑150 (MQ=255)
                          ttGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCAAAATGGCGTGACATGCCCTAGGGCGGTGTTGTACCTGTAGGAGCGGGTTTTACCGCGaat                                                                                                                            >  1:6837574/1‑149 (MQ=255)
                            gcaaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGGGGTGATGTACCTGTAGGAGCGGTTTTACCCGCGAagag                                                                                                                          <  2:1102258/150‑1 (MQ=255)
                              aaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCACAAGGGCGTCACATTCCCTAGGGCGGTGTTGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                                                        >  1:5945873/1‑150 (MQ=255)
                                  caaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGTAAACCGCACCCACAAGGGCGGCACAAGCCCAAGGGCGGGGGTGGACCAGGATGGGCGGGGTTTCCCCCGAAGAGGCCGGc                                                                                                                    >  1:6215678/1‑150 (MQ=255)
                                  caaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGAAAACCCCTCCCACAAGGGCGGCCCAATCCCGAGGGCGGGGATGTACCGGTGGGAGCGGGTTTAACCGCGAAGAGGCCGGc                                                                                                                    >  1:3575423/1‑150 (MQ=255)
                                   aaaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGGTTACCCGCGAAGAGGCCGGca                                                                                                                   >  2:4956688/1‑150 (MQ=255)
                                    aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTTAACCCGCTCCCACAAGGGCGGCACATGCCCTAGGGCGGTGATGTACCGGTTGGAGCGGGTTTTACCCCGCAGAGGCCGGcac                                                                                                                  >  1:5157332/1‑150 (MQ=255)
                                    aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGGAAACCCCGCCCCCAAGGGCGGCCCCGGCCCCAGGGCGGGGGTGTTCCGGGAGGGGCGGGGTTACACCCGAAGAGggcggccc                                                                                                                  >  1:121939/1‑148 (MQ=255)
                                     aaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTACCAAAAAGGCGGCACAAGCCCAAGGGCGGGGATGTACCAGGAGGAGCGGGTTTAACCGCGAAGAGGCCGGcaca                                                                                                                 >  1:8835952/1‑150 (MQ=255)
                                         gctgcTTGGGGTGGTGCGTGGGGTTCCCGGTGTCTGTTGCCTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACACGGGCGTCACTGGCCCTTGGGGGGTGATGTACCTGTCGGTGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCa                                                                                                             <  2:3575423/147‑1 (MQ=37)
                                         gcagcTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAAGCCCGAGGGCGGTGATGTACCAGTAGGAGCGGGTTTACCCGCGAAGTGGCCGGCACAGGCa                                                                                                             >  2:822957/1‑150 (MQ=255)
                                           agcTTGGGGTGGTGCGTGGGGTTCCTGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAac                                                                                                           <  2:6694414/150‑1 (MQ=255)
                                           agcTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCGAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAac                                                                                                           >  2:4740935/1‑150 (MQ=255)
                                                  ggtggtGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTGAACCCCCGCCCACAAGGGCGGCACATGCCCAAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCCCAGGCAACCACGATa                                                                                                    >  1:2798588/1‑150 (MQ=38)
                                                            tGGGGTTCCTGGTGTCTGTTGACTGTGCCAGCCTCTGCGCGGGAAACCCCGCGCCCACACGGGCGTCACTGGCCCTTGGGCGGTTTTGTACCTGTCGGAGCGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCg                                                                                           <  2:2143447/150‑1 (MQ=37)
                                                                gTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTATACCCGGTGCCACAAGGGCGTCACATGGCCTTGGGGGGTGTTGTTGGTGTTGGTGCGGGTTTTCCCGCGAAGTGGGCGGGACAGGGAACGTCGTTAGGCGAAGCgggct                                                                                       >  2:2941831/1‑147 (MQ=255)
                                                                                   ggggCCTGCGTCTTCGGGGGTAACCCGGCTCCCACCGGGGGGTCACTGGCCCTTGGGGGGTTCTTTACGTTTTGGAGGGGGTTCACCCGCTAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGctc                                                                    <  2:4009050/147‑1 (MQ=255)
                                                                                          gCCTCTTCGCGGGTAAACCCCCTCCCACAAGGGCGTCACATTCCCAAGGGCGGTGATGTACCGGTAGGAGCGGGTTTAACCCCGAAGAGGCCGGCACAGGCAAAGACCATAACCCAACCGGCCAGCACCCCCACACGCAGCTCTTTCGta                                                             >  1:7551236/1‑149 (MQ=25)
                                                                                           ccTCTCCGGGGTTAAACCCGCCCCCACACGGGCGACACTTGCCCTCGGGGGGTGATGTCCCTTTAGGAGGGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCt                                                            <  2:2776520/150‑1 (MQ=255)
                                                                                              ctTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCggg                                                                                                                                         >  2:5663207/1‑70 (MQ=255)
                                                                                               ttCGGGGGCAAACCCTCTCCCACAAGGGCGTCACATGCCCTAGGGCGTTTATTTACCTGTAGTAGCGGGTTTACCCGCGAAGTGGCCGGCACATGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGAt                                                        <  2:9135782/150‑1 (MQ=255)
                                                                                               ttCGCGGGTAAACCCGCGCCCCCCAGGGCGTCACATGCCCTAGGGCGGTGAGGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGAt                                                        <  2:5156529/150‑1 (MQ=255)
                                                                                                 ggcgGTTAACCCCGCCCCCACTGGGGCTCCACAGGCCCTTGGGGGGTTTTTTCCCTGTTGGTGGGGTTTAACCGGCTAAGCGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCt                                                      <  1:2937173/149‑1 (MQ=255)
                                                                                                 cgggggAAACCCCGCCCCCACACGGGCGACACAGGCCCTTGGGCGGTTATGTACCTGTCGGAGGGGGTTAACCCGCGAAGTGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCt                                                      <  2:1585001/147‑1 (MQ=255)
                                                                                                  gcgGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTg                                                     >  1:7222607/1‑150 (MQ=255)
                                                                                                  gcgGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTg                                                     >  1:3135993/1‑150 (MQ=255)
                                                                                                   gggtaaaaCCCGCCCCCCCACGGGCGCCACAGGCCCTGGGGGGGTGCTGGACCTGTCGGTGGGTTTTCCCCCGCGAAGTGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGt                                                    <  2:1720215/145‑1 (MQ=255)
                                                                                                     ggaaaccccGCTCCCACACGGGCGTCACTGGCCCTCGGGCGGTGATGTCCCTGTTGGCGCGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCa                                                  <  2:6176391/144‑1 (MQ=38)
                                                                                                     ggTAACCCCGCCCCCACCCGGGGGTCACTTGCCCTTGGGGGGTTTTGTACGTGTTGGTGGGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCa                                                  <  2:6193788/150‑1 (MQ=255)
                                                                                                       tAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgc                                                >  1:6460338/1‑150 (MQ=255)
                                                                                                        aaaCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgcg                                               >  2:8985171/1‑150 (MQ=255)
                                                                                                        aaaCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgcg                                               >  2:4780932/1‑150 (MQ=255)
                                                                                                         aaCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgcgc                                              >  1:5267095/1‑150 (MQ=255)
                                                                                                          aCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt                                             >  2:8714532/1‑150 (MQ=255)
                                                                                                          aCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt                                             >  2:23550/1‑150 (MQ=255)
                                                                                                          aCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt                                             >  2:8271483/1‑150 (MQ=255)
                                                                                                          aCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt                                             >  2:578276/1‑150 (MQ=255)
                                                                                                           cccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCtt                                            >  1:6741777/1‑150 (MQ=255)
                                                                                                           cccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCtt                                            >  2:1154496/1‑150 (MQ=255)
                                                                                                            ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCTCAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                                           >  2:3263958/1‑150 (MQ=255)
                                                                                                            ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                                           >  2:8832685/1‑150 (MQ=255)
                                                                                                            ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                                           >  1:7803188/1‑150 (MQ=255)
                                                                                                              gCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTg                                         >  1:935222/1‑150 (MQ=255)
                                                                                                                tCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTggg                                       >  1:6218933/1‑150 (MQ=255)
                                                                                                                  ccACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGtt                                     >  1:2167764/1‑150 (MQ=255)
                                                                                                                   cacaCGGGCGTCACATGCCCTTGGGGGGTTTTGTACCTTTTGGAGGGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTa                                    <  2:2719906/150‑1 (MQ=38)
                                                                                                                    acaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAg                                   >  1:3381607/1‑150 (MQ=255)
                                                                                                                    acaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAg                                   >  1:7671517/1‑150 (MQ=255)
                                                                                                                    acaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAg                                   >  1:8291531/1‑150 (MQ=255)
                                                                                                                       aGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAgcgc                                >  2:885140/1‑150 (MQ=255)
                                                                                                                        gggCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCa                               >  1:8194928/1‑150 (MQ=255)
                                                                                                                            gTCACATGCCCTAGGGCGGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACACGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGgc                           >  1:1878261/1‑150 (MQ=38)
                                                                                                                               caccgCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGgcgc                         >  1:6346179/5‑150 (MQ=25)
                                                                                                                                caTGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTc                       >  1:7196183/1‑150 (MQ=255)
                                                                                                                                 aTGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCg                      >  2:6096838/1‑150 (MQ=255)
                                                                                                                                 aTGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCg                      >  2:1509588/1‑150 (MQ=255)
                                                                                                                                    cccTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAAc                   >  2:337602/1‑150 (MQ=255)
                                                                                                                                     ccTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACt                  >  2:6168545/1‑150 (MQ=255)
                                                                                                                                     ccTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACt                  >  2:7496089/1‑150 (MQ=255)
                                                                                                                                     ccTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACt                  >  2:4163963/1‑150 (MQ=255)
                                                                                                                                      cTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTc                 >  1:7235896/1‑150 (MQ=255)
                                                                                                                                        gggggggTTATTTACGTGTAGGAGGGGGTTAACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTcga               <  1:5171287/149‑1 (MQ=39)
                                                                                                                                         ggggggTGCGGTCCCTGTCGGCGCGGGTTTCCCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTcgac              <  2:4847087/150‑1 (MQ=38)
                                                                                                                                         gggCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTcgac              >  2:25993/1‑150 (MQ=255)
                                                                                                                                              ggtgtgttACTTGTTGGTGGGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATa          <  2:526622/143‑1 (MQ=255)
                                                                                                                                                   tGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATACGCTg     >  2:9086533/1‑150 (MQ=255)
                                                                                                                                                    gTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATACGCTgg    >  2:7366088/1‑150 (MQ=255)
                                                                                                                                                    gTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATACGCTgg    >  2:4513809/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
TCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATACGCTGGCG  >  NC_002947/4863006‑4863303

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