Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,866,178 T→C 9.1% intergenic (‑135/+102) 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,1780TC 9.1% 709.5 / 42.6 241intergenic (‑135/+102)ligA/PP_4275DNA ligase/cell division protein ZipA
Reads supporting (aligned to +/- strand):  ref base T (144/73);  new base C (14/8);  total (160/81)
Fisher's exact test for biased strand distribution p-value = 8.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.22e-02

Only 100 of 245 total aligned reads displayed.

GGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATACG  >  NC_002947/4866029‑4866324
                                                                                                                                                      |                                                                                                                                                  
ggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGAt                                                                                                                                                    <  2:1222783/150‑1 (MQ=255)
ggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGAt                                                                                                                                                    <  2:5721244/150‑1 (MQ=255)
  tttCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCACCCCCAAAGGCGGGCCAAACCCGAGGGGGGTGgtgt                                                                                                                                                  >  1:1866877/1‑150 (MQ=255)
    tCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTAc                                                                                                                                                <  1:1298002/150‑1 (MQ=255)
    tCGGCGGTCATGGGATGTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTAc                                                                                                                                                <  2:4009338/150‑1 (MQ=255)
      ggcggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATTCCCAAGGGCGGTGATGTACCt                                                                                                                                              >  1:2154637/1‑150 (MQ=255)
         ggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTa                                                                                                                                           <  2:3560490/150‑1 (MQ=255)
             aTGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAAACCCGACCCCCAAGGCGTCACAAGCCCGAGGGCGGGGATGTAACGggaggat                                                                                                                                       >  1:627185/1‑149 (MQ=255)
              tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGc                                                                                                                                      >  2:6418779/1‑150 (MQ=255)
                ggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCgg                                                                                                                                    <  1:3414602/150‑1 (MQ=255)
                ggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCgg                                                                                                                                    <  1:4358259/150‑1 (MQ=255)
                 gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCggg                                                                                                                                   <  1:5138527/150‑1 (MQ=255)
                 gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATCCCCTAGGGCGGTGATGTACCTGTAGGAGCggg                                                                                                                                   <  1:3774603/150‑1 (MQ=255)
                  aTTTGTCTTCTCTTGCACAGCAACAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCGTCGCGGGTAAACCCGCTCCTCCAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGt                                                                                                                                  <  2:6367665/150‑1 (MQ=9)
                    ttGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGttt                                                                                                                                <  1:4285549/150‑1 (MQ=255)
                      gTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTAc                                                                                                                              <  1:5205944/150‑1 (MQ=14)
                      gTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTAc                                                                                                                              <  1:1654525/150‑1 (MQ=255)
                       tcttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTAcc                                                                                                                             <  1:3128315/150‑1 (MQ=255)
                        cttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTAccc                                                                                                                            >  2:3899175/1‑150 (MQ=255)
                          tctcTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCtgc                                                                                                                          >  2:6870507/1‑150 (MQ=255)
                          tctcTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCcgc                                                                                                                          >  2:5014591/1‑150 (MQ=255)
                           ctctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGAAAACCCCGCCCCCAAGGGCGTCACAAGCCCAAGGGCGGGGATGTTCCTGTAGGGGCGGGTTTAACcccg                                                                                                                         >  1:6348725/1‑150 (MQ=255)
                            tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTTAACCCCCGCCCACCAGGGCGGCACCAGCCCCAGGGCGGGGCGGTGACGGGAGGGGCGGGTTTTCCcgcgc                                                                                                                        >  1:3713551/1‑149 (MQ=255)
                            tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCCCGCCCCCAAGGGCGGCACCTGGCCGAGGGCGGGGGTGTAGCTGTTTGTGCGGGTTTTACCGCGa                                                                                                                        >  1:3616102/1‑150 (MQ=255)
                             gtggcacagcaaaAGCGCCGGTTGGGGTGGTGGGTGGGGTTCCCGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACACGGGCGTCACAGGCCCTCGGGCGGTGTTGTACCTGTAGGAGCGGGTTTACCCGCGaa                                                                                                                       <  2:29249/147‑1 (MQ=37)
                              ttGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCGCCCACAAAGGCGGCACATGCCCGAGGGCGGTGGTGTACCGGTAGGAGCGGGTTTAGCCGCGAag                                                                                                                      >  1:5410426/1‑150 (MQ=255)
                               tGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAaga                                                                                                                     <  2:5625618/150‑1 (MQ=255)
                                  aaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                                                  >  2:5363424/1‑150 (MQ=255)
                                   aagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGcc                                                                                                                 >  2:5530654/1‑150 (MQ=255)
                                    agcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTGAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGGGGTGTAGCTGTTGGGGCGGGTTTACCCGCGAAGAGGCCg                                                                                                                >  1:4161455/1‑150 (MQ=38)
                                     gaaaaaGAGCCGGTTGGGGTGGTGGTTGGGGTTCCTGTTGTCTGTTGACTGGCCCAGCCTCTTCGGGGGTAAACCCGCTCCCACATGGGCGTCAAATGCCCTTGGGGGGTGTTGTACGTGTAGGAGGGGGTTTACCCGCGAAGAGGCCgg                                                                                                               <  2:3222657/148‑1 (MQ=255)
                                       aaaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGca                                                                                                             >  2:6211917/1‑150 (MQ=255)
                                       aaaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGACCGGca                                                                                                             >  2:6212159/1‑150 (MQ=255)
                                       aaaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGGAGGAGCGGGTTTACCCGCGAAGAGGCCGGca                                                                                                             >  2:6757503/1‑150 (MQ=255)
                                        aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTTAACCCCCTCCCACAAGGGCGTGACATGCCCGAGGGCGGGGTTGTACCGGTTGGGGCGGGTTTTACCGCGGAGAGGGCGGcac                                                                                                            >  1:5397653/1‑150 (MQ=255)
                                           gagCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAgg                                                                                                         >  2:5579277/1‑150 (MQ=255)
                                                       gtggtgCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCGCCCACCACGGCGGCACATGCCCGAGGGCGGGGAGGTACCGGGAGGAGCGGGTTTAACCGCCCAGAGGCCGGCCCAGGCCACCACGAAAc                                                                                             >  1:6039025/1‑150 (MQ=255)
                                                         ggtgCGTGGGGTGCCTGGTGGCTGTTGACTGTGCCAGCCTCTGCGGGGGAAACCCCGCCCCCACACGGGCGCCACCGCCCCGCGGGGGGTTCTGGACGTTTCGGCGGGGGTTACCCCGGTATTAGGCCGGAACTGGCAACGACGATAccc                                                                                           <  2:5146726/150‑1 (MQ=255)
                                                                     tCCTGGTGTGTTTTGAGTGCCCCCCCCTCTGGGCGGTAAACCCCCCCCCCACCCGGGCGTCAATGGCCTTCGGGGGTTGTTGTACCTTTAGGCGGGTTTTAACCCGCTATGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGc                                                                               <  2:3053278/150‑1 (MQ=255)
                                                                          gtgtCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGcc                                                                          >  2:3275766/1‑150 (MQ=255)
                                                                          gtgggtgttgtCTGTGCCTGCCTCTTCGGGGTTAAACCCGCCCCCACACGGGCTTAACATGCCCTTGGGGGGTTTTTTACCTTTCGTCGGGGGTTTACCCGCTAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGcc                                                                          <  2:4736158/145‑1 (MQ=255)
                                                                          gggtgggTTCACGGTCCCCGCCGCTCCGGGGGAAACCCCGCCCCCACCCGGGCGCCACTGGCCCTCGGGCGGTGTTGAACGTGTCGGCGCGGTTTACCCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGcc                                                                          <  2:4736026/144‑1 (MQ=255)
                                                                              cttgtgccTGGGCCCGCATCTCCGGGGGAAAACCCGCCCCCACATGGGCTACACTTGCCGTCGGGGGGTTTTTTCCCTTTCGTAGGGTTTTAACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCa                                                                      <  2:927890/143‑1 (MQ=255)
                                                                                     cTGTGCCAGCCTCTTCGCGGGTAAACCCGCACCCACAAGGGCGTCACATGCCCAAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGGCGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCACCACCCCCGCAGgctgca                                                               >  2:1424398/1‑149 (MQ=38)
                                                                                         gCCAGCCTCTTCGCGGGTAAAACCCCTACCACAAGGGCGTAACAAGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCCGCACAGGCAACGACGATAACCCAACCGGCCAGCACCCCCTCAGGGTGCTAttt                                                           >  2:1743993/1‑150 (MQ=37)
                                                                                          cccgCCTCTTCGCGGGAAAACCGCCGCCCACTGGGGCGACACAGGCCCTTGGGGGGTGATGTACGTGTAGGCGGGGTTTTCCCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTc                                                          <  1:5003744/147‑1 (MQ=255)
                                                                                           cAGCCTCTTCGCGGGTAAACCCCCGCCCACAAGGGCGTCCCAAGCCCGAGGGCGGTGATGTACCTGGAGGAGCGGGTTTACCCCCGAAGAGGCCGGCCCAGGCAACGACGATACCCGCACCGGCCAGCACGCCCACAGGCTGCTCTTTCg                                                         >  2:3162485/1‑150 (MQ=38)
                                                                                           cAGCCTCTTCGCGGGTAAAACCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCg                                                         >  2:6624904/1‑150 (MQ=255)
                                                                                            aGCCTCTTCGCGGGTAGACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGt                                                        >  2:1179631/1‑150 (MQ=255)
                                                                                            aGCCTCTTCGCGGGTAAACCCGCGCCCACCAGGGCGGCACATGCCCAAGGGCGGTGGTGTACCGGTTGGCGCGGGTTTTCCCGCGAAGAGGCCGGCACCGGCAACGACGATACCCCAAACGGCCAGCCCGCCCGCAGGCGGCTCgttggg                                                        >  2:852081/1‑147 (MQ=255)
                                                                                              ccTCTTCGCGGTAAAACCCGCTCCCACAAGGGCTACAAATGCCCTAGGGCGGTTATGAACCTTTAGTAGCGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCt                                                      <  2:1589638/150‑1 (MQ=255)
                                                                                               ctctTGGGGGTTAAACCCGCTCCCACACGGGCGACACTGGCCCTCGGGGGTTTTTGTACCTTTTGGTGGGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTg                                                     <  2:3581869/150‑1 (MQ=255)
                                                                                                 gggagcgGGTTTACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGa                                                   >  2:3358718/5‑150 (MQ=255)
                                                                                                  ttCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGAt                                                  >  2:3762891/1‑150 (MQ=255)
                                                                                                   gagcgGGTTTACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCAAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATc                                                 >  1:2859675/3‑150 (MQ=255)
                                                                                                     gcgGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTg                                               >  1:413261/1‑150 (MQ=255)
                                                                                                      cgGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGt                                              >  2:5787395/1‑150 (MQ=255)
                                                                                                       gggTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTc                                             >  2:3741937/1‑150 (MQ=255)
                                                                                                           aaaCCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgcg                                         >  2:2471945/1‑150 (MQ=14)
                                                                                                           aaaCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgcg                                         >  2:2841403/1‑150 (MQ=255)
                                                                                                            aaCCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGAGCTGTCAgcgc                                        >  2:4801985/1‑150 (MQ=14)
                                                                                                            aaCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgcgc                                        >  1:757873/1‑150 (MQ=255)
                                                                                                            aaCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgcgc                                        >  1:6798921/1‑150 (MQ=255)
                                                                                                             aCCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGAAGAGGCTGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt                                       >  2:1308870/1‑150 (MQ=14)
                                                                                                             aCCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCACCGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt                                       >  2:2660498/1‑150 (MQ=14)
                                                                                                              cccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCtt                                      >  2:1978477/1‑150 (MQ=255)
                                                                                                              cccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCtt                                      >  2:2503323/1‑150 (MQ=255)
                                                                                                              cccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCtt                                      >  1:1515723/1‑150 (MQ=255)
                                                                                                              cccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCtt                                      >  1:2791204/1‑150 (MQ=255)
                                                                                                              cccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCtg                                      >  2:5988762/1‑149 (MQ=255)
                                                                                                               ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                                     >  2:6342135/1‑150 (MQ=255)
                                                                                                               ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                                     >  1:1183142/1‑150 (MQ=255)
                                                                                                               ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                                     >  1:2251608/1‑150 (MQ=255)
                                                                                                                 gCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTg                                   >  2:1133572/1‑150 (MQ=255)
                                                                                                                  cccccACATGGGCGTCACATGCCCTTGGGGGGTTTTGTACCTGTTGGAGGGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTgg                                  <  2:4694070/148‑1 (MQ=38)
                                                                                                                  cTCCCACATGGGCGTCACATGCCCTTGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTgg                                  <  2:3848474/150‑1 (MQ=255)
                                                                                                                  cTCCCACACGGGCGTCACATGCCCTTGGGCGGTTTTGTACCTGTAGGAGGGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTgg                                  <  2:4109456/150‑1 (MQ=39)
                                                                                                                   tCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTggg                                 >  1:1813709/1‑150 (MQ=255)
                                                                                                                   ccccACACGGGCGTCCCAGGCCCTCGGGGGGTGATGTCCCTGTAGGAGGGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTggg                                 <  2:5402766/149‑1 (MQ=38)
                                                                                                                    cccACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGt                                >  2:3336822/1‑150 (MQ=255)
                                                                                                                      cacaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTa                              >  2:6139540/1‑150 (MQ=255)
                                                                                                                       acaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATAccc                                                                                           >  1:6160648/1‑88 (MQ=255)
                                                                                                                        cctgggCGTCACCGGCCCTCGGGGGGTGATGTACCTGTCGGCGGGGGTTTACCCGCGAAGTGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAgc                            <  1:4525028/147‑1 (MQ=255)
                                                                                                                         aaGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAgcg                           >  2:1423145/1‑150 (MQ=255)
                                                                                                                         aaGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAgcg                           >  2:2255645/1‑150 (MQ=255)
                                                                                                                             gCGTCACTTGCCCTCGGGGGGTTTTGTACCTTTTGGTGGGGGTTACCCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACg                       <  2:1866978/150‑1 (MQ=255)
                                                                                                                               gTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGgc                     >  1:3370557/1‑150 (MQ=255)
                                                                                                                                   caTGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACTACGATACCCGAACCTGCCAGCACGCCCTCAGGCTGCTCTTACGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGAGCTc                 >  1:562986/1‑150 (MQ=255)
                                                                                                                                    aTGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCg                >  1:5892758/1‑150 (MQ=255)
                                                                                                                                     tGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGa               >  2:2959593/1‑150 (MQ=255)
                                                                                                                                     tGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGa               >  1:6179352/1‑150 (MQ=255)
                                                                                                                                       cccTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAAc             >  1:3063616/1‑150 (MQ=255)
                                                                                                                                       cccTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAAc             >  1:2724200/1‑150 (MQ=255)
                                                                                                                                       accTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAAc             <  2:6397405/149‑1 (MQ=255)
                                                                                                                                       accTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAAc             <  2:2429020/149‑1 (MQ=255)
                                                                                                                                        ccTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACt            >  2:1267560/1‑150 (MQ=255)
                                                                                                                                         cTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTc           >  1:4708821/1‑150 (MQ=255)
                                                                                                                                          tAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTcg          >  1:1220078/1‑150 (MQ=255)
                                                                                                                                           aGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTcga         >  1:3815223/1‑150 (MQ=255)
                                                                                                                                            ggggggTGATGTACCTTTTGGAGGGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTcgac        <  2:4834599/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                  
GGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGACGATACG  >  NC_002947/4866029‑4866324

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: