Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,863,150 A→T 8.5% intergenic (‑134/+103) 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,1500AT 8.5% 541.2 / 5.7 221intergenic (‑134/+103)ligA/PP_4275DNA ligase/cell division protein ZipA
Reads supporting (aligned to +/- strand):  ref base A (120/70);  new base T (8/12);  total (136/85)
Fisher's exact test for biased strand distribution p-value = 5.44e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 225 total aligned reads displayed.

GGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGAC  >  NC_002947/4863001‑4863291
                                                                                                                                                       |                                                                                                                                             
gggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGa                                                                                                                                               <  2:6452845/150‑1 (MQ=255)
 ggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGAt                                                                                                                                              <  2:411699/150‑1 (MQ=255)
 ggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGAt                                                                                                                                              <  1:2254990/150‑1 (MQ=255)
 ggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGAt                                                                                                                                              <  2:7042405/150‑1 (MQ=255)
   tttCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGt                                                                                                                                            <  2:6742548/150‑1 (MQ=255)
   tttCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAAGCCCTAGGGCGGGGATGt                                                                                                                                            >  1:4434709/1‑150 (MQ=255)
         cggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGt                                                                                                                                      <  1:2558835/150‑1 (MQ=255)
         cggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGt                                                                                                                                      <  1:6210159/150‑1 (MQ=255)
          ggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTa                                                                                                                                     <  1:308284/150‑1 (MQ=255)
           gTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAg                                                                                                                                    <  2:3377059/150‑1 (MQ=255)
           gTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAg                                                                                                                                    <  2:431514/150‑1 (MQ=255)
              aTGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCGAGGGCGGTGATGTACCGGgaggag                                                                                                                                 >  2:8855564/1‑150 (MQ=255)
               tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGc                                                                                                                                >  2:7156682/1‑150 (MQ=255)
               tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCGCCCACAAAGGCGACACATGCCCAAGGGCGGTGATGTACCGGGAGGAGc                                                                                                                                >  2:3041973/1‑150 (MQ=255)
               tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAAACCCACCCACAAAGGCGTCACAATCCCAAGGGCGGGGATGTAACAGTATGAGc                                                                                                                                >  1:4649316/1‑150 (MQ=255)
                 ggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTTAACCCCCGCCCACAAAGGCGGCACAAGCCCAAGGGCGGGGTTGTAACAGTAGGAGCgg                                                                                                                              >  1:8420756/1‑150 (MQ=255)
                 ggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCCCGCCCACCACGGCGTGACAAGCCCTACGGCGGGGATGTAACAGTATGAGCgg                                                                                                                              >  1:8421440/1‑150 (MQ=255)
                  gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCggg                                                                                                                             <  1:7132669/150‑1 (MQ=14)
                  gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCggg                                                                                                                             <  2:1857871/150‑1 (MQ=255)
                  gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCggg                                                                                                                             <  2:4539692/150‑1 (MQ=255)
                  gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCCTGAGGGCGTGCTGGCCGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCgg                                                                                                                              <  1:2046957/150‑1 (MQ=37)
                   aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGt                                                                                                                            <  2:6166106/150‑1 (MQ=255)
                      tGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTa                                                                                                                         <  1:3529670/150‑1 (MQ=255)
                       gTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTAc                                                                                                                        <  2:1049204/150‑1 (MQ=14)
                        tcttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTAcc                                                                                                                       <  1:388348/150‑1 (MQ=14)
                        tcttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTAcc                                                                                                                       >  2:7743231/1‑150 (MQ=255)
                         cttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTAccc                                                                                                                      <  2:2720291/150‑1 (MQ=255)
                         cttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCTCCCACAACGGCGTCACAATCCCAAGGGCGGTGATGTAACAGTAGGAGCGGGTTTAacc                                                                                                                      >  1:7037423/1‑150 (MQ=255)
                          ttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCcg                                                                                                                     <  1:9027265/150‑1 (MQ=14)
                          ttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCcg                                                                                                                     >  2:6822529/1‑150 (MQ=255)
                          ttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCcg                                                                                                                     >  1:3744103/1‑150 (MQ=255)
                             tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGa                                                                                                                  >  2:6566248/1‑150 (MQ=255)
                             tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGa                                                                                                                  >  2:2730543/1‑150 (MQ=255)
                             tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGa                                                                                                                  >  2:7349459/1‑150 (MQ=255)
                             tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCAAATGCCCTAGGGCGGTGATTTACCTGTAGGAGCGGGGTTAACCGCGa                                                                                                                  >  1:2075651/1‑150 (MQ=255)
                             tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCACAATGGGGTCACATGCCCTAGGGCGGTGATGTTGCTGTTGGTGGGGGTTTTCCCGCGa                                                                                                                  >  1:8903979/1‑150 (MQ=255)
                              ctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGaa                                                                                                                 <  2:2758813/150‑1 (MQ=11)
                               ttGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGAag                                                                                                                >  2:7644688/1‑150 (MQ=255)
                               ttGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGGCACATGCCCGATGGCGGTGCGGTACCGGGGGGCGCGGGTTTACCCGCGAag                                                                                                                >  1:4375027/1‑150 (MQ=38)
                                 gcaaagcaaaATAGCCGCTTGGGGTGGTGCGGGGGGTTCCCGGGGTCTGTTGACTGTGCCAGCCTCTGCGCGGGTAACCCCGCCCCCACACGGGCGTCACATGCCCTCGGGGGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAagag                                                                                                              <  2:1102353/150‑1 (MQ=255)
                                 gcaaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGGGATGTACCTGTAGGAGCGGGTTTACCCGCGAagag                                                                                                              >  2:3586815/1‑150 (MQ=255)
                                   aaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                                            >  2:6076299/1‑150 (MQ=255)
                                   aaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGGGATGTACCGGGAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                                            >  1:6695053/1‑150 (MQ=255)
                                       caaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGAAAACCCCTCCCACAAGGGCGGCCCAATCCCGAGGGCGGGGATGTACCGGTGGGAGCGGGTTTAACCGCGAAGAGGCCGGc                                                                                                        >  1:3575423/1‑150 (MQ=255)
                                       caaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGAAAACCCCGACCACAACGGCGTCACAAGCCCGAGGGGGGGGGGGTAACGGTGGGGGCGGGGTTAACCGCGAAGcggccggc                                                                                                        >  1:6216209/1‑150 (MQ=255)
                                        aaaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGGTTACCCGCGAAGAGGCCGGca                                                                                                       >  2:4956688/1‑150 (MQ=255)
                                         aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGGAAACCCCGCCCCCAAGGGCGGCCCCGGCCCCAGGGCGGGGGTGTTCCGGGAGGGGCGGGGTTACACCCGAAGAGggcggccc                                                                                                      >  1:121939/1‑148 (MQ=255)
                                         aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGAAACCCCCTCCCCAAAGGGCGGCACATGCCCAAGGGCGGGGATTTTCCTGTAGGAGGGGGGGTACCCGCCAAGACGCCGGcgc                                                                                                      >  1:5157373/1‑148 (MQ=255)
                                              gcagcTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAAGCCCGAGGGCGGTGATGTACCAGTAGGAGCGGGTTTACCCGCGAAGTGGCCGGCACAGGCa                                                                                                 >  2:822957/1‑150 (MQ=255)
                                               cagcTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCGCCCCAAAGGGCGTCACAATGCCTAGGGCGGTGTTGTTACGGTATGAGCGGGTTTAACCGCGAAGAGGGCGGCACAGGcca                                                                                                >  1:6350033/1‑148 (MQ=255)
                                                agcTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCGAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAac                                                                                               >  2:4740935/1‑150 (MQ=255)
                                                             gCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCCCGCCCCCAACGGCGGCACAAGCCCAAGGGCCGGGAAGTACCGGTAGGAGCGGGGTTACCCCCCAAGAGGCCCGCACACGCCACCACGAAACCCCaa                                                                                  >  1:5234041/1‑150 (MQ=255)
                                                                 tGGGGTTCCTGGTGTCTGTTGACTGTGCCAGCCTCTGCGCGGGAAACCCCGCGCCCACACGGGCGTCACTGGCCCTTGGGCGGTTTTGTACCTGTCGGAGCGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCg                                                                               <  2:2143447/150‑1 (MQ=37)
                                                                     gTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGCTACCCGAAACGGCCa                                                                           >  2:3216792/1‑150 (MQ=255)
                                                                         cAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGGCACATGCCCTAGGGCGGTGATGTAACGGGAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCAc                                                                       >  2:502774/1‑150 (MQ=255)
                                                                                            ccAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTc                                                    <  2:3526456/150‑1 (MQ=255)
                                                                                               gCCTCTTCGCGGGTAAACCCCCTCCCACAAGGGCGTCACATTCCCAAGGGCGGTGATGTACCGGTAGGAGCGGGTTTAACCCCGAAGAGGCCGGCACAGGCAAAGACCATAACCCAACCGGCCAGCACCCCCACACGCAGCTCTTTCGta                                                 >  1:7551236/1‑149 (MQ=25)
                                                                                                ccTCTTCGGGGTTAACCCCGCCCCCCCATGGGCGCAACATGCCCTTGGGCGTTGATGTACCTTTAGGAGGGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCt                                                <  2:2776787/150‑1 (MQ=255)
                                                                                                 ctctTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTAACCGCGAAGAGGCCGGCACAGGCAACGACGATAACCGAACCGGCCACCACGCCCTCAGGCTGCTCTTTCGTCTg                                               >  1:6072460/1‑150 (MQ=255)
                                                                                                   ctTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCggg                                                                                                                             <  1:5663207/70‑1 (MQ=255)
                                                                                                    ttCGGGGGCAAACCCTCTCCCACAAGGGCGTCACATGCCCTAGGGCGTTTATTTACCTGTAGTAGCGGGTTTACCCGCGAAGTGGCCGGCACATGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGAt                                            <  2:9135782/150‑1 (MQ=255)
                                                                                                    ggagcgGGTTTACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGAt                                            >  2:1269548/4‑150 (MQ=255)
                                                                                                      cgggggAAACCCCGCCCCCACACGGGCGACACAGGCCCTTGGGCGGTTATGTACCTGTCGGAGGGGGTTAACCCGCGAAGTGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCt                                          <  2:1585001/147‑1 (MQ=255)
                                                                                                       gcgGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTg                                         >  1:7222634/1‑150 (MQ=255)
                                                                                                        cgGGTAAACCCGCTCCCACATGGGCGTCACTTGCCCTTGGGCGGTGTTGTACCTGTTGGAGCGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGt                                        <  2:5047520/150‑1 (MQ=38)
                                                                                                         gggTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTc                                       >  1:7209063/1‑150 (MQ=255)
                                                                                                          gttAAACCCGCCCCCACCCGGGCGTCACTGGCCCTTGGGGGGTGTTGTACCTGTCGGAGGGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCa                                      <  2:2677690/148‑1 (MQ=255)
                                                                                                          ggTAACCCCGCTCCCACACGGGCGTCACATGCCCTAGGGGGGTGATGTACCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCa                                      <  2:8709801/150‑1 (MQ=255)
                                                                                                          ggTAACCCCGCCCCCACCCGGGGGTCACTTGCCCTTGGGGGGTTTTGTACGTGTTGGTGGGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCa                                      <  2:6193788/150‑1 (MQ=255)
                                                                                                            tAAACCCGCTCCCCCACGGGCGTCACAGGCCCTTGGGCGGTGATGTACCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgc                                    <  2:5157880/150‑1 (MQ=39)
                                                                                                             aaaCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgcg                                   >  2:3305422/1‑150 (MQ=255)
                                                                                                             aaaCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgcg                                   >  1:7314507/1‑150 (MQ=255)
                                                                                                             aaaCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgcg                                   >  2:4780932/1‑150 (MQ=255)
                                                                                                              aaCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgcgc                                  >  1:5267026/1‑150 (MQ=255)
                                                                                                               aCCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt                                 >  2:8760652/1‑150 (MQ=14)
                                                                                                               aCCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt                                 >  2:5609756/1‑150 (MQ=14)
                                                                                                               aCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt                                 >  2:578276/1‑150 (MQ=255)
                                                                                                               aCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt                                 >  2:8271483/1‑150 (MQ=255)
                                                                                                               aCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt                                 >  2:23550/1‑150 (MQ=255)
                                                                                                                 ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                               >  1:7803188/1‑150 (MQ=255)
                                                                                                                 ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                               >  2:6489889/1‑150 (MQ=255)
                                                                                                                   gCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTg                             >  1:935222/1‑150 (MQ=255)
                                                                                                                   gCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTg                             >  2:4243664/1‑150 (MQ=255)
                                                                                                                    cTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTgg                            >  2:9028512/1‑150 (MQ=255)
                                                                                                                     tCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTggg                           >  1:3386492/1‑150 (MQ=255)
                                                                                                                     tCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTggg                           >  1:6250696/1‑150 (MQ=255)
                                                                                                                        cacaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAgg                                                              <  1:2877530/112‑1 (MQ=255)
                                                                                                                        cacaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAgg                                                              >  2:2877530/1‑112 (MQ=255)
                                                                                                                         acaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAg                       >  1:7671517/1‑150 (MQ=255)
                                                                                                                         acaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAg                       >  1:8291531/1‑150 (MQ=255)
                                                                                                                            aGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAgcgc                    >  2:885040/1‑150 (MQ=255)
                                                                                                                             gggCGTCACATGCCCTAGGGGGGTTTTGTCCCTTTAGGAGCGGTTTAACCCGCGATGTGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCa                   <  1:6822529/150‑1 (MQ=255)
                                                                                                                                 gTCACATGCCCTAGGGCGGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACACGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGgc               >  1:1878261/1‑150 (MQ=38)
                                                                                                                                     caggCCCTTGGGGGGTTTTTTACCTGTTGGTGGGGTTTTACCCGCGAAGTGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTc           <  1:8241344/147‑1 (MQ=255)
                                                                                                                                      aTGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCg          >  2:1509588/1‑150 (MQ=255)
                                                                                                                                         cccTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAAc       >  1:6538641/1‑150 (MQ=255)
                                                                                                                                          ccTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACt      >  2:7496089/1‑150 (MQ=255)
                                                                                                                                          ccTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACt      >  2:6168601/1‑150 (MQ=255)
                                                                                                                                              ggggggTGCGGTCCCTGTCGGCGCGGGTTTCCCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTcgac  <  2:4847087/150‑1 (MQ=38)
                                                                                                                                              gggCGGTGATTTACCTGTAGGAGAGGGTTTACCCGAGAAGAGGCCGGCACAGGCAACTACGATACCAGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGACTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTccac  >  2:1966984/1‑150 (MQ=39)
                                                                                                                                                       |                                                                                                                                             
GGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGAC  >  NC_002947/4863001‑4863291

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: