Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,866,177 A→T 10.2% 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,866,1770AT10.2% 479.0 / 5.2 196intergenic (‑134/+103)ligA/PP_4275DNA ligase/cell division protein ZipA
Reads supporting (aligned to +/- strand):  ref base A (99/63);  new base T (6/12);  total (120/76)
Fisher's exact test for biased strand distribution p-value = 4.12e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 201 total aligned reads displayed.

GGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATACG  >  NC_002947/4866028‑4866324
                                                                                                                                                      |                                                                                                                                                   
gggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGa                                                                                                                                                     <  2:1146826/150‑1 (MQ=255)
gggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGa                                                                                                                                                     <  1:891876/150‑1 (MQ=255)
gggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAACGGCGTCACATGCCCAAGGGCGGTGa                                                                                                                                                     >  1:5466854/1‑150 (MQ=255)
 ggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGAt                                                                                                                                                    <  2:1627895/150‑1 (MQ=255)
  gTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCACCCACAAGGGCGTCACAAGCCCAAGGGCGGTGATg                                                                                                                                                   >  1:5315070/1‑150 (MQ=255)
   tttCGGCGGTCATTGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGt                                                                                                                                                  <  2:4199898/150‑1 (MQ=255)
     tCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTAc                                                                                                                                                >  2:1819183/1‑150 (MQ=255)
     tCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTAc                                                                                                                                                <  1:4057820/150‑1 (MQ=255)
     tCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTAc                                                                                                                                                <  1:2378884/150‑1 (MQ=255)
       ggcggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCt                                                                                                                                              <  2:514387/150‑1 (MQ=255)
         cggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGt                                                                                                                                            <  2:2680921/150‑1 (MQ=255)
           gTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAg                                                                                                                                          <  1:1986970/150‑1 (MQ=255)
             cATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGa                                                                                                                                        <  2:3504055/150‑1 (MQ=255)
              aTGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCCCTCCCACAAGGGCCTCACAAGCCCGAGGGCGGTGATGTACCAGTAGGAg                                                                                                                                       >  2:1390648/1‑150 (MQ=255)
               tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTAc                                                                                                                                                <  2:4430008/140‑1 (MQ=255)
               tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTAc                                                                                                                                                >  1:4430008/1‑140 (MQ=255)
               tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGc                                                                                                                                      <  1:1881051/150‑1 (MQ=255)
                gggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCg                                                                                                                                     >  1:1684614/1‑150 (MQ=255)
                gggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCg                                                                                                                                     >  2:2464047/1‑150 (MQ=255)
                gggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACCGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCg                                                                                                                                     <  2:338599/150‑1 (MQ=255)
                 ggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCgg                                                                                                                                    <  2:1684614/150‑1 (MQ=255)
                 ggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCgg                                                                                                                                    <  2:3934669/150‑1 (MQ=255)
                 ggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCgg                                                                                                                                    <  2:236643/150‑1 (MQ=255)
                  gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCggg                                                                                                                                   <  2:373682/150‑1 (MQ=14)
                  gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCggg                                                                                                                                   <  1:5591075/150‑1 (MQ=14)
                  gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCggg                                                                                                                                   <  1:1857016/150‑1 (MQ=14)
                  gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCggg                                                                                                                                   <  2:1407105/150‑1 (MQ=255)
                   atgtGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGt                                                                                                                                  <  2:4374554/147‑1 (MQ=14)
                   aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGt                                                                                                                                  <  2:689513/150‑1 (MQ=14)
                    tttGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGtt                                                                                                                                 <  1:570961/150‑1 (MQ=14)
                    tttGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGtt                                                                                                                                 <  1:3730227/150‑1 (MQ=14)
                    tttGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGtt                                                                                                                                 <  2:2871169/150‑1 (MQ=14)
                     ttGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCACAAGGGCGTCACATGCCCGAGGGCGGTGATGTACCGGTAGGAGCGGGttt                                                                                                                                >  1:1338077/1‑150 (MQ=255)
                     ttGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCGCCCAAAAACGCGGCACAATCCCAAAGGCGGGGTGGTTCCAGTGGGAGCGGGttt                                                                                                                                >  1:4752425/1‑150 (MQ=255)
                      tGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTa                                                                                                                               <  1:2956695/150‑1 (MQ=14)
                         cttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTAccc                                                                                                                            <  2:1218971/150‑1 (MQ=255)
                             tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGt                                                                                                                                            >  2:4063315/1‑130 (MQ=255)
                             tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGa                                                                                                                        >  2:4080194/1‑150 (MQ=255)
                             tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGa                                                                                                                        >  2:3223702/1‑150 (MQ=255)
                             tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGa                                                                                                                        >  2:4678001/1‑150 (MQ=255)
                               ttGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTAACCGCGAag                                                                                                                      >  2:1997557/1‑150 (MQ=255)
                                tGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAACCTCTTCGCGGGTGAACCCGCTCCCCCAAGGGCGGCACATGCCCGAGGGCGGGGGGGTACCAGGGGGTGCGGGGTTAACCGCGCaga                                                                                                                     >  1:3107760/1‑150 (MQ=255)
                                 ccaaagcaacAGCCCCGGTTGGGGTGGTGCGTGGGGTTCCCGGTGTGTGTTGAGTGTGCCAGCCTCTTCGCGGGTAACCCCGCCCCCACACGGGCGTCACTTGCCGTTGGGCGGTGTTTTACCTTTTGTTGGGGGTTTACCCGCGAagag                                                                                                                    <  2:1860643/149‑1 (MQ=255)
                                   aaagcaaaAGAGCTGCTTGGGGTGGTGGGTGGGGTTCCTGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                                                  <  2:5510218/150‑1 (MQ=255)
                                       caaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGTAACCCCCACCCACCAGGGCGGCCCATGCCCGAGGGCGGGGGTGGTCCGGGTGGTGCGGGTTTACCCGCGAAGcggccggc                                                                                                              >  1:4810823/1‑150 (MQ=255)
                                           atagcagcTTGGGTTGGTGGGTGGTTTTCCTTGTGTATGTTCAATGTGCAAGCCTCTTCGCGGTAAAACCCGCTCCAAAAAGGGCGTAACAGGCCATAGGGCGTTTATTTACCTGTAGGAGCGGTTTAACCCGCGAAGAGGCCGGCACAg                                                                                                          <  2:4866954/148‑1 (MQ=255)
                                            gagCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGGGTCACATGGCCTAGGGCGGTGTTGTTGCTGTTGGAGCGGGTTTTGCCGCGTAGAGggcggcgccgg                                                                                                         >  2:4955296/1‑145 (MQ=37)
                                               cagcTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGATCGGGTTTACCCGCGAAGAGGCCGGCACAGGCaa                                                                                                      >  2:1413544/1‑150 (MQ=255)
                                                  cTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGTGTAACCCCGCTCCCACACGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTTGGAGGGGGTTTACCCGCGAAGAGGCCGGCACAGGCAacga                                                                                                   <  2:5208188/150‑1 (MQ=38)
                                                        gtggtgCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCCCTCCCAAAAGGGCGGCACAAGCCCGAGGGCGGGGGTGTACCGGTAGGAGCGGGTGTAACCCCGAAGAGGCCGGCACAGGCAAACACGATaa                                                                                             >  1:416685/1‑149 (MQ=37)
                                                         tggtgCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTCCGCGGTAAAACCCGCACCCACAAGGGCGTCACAGGCCCTCGGGCGGTGATGTGCCTGTTGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATAcc                                                                                            <  2:2758662/150‑1 (MQ=38)
                                                           gtgCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCCCGAAGAGGCCGGCCCAGGCAACGACGATACCCg                                                                                          >  2:1065121/1‑150 (MQ=255)
                                                               gTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCCCAAGGGCGTCACATGCCCGAGGGCGGTGATGTACCGGTAGGCGCGGGGTTAACCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAAcc                                                                                      >  2:2928641/1‑150 (MQ=39)
                                                                  gggTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGc                                                                                   >  2:3306814/1‑150 (MQ=255)
                                                                    gTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGGTGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACCACGATACCCGAACCGGCCa                                                                                 >  2:3355551/1‑150 (MQ=255)
                                                                       ccAGGTGTGTTTTGACTGTGCCCGCCTCTGCGGGGGTAAACCCGCTCCCACAAGGGCGTCACAGGCCCTCGGGGGTTGTTGTACCTTTAGGAGGGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCa                                                                              <  2:4550460/150‑1 (MQ=37)
                                                                               cggTTGACTGTGCCAGCCGCGCCGGGGGTAACCCCGCCCCCACACGGGGGTCACAGGCCCTCGGGGGGTGCTGTCCCTGTAGGAGCGGGTTAACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCa                                                                      <  2:3868911/148‑1 (MQ=25)
                                                                                tGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAAGCCCTAGGGCGGGGGTGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGTACCGGCCAGCACGCCCACAg                                                                     >  2:2183293/1‑150 (MQ=38)
                                                                                tGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCACCCACCAGGGCGTCAAAATCCCAAGGGCGGTGATGTACCGGTTGGAGCGGGGTTACCCGCGTAGAGGCCGGCACAGGCCACGACGATACCCGAACCCGCCAGCACCCCCTCAg                                                                     >  2:1186441/1‑150 (MQ=255)
                                                                                           cccgCCTCTGCGCGGGTAACCCCGCTCCCACACGGGCGTCACAGGCCCTCGGGCGGTGTTGTACCTGTCGGCGCGGGTTTCCCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTc                                                          <  2:3132252/147‑1 (MQ=37)
                                                                                            cAGCCTCTTCGCGGGTAAACCCGCTACCCCCAGGGCCGCACATGCCCAAGGGCGGTGATTTACCTGTAGGAGCGGGGTTACCCGCGAAGAGGCCGGCCCCGGCCACCACCATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCATTCg                                                         >  2:68914/1‑150 (MQ=37)
                                                                                                  gggagcgGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGa                                                   >  2:186662/5‑150 (MQ=255)
                                                                                                   ttCGCGGGTAAACCCGCTCCCCCACGGGCGTCACATGCCCTCGGGCGGTGATGCACCGGTCGGCGCGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGAt                                                  <  2:5356515/150‑1 (MQ=255)
                                                                                                      gcgGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTg                                               >  1:111269/1‑150 (MQ=255)
                                                                                                         ggTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCa                                            >  1:1670984/1‑150 (MQ=255)
                                                                                                          gTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGGTCTGTCAg                                           >  2:3172371/1‑150 (MQ=255)
                                                                                                           tAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgc                                          >  1:3601155/1‑150 (MQ=255)
                                                                                                             aaCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgcgc                                        >  2:5276980/1‑150 (MQ=255)
                                                                                                               cccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCtt                                      >  1:1406769/1‑150 (MQ=255)
                                                                                                               cccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCtt                                      >  1:5081590/1‑150 (MQ=255)
                                                                                                                ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                                     >  2:3616940/1‑150 (MQ=255)
                                                                                                                ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                                     >  2:745018/1‑150 (MQ=255)
                                                                                                                ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                                     >  1:2293392/1‑150 (MQ=255)
                                                                                                                ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                                     >  1:3781522/1‑150 (MQ=255)
                                                                                                                ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                                     >  1:4858759/1‑150 (MQ=255)
                                                                                                                 cGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCt                                    >  2:629035/1‑150 (MQ=255)
                                                                                                                 cGCCCCCCCATGGGCTTCACATGCCCTTGGGGGTTTTTGTACCTGTTGGAGGGGGTTAACCCGCGAAGTGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCt                                    <  1:5082440/150‑1 (MQ=255)
                                                                                                                   cTCCCACACGGGGGTAACATGCCATTGGGCGGTGATGTACCTGTAGGTGTTGGTTTACCCGCTAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTgg                                  <  2:3759875/150‑1 (MQ=255)
                                                                                                                   cTCCCACACGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTgg                                  <  2:4609547/150‑1 (MQ=255)
                                                                                                                      ccACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTCGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCATACGTCTGGATCTGTCAGCGCTTCTGGGtt                               >  2:1606317/1‑150 (MQ=255)
                                                                                                                      ccACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGtt                               >  2:2276717/1‑150 (MQ=255)
                                                                                                                      ccACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGtt                               >  1:1300492/1‑150 (MQ=255)
                                                                                                                       cacaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTa                              >  1:1212099/1‑150 (MQ=255)
                                                                                                                         caAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAgc                            <  2:5036374/150‑1 (MQ=255)
                                                                                                                          aaGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAgcg                           >  2:4012599/1‑150 (MQ=255)
                                                                                                                            gggCGTCACTGGCCCTTGGGGGGTGTTGTACTTGTTGGCGGGGGTTACCCCGCTACGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCa                         <  1:3163326/150‑1 (MQ=255)
                                                                                                                             ggCGTCACATGCCCTAGGGCGGTGATGTACCTGTGGGAGCGGGTTTACCCGCGAAGAgg                                                                                                                   <  1:3793121/59‑1 (MQ=255)
                                                                                                                             ggCGTCACATGCCCTAGGGCGGTGATGTACCTGTGGGAGCGGGTTTACCCGCGAAGAgg                                                                                                                   >  2:3793121/1‑59 (MQ=255)
                                                                                                                             ggCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCAc                        >  2:1449089/1‑150 (MQ=255)
                                                                                                                                ttCACAGGCCCTTGGGCGGTGTTGTACCTGTAGGCGCGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGgc                     <  2:2737229/149‑1 (MQ=39)
                                                                                                                                     aTGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCg                >  1:4290315/1‑150 (MQ=255)
                                                                                                                                     aTGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCg                >  1:195867/1‑150 (MQ=255)
                                                                                                                                       gCCCTTGGGCGGTTTTTTACCTGTTGGAGGGGTTTTACCCGCGAAGTGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGaa              <  1:2786143/150‑1 (MQ=255)
                                                                                                                                       gCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGaa              <  1:3498761/150‑1 (MQ=255)
                                                                                                                                         ccTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTg                                                     <  1:4090754/109‑1 (MQ=255)
                                                                                                                                         ccTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTg                                                     >  2:4090754/1‑109 (MQ=255)
                                                                                                                                           tAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTcg          >  1:3956784/1‑150 (MQ=255)
                                                                                                                                            gggggggTGTTGCACCTGTCGGTGCGGTTTACCCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTcga         <  1:4959530/149‑1 (MQ=255)
                                                                                                                                             gggCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTcgac        >  1:5032362/1‑150 (MQ=255)
                                                                                                                                                  ggTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATa    >  1:4551642/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                   
GGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATACG  >  NC_002947/4866028‑4866324

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: