Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,866,181 A→C 5.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,866,1810AC 5.2% 882.3 / 23.8 285intergenic (‑138/+99)ligA/PP_4275DNA ligase/cell division protein ZipA
Reads supporting (aligned to +/- strand):  ref base A (175/88);  new base C (5/12);  total (185/100)
Fisher's exact test for biased strand distribution p-value = 3.20e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 286 total aligned reads displayed.

TTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATACGCTGGC  >  NC_002947/4866032‑4866329
                                                                                                                                                       |                                                                                                                                                    
 tCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTAc                                                                                                                                                     >  2:3623808/1‑150 (MQ=255)
   ggcggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCt                                                                                                                                                   <  1:8460827/150‑1 (MQ=255)
   ggcggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCt                                                                                                                                                   <  1:5139562/150‑1 (MQ=255)
    gcggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTg                                                                                                                                                  <  1:5238377/150‑1 (MQ=255)
      ggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTa                                                                                                                                                <  2:3333940/150‑1 (MQ=255)
            gggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCACCCCCAAGGGCGTCACAAGCCCGAGGGCGGTGATGTAACGGTGGGAGCg                                                                                                                                          >  1:3863238/1‑150 (MQ=255)
             ggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCgg                                                                                                                                         <  1:1350809/150‑1 (MQ=255)
             ggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCgg                                                                                                                                         <  2:3108592/150‑1 (MQ=255)
             ggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCgg                                                                                                                                         <  2:3442182/150‑1 (MQ=255)
             ggATTTGTCGTCTCTTGCAAAGCAAAAGAGCAGCGTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCgg                                                                                                                                         <  2:8228703/150‑1 (MQ=255)
              gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCggg                                                                                                                                        <  1:1141753/150‑1 (MQ=255)
              gATGTGTCTTCTCTTGCAAAGCAAACGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCggg                                                                                                                                        <  1:6416801/150‑1 (MQ=14)
               aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGt                                                                                                                                       <  1:1557428/150‑1 (MQ=14)
               aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGAGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGt                                                                                                                                       <  2:7849033/150‑1 (MQ=255)
                tttGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGtt                                                                                                                                      <  2:6629678/150‑1 (MQ=255)
                 ttGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGttt                                                                                                                                     <  1:2387181/150‑1 (MQ=255)
                 ttGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGAAACCCGCGCCCACAACCGCGAGACAAGCCCGAGGGCGGTGATGGACCGGGTGGAGCGGGttt                                                                                                                                     >  1:6001709/1‑150 (MQ=255)
                   gTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTAc                                                                                                                                   <  1:2277082/150‑1 (MQ=255)
                      ttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCACCCACAAGCGCGTCCCCAGCCCGAGGGCGGGGGTGTACCGGGAGGGGCGGGTTTTCCcg                                                                                                                                >  1:3101754/1‑150 (MQ=255)
                      ttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCACCCACAAGGGCGTCACATCCCCAAGGGCCGTGAAGTACCAGTAGGAGCCGGTTTACCcg                                                                                                                                >  1:1428613/1‑150 (MQ=255)
                       tctcTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTTAAAACCCTCCCACAACGGCGGCCAATGCCCGAGGGGGGGGATGGAACTGTATGAGCGGGGTTAacccc                                                                                                                               >  1:3726759/1‑148 (MQ=255)
                        ctctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGGAGGAGCGGGTTTACCcgcg                                                                                                                              >  2:7987282/1‑150 (MQ=255)
                         tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGa                                                                                                                             >  2:5647955/1‑150 (MQ=255)
                         tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGa                                                                                                                             >  2:2004013/1‑150 (MQ=255)
                         tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCGAGGGCGGGGCTGTACCGGTAGGAGCGGGTTTACCCGCGa                                                                                                                             >  2:2981526/1‑150 (MQ=255)
                         tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAAGCCCTAGGGCGGGGATGTACCTGTAGGAGCGGGTTTACCCGCGa                                                                                                                             >  1:5827596/1‑150 (MQ=255)
                         tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCACCCCCCACGGCGTAACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGa                                                                                                                             >  2:1239883/1‑150 (MQ=255)
                         tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCACCCACCAGGGCGTCCCATGCCCTAGGGCGGTGTTGTACCGGTAGGAGCGGGGTTAACCGCGa                                                                                                                             >  2:2360588/1‑150 (MQ=255)
                           ttGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAag                                                                                                                           >  2:3160730/1‑150 (MQ=255)
                           ttGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTAACCGCGAag                                                                                                                           >  2:7603638/1‑150 (MQ=255)
                           ttGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCAAAATGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGAag                                                                                                                           >  1:8676631/1‑150 (MQ=255)
                            tGCAACGCAAATGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGAGTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGGGGGTTTACCCGCGAaga                                                                                                                          <  2:5948622/150‑1 (MQ=255)
                             gcaaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAagag                                                                                                                         >  2:5131892/1‑150 (MQ=255)
                              caaagcaaaTTAGCTGCTTGGGGTGGTGGTTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACATGGGCGTCACTTGCCCTTGGGCGGTGTTGTACCTGTAGGAGCGGTTTTACCCGCGAAGAgg                                                                                                                        <  2:4064395/150‑1 (MQ=37)
                              caaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAgg                                                                                                                        >  2:8371717/1‑150 (MQ=255)
                                aagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGGCACATGCCCAAGGGCGGGGATGTACCGGGAGGAGCGGGTTTACCCGCGAAGAGGcc                                                                                                                      >  1:5723753/1‑150 (MQ=38)
                                 agcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACACGCTCCCACAACGGCGTCACAAGCCCGAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCg                                                                                                                     >  2:1653689/1‑150 (MQ=38)
                                  gcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAGCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCgg                                                                                                                    <  2:6980822/150‑1 (MQ=255)
                                   caaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGc                                                                                                                   >  2:6439894/1‑150 (MQ=255)
                                   caaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGAAACCCGCTCCCACCAGGGCGTCACATGCCCAAGGGCGGGGCTGTACCTGTAGGAGCGGGGTTAACCGCGAAGAGGCCGGc                                                                                                                   >  1:7931511/1‑150 (MQ=255)
                                      aaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTTAAACCCCGGCCCCAAGGGCGTCACAAGCCCGAGGGCGGGGGTGTAACGGTAGGAGCGGGGTTTCCCGCGGAGAGGCCGGcaca                                                                                                                >  1:6285375/1‑150 (MQ=255)
                                              cTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAGGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCGGCACAGGCAacga                                                                                                        <  2:3559328/150‑1 (MQ=255)
                                              cTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCACCCACAAAGGCCTCACAAGCCCAACCGCGGGGGTGGACAAGGAGGAGCGGGGTTACCCCCCAAAACGCCCGCACACGCAaccc                                                                                                        >  1:5597686/1‑148 (MQ=255)
                                               ttGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCCCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAacgac                                                                                                       >  2:3916666/1‑150 (MQ=255)
                                                       gtgCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCg                                                                                               >  2:4182303/1‑150 (MQ=255)
                                                         gCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCGCCCCCAAGGGCGTCACATGCCCTAGGGCGGGGCTGTACCGGGAGGAGCGGGGTTACCCGCGAAGAGGCCGGCCCAGGCAACCACGAAAACCGaa                                                                                             >  1:538323/1‑150 (MQ=38)
                                                                      ggggtCTGTTGACTGGGCCAGCCCCTTCGCGGGTAACCCCGCTCCCACACGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGc                                                                                <  2:4999784/147‑1 (MQ=38)
                                                                       gtgtCTGTTGAGTGGGCCCGCCTCTGCGCGGGAAACCCCGCCCCCCCCCGGGCGTCACAGGCCCTCGGGGGTTGATGTACCTTTCGGCGGGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGcc                                                                               <  2:8286794/150‑1 (MQ=255)
                                                                            tGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCGAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGcccgccg                                                                          >  2:490270/1‑146 (MQ=255)
                                                                             gTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAgg                                                                         >  2:3582513/1‑150 (MQ=255)
                                                                                gccggggCCCGCCTGTTCGCGGGTAACCCCGCTCCCCCACGGGCGGCACAGGCCCTTGGGGGGTGTTGTCCCTGTCGGCGGGGGTTTACCCGCGACGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGgctg                                                                      <  2:8458256/144‑1 (MQ=255)
                                                                                   tgtgCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAAGCCCTAGGGCGGTGATGTACCTGGAGGAGCGGGTTTACCCGCGAAGAGGCCGGCCCAGGCAACGACGATACCCGAACCCGCCAGCACGCCCACAGgctgcac                                                                   >  2:1131581/1‑148 (MQ=38)
                                                                                       cccgCCTCTTCGCGGGTAACCCCGCTCCCACACGGGCTTCACATGCCCTCGGGCGGTGATGTACATGTTGGCGCGGGTTAACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTc                                                               <  2:718824/147‑1 (MQ=37)
                                                                                       ccAGCCTCTTCGCGGGAAAACCCGCTCCCACACGGGCGTCACAGGCCCTCGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTc                                                               <  2:6072111/150‑1 (MQ=255)
                                                                                          gCCTCTTCGCGGGTAAAACCCCGCCCCCAAAGGCGTAAAAAGCCCCAGGGCGGGGATGTTACGGTAGGAGCGGGTTTAACCCCGCACAGGCCGGCACAGGCAACGACGATACCCCAAACGGCCCGCACGCCCGCAGGCTGCTCTTTCgga                                                            >  2:1573372/1‑148 (MQ=255)
                                                                                            ctctTCGCGGGTAAACCCCCTCCCACAAGGGCGTCACATGCCCAAGGGCCGGGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCGCAGGCTGCTCGgtcgtcag                                                          >  2:8615349/1‑148 (MQ=38)
                                                                                             tctTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATTTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTgg                                                         >  1:7187333/1‑150 (MQ=255)
                                                                                                gagcgGGTTTACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATc                                                      >  1:8379970/3‑150 (MQ=255)
                                                                                                 cgcgGTTAAACCCGCCCCCACACGGGCGCCACATGCCCTGGGGGGGTGATGTACCTGTAGGCGGGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCt                                                     <  2:4665209/150‑1 (MQ=38)
                                                                                                  gcgGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTg                                                    >  2:7384309/1‑150 (MQ=255)
                                                                                                    gggTAAACCCCGCCCCCACAGGGGCGCCCCTGGCCCTGGGGGGTTGTTTTACGTTTCGGGGGGGTTTAACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGt                                                   <  1:7108743/150‑1 (MQ=255)
                                                                                                        aaaCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgcg                                              >  1:8292957/1‑150 (MQ=255)
                                                                                                          aCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt                                            >  2:7483670/1‑150 (MQ=255)
                                                                                                          aCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt                                            >  1:4827875/1‑150 (MQ=255)
                                                                                                          aCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt                                            >  2:4145373/1‑150 (MQ=255)
                                                                                                             ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCTGGCACAGTCAaca                                                                                                         >  1:2619154/1‑86 (MQ=255)
                                                                                                             ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCTGGCACAGTCAaca                                                                                                         <  2:2619154/87‑2 (MQ=255)
                                                                                                             ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                                          >  2:4416727/1‑150 (MQ=255)
                                                                                                             ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                                          >  1:5434214/1‑150 (MQ=255)
                                                                                                             ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                                          >  1:6902267/1‑150 (MQ=255)
                                                                                                                cccccACACGGGCTTCACTTGCCCTGGGGGGGTTATGTACCTGTTGGCGGGGTTTTCCCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTgg                                       <  2:4798767/148‑1 (MQ=255)
                                                                                                                cTCCCACAAGGGCGTCACTTGCCCTTGGGCGGTGATGTACCTGTTGGTGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTgg                                       <  2:7338057/150‑1 (MQ=255)
                                                                                                                 tCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTggg                                      >  1:1692595/1‑150 (MQ=255)
                                                                                                                  cccACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGt                                     >  1:952692/1‑150 (MQ=255)
                                                                                                                   ccACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGtt                                    >  1:3972499/1‑150 (MQ=255)
                                                                                                                   ccACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGtt                                    >  1:2979389/1‑150 (MQ=255)
                                                                                                                    cacaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTa                                   >  2:1074164/1‑150 (MQ=255)
                                                                                                                     acaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAg                                  >  2:2186747/1‑150 (MQ=255)
                                                                                                                      accgggCGTCACTTGCCCTCGGGGGGTGTTGTACGTTTCGGCGGGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAgc                                 <  2:3988604/147‑1 (MQ=255)
                                                                                                                       aaGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAgcg                                >  1:7509710/1‑150 (MQ=255)
                                                                                                                       aaGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAgcg                                >  1:3564883/1‑150 (MQ=255)
                                                                                                                       aaGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAgcg                                >  2:5537618/1‑150 (MQ=255)
                                                                                                                       aaGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAgcg                                >  1:7509363/1‑150 (MQ=255)
                                                                                                                              tCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGgcg                         >  2:7437278/1‑150 (MQ=255)
                                                                                                                                acaTGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCt                       >  1:7319620/1‑150 (MQ=255)
                                                                                                                                acaTGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCt                       >  2:915592/1‑150 (MQ=255)
                                                                                                                                 ccggCCCTTGGGGGGTTTTGAACCTGTAGGCGGGGTTTTACCCGCGATGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTc                      <  1:3623808/147‑1 (MQ=255)
                                                                                                                                 caTGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTc                      >  1:7623656/1‑150 (MQ=255)
                                                                                                                                    gCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCTGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGTATCTGTCAGCGCTTCTGGGTTAGAGCACGGCGCTCGaa                   >  1:4542385/1‑150 (MQ=255)
                                                                                                                                    gCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGaa                   <  2:8106183/150‑1 (MQ=255)
                                                                                                                                      ccTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACt                 >  1:4506519/1‑150 (MQ=255)
                                                                                                                                         aGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTcga              >  2:3361028/1‑150 (MQ=255)
                                                                                                                                          gggCGGTTATGAACATTAAGGAGCGGGTTAACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTcgac             <  2:1431714/150‑1 (MQ=39)
                                                                                                                                          gggCGGTGTTGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTcgac             <  2:8655345/150‑1 (MQ=255)
                                                                                                                                            gCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTa                                   >  2:4953213/1‑126 (MQ=255)
                                                                                                                                               ggTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATa         >  2:7051454/1‑150 (MQ=255)
                                                                                                                                                  gATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTTGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATACGc      >  2:1850383/1‑150 (MQ=255)
                                                                                                                                                  gATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATACGc      >  2:5890214/1‑150 (MQ=255)
                                                                                                                                                   aTGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATACGCt     >  1:5774160/1‑150 (MQ=255)
                                                                                                                                                     gTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATACGCTgg   >  1:146166/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                    
TTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATACGCTGGC  >  NC_002947/4866032‑4866329

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: