Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,866,156 G→A 8.7% intergenic (‑113/+124) 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,1560GA 8.7% 787.9 / 47.9 272intergenic (‑113/+124)ligA/PP_4275DNA ligase/cell division protein ZipA
Reads supporting (aligned to +/- strand):  ref base G (143/98);  new base A (7/16);  total (154/118)
Fisher's exact test for biased strand distribution p-value = 1.38e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.17e-02

Only 100 of 272 total aligned reads displayed.

GGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACG  >  NC_002947/4866008‑4866303
                                                                                                                                                    |                                                                                                                                                    
ggCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGt                                                                                                                                                     >  2:7291824/1‑150 (MQ=255)
 gCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAAACCCACCCCCAACGGCGTc                                                                                                                                                    >  1:6339936/1‑150 (MQ=255)
       aGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACAAGc                                                                                                                                              >  1:5057257/1‑150 (MQ=255)
       aGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCACCCACAACGGCGTCACATGc                                                                                                                                              >  1:1318917/1‑150 (MQ=255)
          tCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCt                                                                                                                                           <  1:1938231/150‑1 (MQ=255)
           cGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTa                                                                                                                                          <  2:744911/150‑1 (MQ=255)
           cGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTa                                                                                                                                          <  2:3584494/150‑1 (MQ=255)
            gAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAg                                                                                                                                         >  2:4719221/1‑150 (MQ=255)
             aGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAgg                                                                                                                                        <  1:4546889/150‑1 (MQ=255)
             aGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAgg                                                                                                                                        <  2:3255496/150‑1 (MQ=255)
             aGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAgg                                                                                                                                        >  2:2328439/1‑150 (MQ=255)
              ggATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATTACCTAggg                                                                                                                                       >  1:6661778/1‑150 (MQ=255)
              ggATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAggg                                                                                                                                       >  1:6059555/1‑150 (MQ=255)
                aTACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCg                                                                                                                                    <  1:313764/150‑1 (MQ=255)
                aTACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACCGCCCTAGGGcat                                                                                                                                    <  1:1704752/150‑3 (MQ=255)
                  aCGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGt                                                                                                                                  <  2:3045008/150‑1 (MQ=255)
                   cGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCGCCCCCAAGGGCGTCACATGCCCAAGGGCgggg                                                                                                                                 >  1:1031847/1‑148 (MQ=255)
                    gggTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGa                                                                                                                                <  2:4867584/150‑1 (MQ=255)
                      gTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATg                                                                                                                              <  2:5517907/150‑1 (MQ=255)
                       tttCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGt                                                                                                                             >  2:6294639/1‑150 (MQ=255)
                         tCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTAc                                                                                                                           <  2:4620988/150‑1 (MQ=255)
                         tCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTAc                                                                                                                           <  1:4000368/150‑1 (MQ=255)
                           ggcggTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTAccg                                                                                                                         >  2:4961966/1‑149 (MQ=255)
                              ggTCATGGGCTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTa                                                                                                                      <  2:5164120/150‑1 (MQ=255)
                               gTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTACCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAg                                                                                                                     >  2:3381757/1‑150 (MQ=255)
                                 cATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGa                                                                                                                   <  1:3610516/150‑1 (MQ=255)
                                 cATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGa                                                                                                                   <  1:3333289/150‑1 (MQ=255)
                                  aTGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATTCCCTAGGGCGGTGATGTACCGGTAGGAg                                                                                                                  >  2:1576580/1‑150 (MQ=255)
                                  aTGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAg                                                                                                                  >  2:2375733/1‑150 (MQ=255)
                                  aTGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCACAAGGGCGTCAAATTCCCTAGGGCGGTGATGTACCTGTAGGAg                                                                                                                  >  2:1870869/1‑150 (MQ=255)
                                   tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGccc                                                                                                                                            >  1:2662324/1‑124 (MQ=255)
                                   tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGc                                                                                                                 >  2:428755/1‑150 (MQ=255)
                                   tGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCGCCCACAAGGGCGTCACAAGCCCAAAGGCGGTGGTGTACCGGGTGGAGc                                                                                                                 >  1:1529225/1‑150 (MQ=255)
                                     ggATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCgg                                                                                                               <  1:651584/150‑1 (MQ=255)
                                      gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCggg                                                                                                              <  1:7047935/150‑1 (MQ=255)
                                      gATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCggg                                                                                                              <  2:2063718/150‑1 (MQ=255)
                                       aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGt                                                                                                             <  1:1669717/150‑1 (MQ=14)
                                       aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGt                                                                                                             <  2:6534468/150‑1 (MQ=255)
                                       aTTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCTCCCCCAACCGCCGCAAAAGCCCAACGGCCGTGATGTACCAGGAGGCGCgggg                                                                                                             >  1:7095497/1‑149 (MQ=255)
                                        tttGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCGGCCTCTTCGCGGGTAAGCCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGtt                                                                                                            <  2:5286537/150‑1 (MQ=14)
                                          tGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTa                                                                                                          <  2:4409553/150‑1 (MQ=255)
                                             cttctCTTGCAAAGCAACAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTAccc                                                                                                       <  2:6663488/150‑1 (MQ=255)
                                             cttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTAccc                                                                                                       <  2:1996587/150‑1 (MQ=14)
                                             cttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTAccc                                                                                                       <  2:6898864/150‑1 (MQ=14)
                                             cttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCGCCCACCAGGGCCTCACAAGCCCGAGGGCGGGGATGTACCGGTAGGAGCGGGTTTAacc                                                                                                       >  1:5292847/1‑150 (MQ=255)
                                              ttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCcg                                                                                                      <  1:5554980/150‑1 (MQ=14)
                                              ttctCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCcg                                                                                                      <  1:2375733/150‑1 (MQ=255)
                                               tctcTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCCCGCCCACAAGGGCCGCACAGGCCCGACGGCGGTGGTGTACCGGGAGGGGCGGGTTTACCcgc                                                                                                     >  1:2247005/1‑150 (MQ=255)
                                                  ctTGCAACGCAAACGAGCAGCTTGGGGTGGTGCGTGGGGTCCCAGGTGTCTGTTGACTGTGCCAGCCACTTCGCGGGTACACCCGCGCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGaa                                                                                                  <  2:6124606/150‑1 (MQ=39)
                                                   ttGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTa                                                                                                                      <  2:5258749/129‑1 (MQ=255)
                                                     gcaaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACCAGGGCGTCACAAGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAagag                                                                                               >  2:3730663/1‑150 (MQ=255)
                                                       aaagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGc                                                                                             >  2:4787678/1‑150 (MQ=255)
                                                        aagcaaaAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAACCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGTTTTACCCGCGAAGAGGcc                                                                                            <  2:1737234/150‑1 (MQ=255)
                                                            aaaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGca                                                                                        >  2:3875432/1‑150 (MQ=255)
                                                            aaaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCGCCCCCAAAGGCGGCACAAGCCCAACGGGGGTGGTGGACCGGGAGGGGCGGGGTTACCCGCGGAGAAGCCGGcc                                                                                        >  1:5718216/1‑149 (MQ=255)
                                                             aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGGGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGccc                                                                                       >  2:6374150/1‑148 (MQ=255)
                                                             aaaGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCCCTCCCACAAGGGCGTCACATGCCCTAGGGCGGGGATGTACCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGcac                                                                                       >  2:6373883/1‑150 (MQ=255)
                                                                gagCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAAACCGCCCCCCCCAGGGCGTCCCAAGCCCGAGGGCGGGGATGTTCCGGTTGGAGCGGGGTTTCCCGCGAAGAGGCCGGCACAgg                                                                                    >  1:5246754/1‑150 (MQ=255)
                                                                      gttGGGGTGGTGGGGGGGGTTCCCGGTGTCTGTTGATTGTGCCCGCCTCTGCGGGGGTAAACCCGCGCCCACCGGGGCGTCACAGGCCCTTGGGCGGTTTTGTACGTGTTGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAacga                                                                              <  2:2230924/149‑1 (MQ=25)
                                                                                                    tGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGTTGTACCGGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAg                                                >  2:1385352/1‑150 (MQ=255)
                                                                                                             tgCCAGCCTCTTCGCGGGTAAACCCCCTCCCACAAGGGCGTCCCATGCCCGAGGGCGGGGCGGTCCCGGTAGGAGCGGGTGTACCCGCGAAGAGGCCGGCACCGGCAACGCCGATACCCGAACCGGCCCGCACCCCCTCCGGCTGctcgt                                       >  2:5938982/1‑148 (MQ=255)
                                                                                                             ggCCAGCCTCTTGGCGGGTAACCCCGCGCCCACATGGGCGTCACTTGCCCTCGGGGGGTGTTGTACCTGTAGGCGCGGGTTAACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCtt                                       <  2:6983747/149‑1 (MQ=37)
                                                                                                               ccAGCCTCTTCGCGGGTAAACCCGCTCCCACACGGGCTTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTc                                     <  2:4265230/150‑1 (MQ=255)
                                                                                                               ccAGCCTCGCCGCGGGTAACCCCGCCCCCACACGGGCGTCACTGGCCCTAGGGCGGGCTTGTACCTGAAGGCGCGGGTTTACCCGCGAAGAGGCCGGCGCAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCtttt                                     <  2:5701485/150‑2 (MQ=37)
                                                                                                                    ctctTCGGGGGAAACCCCGCCCCCACACGGGGGGCACAGGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTg                                <  2:5874291/150‑1 (MQ=38)
                                                                                                                    ctctTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCt                                                                                                                         <  1:6973387/61‑1 (MQ=255)
                                                                                                                       ttCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGAt                             >  2:4167388/1‑150 (MQ=255)
                                                                                                                       ttCGCGGGTAAACCCGCTACCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGTGCGGGTTTAACCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTGGTCTGGAt                             >  1:1105768/1‑150 (MQ=255)
                                                                                                                         cgcgGGAAACCCCGCCCCCACAGGGGCGCCACCGGCCCTTGGGGGGTGTTGTACCTTTAGGCGCGGTTTCACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCt                           <  1:7410437/150‑1 (MQ=255)
                                                                                                                            gggTAACCCCGCTCCCACACGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTc                        <  2:1131271/150‑1 (MQ=255)
                                                                                                                            gggTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTc                        >  2:3634070/1‑150 (MQ=14)
                                                                                                                             ggTAAACCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCa                       >  2:4370130/1‑150 (MQ=14)
                                                                                                                               tAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgc                     >  1:6605232/1‑150 (MQ=255)
                                                                                                                                 aaCCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAgcgc                   >  2:2613010/1‑150 (MQ=14)
                                                                                                                                  aCCCGCTCCTACAGGTACATCACCGCCCTAGGGCATGTGACGCCCTTGTGGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAg                      <  1:4822915/146‑1 (MQ=14)
                                                                                                                                  aCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt                  >  2:2453194/1‑150 (MQ=255)
                                                                                                                                  aCCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCt                  >  1:5250920/1‑150 (MQ=255)
                                                                                                                                   cccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCtt                 >  2:5361791/1‑150 (MQ=255)
                                                                                                                                    ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCTGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTc                                     >  1:5637918/1‑129 (MQ=255)
                                                                                                                                    ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCTGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTc                                     <  2:5637918/129‑1 (MQ=255)
                                                                                                                                    ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                >  1:2841080/1‑150 (MQ=255)
                                                                                                                                    ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                >  1:4365148/1‑150 (MQ=255)
                                                                                                                                    ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                >  1:4222071/1‑150 (MQ=255)
                                                                                                                                    ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTc                >  2:2968577/1‑150 (MQ=255)
                                                                                                                                    ccGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCGTc                >  2:2816622/1‑150 (MQ=255)
                                                                                                                                     cGCTCCCCCACGGGCGTCACTTGCCCTCGGGGGGTGTTGTACCTGTCGGCGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCt               <  2:1031847/150‑1 (MQ=38)
                                                                                                                                      gCTCCCACACGGGCGTCACTTGCCGTTGGGGGGTGTTGTACCTGTTGGAGGGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTg              <  2:5292847/150‑1 (MQ=38)
                                                                                                                                          ccACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGtt          >  1:6124851/1‑150 (MQ=255)
                                                                                                                                          ccACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGtt          >  1:3204002/1‑150 (MQ=255)
                                                                                                                                          ccACAAGGGCGTCACATGCCCTAAGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCGCGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGtt          >  1:7307717/1‑150 (MQ=255)
                                                                                                                                           cccacgggCGTCACAGGCCCTTGGGGGGTGTGGTACCTGTCGGAGGGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTa         <  2:5010035/145‑1 (MQ=38)
                                                                                                                                            acaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAg        >  1:3759033/1‑150 (MQ=255)
                                                                                                                                            acaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAg        >  1:1862327/1‑150 (MQ=255)
                                                                                                                                            acaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAg        >  1:1425537/1‑150 (MQ=255)
                                                                                                                                            acaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAg        >  1:6573562/1‑150 (MQ=255)
                                                                                                                                            acaAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAg        >  1:1205274/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)
                                                                                                                                                gggCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCa    <  1:5495862/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGCACGCAGTTCGAGGATACGGGTTTCGGCGGTCATGGGATTTGTCTTCTCTTGCAAAGCAAAAGAGCAGCTTGGGGTGGTGCGTGGGGTTCCAGGTGTCTGTTGACTGTGCCAGCCTCTTCGCGGGTAAACCCGCTCCCACAAGGGCGTCACATGCCCTAGGGCGGTGATGTACCTGTAGGAGCGGGTTTACCCGCGAAGAGGCCGGCACAGGCAACGACGATACCCGAACCGGCCAGCACGCCCTCAGGCTGCTCTTTCGTCTGGATCTGTCAGCGCTTCTGGGTTAGCGCACG  >  NC_002947/4866008‑4866303

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: