Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,866,173 G→T 8.8% intergenic (‑130/+107) 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,1730GT 8.8% 826.2 / 50.2 260intergenic (‑130/+107)ligA/PP_4275DNA ligase/cell division protein ZipA
Reads supporting (aligned to +/- strand):  ref base G (145/90);  new base T (7/16);  total (154/106)
Fisher's exact test for biased strand distribution p-value = 6.50e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.43e-01

Only 100 of 265 total aligned reads displayed.

ATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGAC  >  NC_002947/4866024‑4866318
                                                                                                                                                       |                                                                                                                                                  
aTACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCg                                                                                                                                                    <  1:313764/150‑1 (MQ=255)
  aCGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGt                                                                                                                                                  <  2:3045008/150‑1 (MQ=255)
   cGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCGCCCCCAAGGGCGTCACATGCCCAAGGGCgggg                                                                                                                                                 >  1:1031847/1‑148 (MQ=255)
    gggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGa                                                                                                                                                <  2:1812187/150‑1 (MQ=255)
    gggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGa                                                                                                                                                >  2:260641/1‑150 (MQ=255)
    gggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGa                                                                                                                                                <  2:237225/150‑1 (MQ=255)
    gggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGa                                                                                                                                                <  2:4867584/150‑1 (MQ=255)
    gggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCGCCCACAAGGGCGTCACAAGCCCGATGGCGGTGa                                                                                                                                                >  1:2230924/1‑150 (MQ=255)
    gggTTTCGGCGGTCATGGGATTTGTCTTCTCGTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGa                                                                                                                                                <  2:5233491/150‑1 (MQ=255)
     ggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGAt                                                                                                                                               <  1:6599740/150‑1 (MQ=255)
     ggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCAAGGGCGGTGAt                                                                                                                                               >  2:2247172/1‑150 (MQ=255)
      gTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATg                                                                                                                                              <  2:5517907/150‑1 (MQ=255)
      gTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATg                                                                                                                                              <  1:4298343/150‑1 (MQ=255)
       tttCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGt                                                                                                                                             >  2:6294639/1‑150 (MQ=255)
       tttCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGt                                                                                                                                             <  2:331837/150‑1 (MQ=255)
         tCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTAc                                                                                                                                           >  2:5756205/1‑150 (MQ=255)
         tCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTAc                                                                                                                                           <  1:1670232/150‑1 (MQ=255)
         tCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTAc                                                                                                                                           <  2:4620988/150‑1 (MQ=255)
            gcggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTg                                                                                                                                        <  1:3666508/150‑1 (MQ=255)
               gTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTACCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAg                                                                                                                                     >  2:3381757/1‑150 (MQ=255)
                   tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGc                                                                                                                                 >  2:428755/1‑150 (MQ=255)
                   tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCGCCCACAAGGGCGTCACAAGCCCAAAGGCGGTGGTGTACCGGGTGGAGc                                                                                                                                 >  1:1529225/1‑150 (MQ=255)
                     ggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCgg                                                                                                                              <  1:651584/150‑1 (MQ=255)
                     ggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCgg                                                                                                                              <  2:363375/150‑1 (MQ=255)
                      gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCGGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCggg                                                                                                                             <  1:1146535/150‑1 (MQ=14)
                      gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCggg                                                                                                                             <  1:4680198/150‑1 (MQ=255)
                      gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCggg                                                                                                                             <  2:5286237/150‑1 (MQ=255)
                      gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCggg                                                                                                                             <  1:5673954/150‑1 (MQ=255)
                      gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCggg                                                                                                                             <  1:7047935/150‑1 (MQ=255)
                       aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGt                                                                                                                            <  2:6534468/150‑1 (MQ=255)
                       aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCTCCCCCAACCGCCGCAAAAGCCCAACGGCCGTGATGTACCAGGAGGCGCgggg                                                                                                                            >  1:7095497/1‑149 (MQ=255)
                        tttGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGtt                                                                                                                           <  2:4473114/150‑1 (MQ=14)
                          tGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTa                                                                                                                         <  2:4409553/150‑1 (MQ=255)
                             cttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTAccc                                                                                                                      <  2:1996587/150‑1 (MQ=14)
                              ttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATTCCCGAGGGCGGTGATGTACCTGGAGGAGCGGGTTTAaccc                                                                                                                     >  2:1582920/1‑149 (MQ=255)
                               tctcTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCcgc                                                                                                                    <  2:3656857/150‑1 (MQ=255)
                                 tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGa                                                                                                                  >  2:5327317/1‑150 (MQ=255)
                                 tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCAGCCACAAAGGCGGCACAAGCCCGAGGGGGGTGGTGGACCTGTAGGAGCGGGGTTAcccccca                                                                                                                  >  1:3325181/1‑146 (MQ=255)
                                 tctTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGTAACCCCCGCCCCCAACGGCGGCACCAGCCCAAGGGCGGTGGTGTTACTGTAGGGGCGGGGTTTAccgccc                                                                                                                   >  1:2604037/1‑146 (MQ=255)
                                   ttGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTa                                                                                                                                      >  1:5258749/1‑129 (MQ=255)
                                   ttGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCCCTCCCCCAAGGGCGTCACAATCCCAAGGGCGGTGCTGTTCCGGGGGGGGCGGGGTTACCCGCGacg                                                                                                                >  1:3909616/1‑148 (MQ=255)
                                     gcaaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCTCCCACAAGGGCGTCCCATGCCCGAGGGCGGGGTTGTACCGGGTGGAGCGGGTTTAACCGCGAagag                                                                                                              >  1:2263739/1‑150 (MQ=255)
                                        aagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCCCTCCCACCAGGGCGTCACCAGCCCTCGGGGGGTGAGGTACCTGGAGGGGCGGGTGTAGCCGCGCAGAGGcc                                                                                                           >  1:4131866/1‑150 (MQ=255)
                                            aaaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGGTAAACCCCGACCACAAAGGCGGCACATGCCCGAGGGCGGTGTTTTTCCGGTGGGTGCGGGTTTTCCCGCGGAGAGGCCGGcc                                                                                                       >  1:5718366/1‑149 (MQ=255)
                                            aaaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAAGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTTACCCCGAAGAGGCCGGca                                                                                                       >  1:5385928/1‑150 (MQ=255)
                                             aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGGGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGccc                                                                                                      >  2:6374150/1‑148 (MQ=255)
                                             aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCCCACACAACGGCGGCACAAGCCCGAGGGCGGGGGTGTTCCGGTAGGAGCGGGGTTTCCCCCGCAAAGGGCGGgac                                                                                                      >  1:6834177/1‑150 (MQ=255)
                                              aaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGcaca                                                                                                     >  2:4765107/1‑150 (MQ=255)
                                                gagCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCCCCCCCCAGGGCGTCCCAAGCCCGAGGGCGGGGATGTTCCGGTTGGAGCGGGGTTTCCCGCGAAGAGGCCGGCACAgg                                                                                                   >  1:5246754/1‑150 (MQ=255)
                                                                             tggTGTGTGTTGACTGTGCCCGCCTCTTCGGGGTTAACCCCGCTCCCACAGGGGCGTCACTTGCCCTTGGGGGTTGTTGTCCGTTTCGGTGGGGGTTAACCCGCGAAGCGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACg                                                                      <  2:6339936/149‑1 (MQ=255)
                                                                               gtgtCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTGGGAGCGGGTTTACCCGCGAAGAGGCCGGCACCGGCAACGACGATACCCGAACCGGCCAGCACGcc                                                                    >  2:213319/1‑150 (MQ=255)
                                                                                tgtCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCCCTCCCACAAGGGCGTCACAAGCCCGAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGGAGAGGCCGGCACAGGCAACGACGATAACCGAACCGGCCAGCACGccc                                                                   >  2:71337/1‑150 (MQ=38)
                                                                                  tCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTc                                                                 >  2:6590951/1‑150 (MQ=255)
                                                                                    tGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGTTGTACCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAg                                                               >  2:1385352/1‑150 (MQ=255)
                                                                                    tGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGGTGTACCTGTGGGAGCGGGTTTACCCGCGCAGAGGCCGGCACAGGCAACGACGATACCCCAACCGGCCAGCACGCCCACAg                                                               >  2:1903773/1‑150 (MQ=38)
                                                                                         ccTGTGCCCGCCTCTCCGCGGTTAAACCCGCTCCCACAAGGGCTTCACAGGCCCTCGGGCGGTGATGTACCTGTAGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGgctgc                                                          <  2:1103886/149‑1 (MQ=37)
                                                                                           tgtgCCAGCCTCTTCGCGGGTAAAACCCCTCCCACAAGGGCGTCACATGCCCTAGGGCGGGGATGTACCGGTTGGAGCGGGTTTACCCGCGCAGAGGCCGGCACCGGCACCGCCGATACCCGAACCGGCCACCACGCCCGCAGGCTGctc                                                        >  2:5692651/1‑150 (MQ=255)
                                                                                             tgCCAGCCTCTTCGCGGGTAAACCCCCTCCCACAAGGGCGTCCCATGCCCGAGGGCGGGGCGGTCCCGGTAGGAGCGGGTGTACCCGCGAAGAGGCCGGCACCGGCAACGCCGATACCCGAACCGGCCCGCACCCCCTCCGGCTGctcgt                                                      >  2:5938982/1‑148 (MQ=255)
                                                                                               ccAGCCTCTTCGCGGGTAAACCCGCTCCCACACGGGCTTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTc                                                    <  2:4265230/150‑1 (MQ=255)
                                                                                                    ctctTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCt                                                                                                                                         >  2:6973387/1‑61 (MQ=255)
                                                                                                     tctTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATTCCCAAGGGCGGTGATGTACCTGGAGGAGCGGGTTTACCCGCCGAGAGGGCGGCACAGGCCACGACGATACCCCAACCGGCCCGCACGCCCACACGCTGCTCTTTCGTCTgg                                              >  1:5189868/1‑150 (MQ=38)
                                                                                                       ttCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGAt                                            >  2:7336640/1‑150 (MQ=255)
                                                                                                       ttCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGAt                                            >  2:5289807/1‑150 (MQ=255)
                                                                                                         cgcgGGAAACCCCGCCCCCACAGGGGCGCCACCGGCCCTTGGGGGGTGTTGTACCTTTAGGCGCGGTTTCACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCt                                          <  1:7410437/150‑1 (MQ=255)
                                                                                                          gcgGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTg                                         >  1:2894367/1‑150 (MQ=255)
                                                                                                          gcgGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTg                                         >  2:6780000/1‑150 (MQ=255)
                                                                                                              ggTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCa                                      >  2:4370130/1‑150 (MQ=14)
                                                                                                                 aaaCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgcg                                   >  1:3500773/1‑150 (MQ=255)
                                                                                                                  aaCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgcgc                                  >  2:516368/1‑150 (MQ=255)
                                                                                                                   aCCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAg                                     <  1:4822915/146‑1 (MQ=14)
                                                                                                                   aCCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAg                                     >  2:4822915/1‑146 (MQ=14)
                                                                                                                   aCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCTGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt                                 >  2:2452898/1‑150 (MQ=255)
                                                                                                                   aCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt                                 >  2:4057855/1‑150 (MQ=255)
                                                                                                                    cccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCtt                                >  2:5361791/1‑150 (MQ=255)
                                                                                                                    cccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCtt                                >  2:142269/1‑150 (MQ=255)
                                                                                                                     ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCTGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTc                                                    <  2:5637918/129‑1 (MQ=255)
                                                                                                                     ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                               >  2:2968577/1‑150 (MQ=255)
                                                                                                                     ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCGTc                               >  2:2816622/1‑150 (MQ=255)
                                                                                                                      cGCTCCCCCACGGGCGTCACTTGCCCTCGGGGGGTGTTGTACCTGTCGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCt                              <  2:1031847/150‑1 (MQ=38)
                                                                                                                      cGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCt                              >  2:5874860/1‑150 (MQ=255)
                                                                                                                           ccACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGtt                         >  1:6124851/1‑150 (MQ=255)
                                                                                                                           ccACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGtt                         >  1:5614534/1‑150 (MQ=255)
                                                                                                                           ccACAAGGGCGTCACATGCCCTAAGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCGCGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGtt                         >  1:7307717/1‑150 (MQ=255)
                                                                                                                            cccacgggCGTCACAGGCCCTTGGGGGGTGTGGTACCTGTCGGAGGGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTa                        <  2:5010035/145‑1 (MQ=38)
                                                                                                                             acaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAg                       >  1:1425537/1‑150 (MQ=255)
                                                                                                                             acaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAg                       >  1:1862327/1‑150 (MQ=255)
                                                                                                                              ccggggCGCCACTTGCCCTTGGGGGGTTTTGTACGTGTAGGTGGGGTTTTCCCCGCGACGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAgc                      <  2:1318917/147‑1 (MQ=255)
                                                                                                                              cacgggCTTAACATGCCCTTGGGCGGTGATTTACCTGTAGGAGCGGTTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAgc                      <  2:3013270/147‑1 (MQ=39)
                                                                                                                              caAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgcgc                                  <  1:898772/138‑1 (MQ=255)
                                                                                                                                 gggCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCAGAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCa                   >  2:49056/1‑150 (MQ=255)
                                                                                                                                 gggCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCa                   <  1:5495862/150‑1 (MQ=255)
                                                                                                                                      tCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGgcg              >  2:374830/1‑150 (MQ=255)
                                                                                                                                        acaTGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCt            >  1:1252357/1‑150 (MQ=255)
                                                                                                                                         caTGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCt                                          >  2:6195709/1‑119 (MQ=255)
                                                                                                                                         caTGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCt                                          <  1:6195709/119‑1 (MQ=255)
                                                                                                                                         caTGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTc           <  1:5556070/150‑1 (MQ=255)
                                                                                                                                          aTGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCg          >  1:6638904/1‑150 (MQ=255)
                                                                                                                                           tGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGa         >  1:4109090/1‑150 (MQ=255)
                                                                                                                                             ccatatggCGGTGATGTACCTGAAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGCCCAGAACTCCCACACGCTGCTCTTTAGACTGGATCTGTCAGAGCATCTGGGTTAGCGAACGCCGCTCTAAc       >  2:4074431/7‑150 (MQ=25)
                                                                                                                                                 aGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTcga   >  1:1708926/1‑150 (MQ=255)
                                                                                                                                                       |                                                                                                                                                  
ATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACGGCGCTCGAACTCGAC  >  NC_002947/4866024‑4866318

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: