Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,043,0300AG53.1% 27.3 / 43.1 32N206D (AAT→GAT) srkACpx stress response Thr/Ser protein kinase; MazF antagonist protein
Reads supporting (aligned to +/- strand):  ref base A (13/2);  new base G (17/0);  total (30/2)
Fisher's exact test for biased strand distribution p-value = 2.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-06

CACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGA  >  NC_000913/4042881‑4043175
                                                                                                                                                     |                                                                                                                                                 
cACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGaa                                                                                                                                                  <  2:870836/151‑1 (MQ=255)
                         acTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGGCTGCCACCCCGGGGATATTTTCTGGCCCGGTGGGCCAAt                                                                                                                          >  2:58874/1‑150 (MQ=255)
                             aTACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGGGGCGACATTGAGACTGCCCCGCCGGGGATATTCTCTGgcgggcggggcctcagt                                                                                                                        >  2:328375/1‑134 (MQ=255)
                               aCCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCTACATGGGGACTGCCACGCCGGGGATATTCTCTGGCGCGATGGTCCAAtgtttgtt                                                                                                                   >  1:297669/1‑151 (MQ=255)
                                   tCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGAATTTACCGTTCTGCGGCTACATGGAGGCTGCCACGCCCGGGAGATTCTCTGGGGCGCTGGGCCAAtggttgtt                                                                                                                   >  2:302486/1‑147 (MQ=255)
                                    ccGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTTTCTGGCGCGATGGTCCAATGTTTGTTGATCt                                                                                                              >  1:52964/1‑151 (MQ=255)
                                            aaaGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCCTTCTGCGGCTACATGGGGACTGCCACGCCGGGGATATTTTTTGGCGCGATGGGCCAATGTTTGTTGATCTGgatg                                                                                                         >  1:751787/1‑148 (MQ=255)
                                            aaaGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTTTTCGGGCACATGGAGACTGCCACGCCGGGGATATTCTCTGGCGCGGTGGTCCCATGTTTTTTGATCTTggtgat                                                                                                       >  2:15227/1‑150 (MQ=255)
                                            aaaGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCCGTCTGCGGCTACATGGAGACTTCCCCGCCGGGGATATTTTCTTGCGCGATGGTCCCATTTTTTTTGAt                                                                                                                >  2:477630/1‑141 (MQ=255)
                                                       ccTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCCACATGGAGACTGCCACGCCGGGGATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGt                                                                                           >  1:814462/1‑151 (MQ=255)
                                                        cTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTAACGTTCTGCGGCTACCTGGAGACTGCCACGCCGGGGATATTTTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCCCGTAAAgg                                                                                            >  2:661369/1‑149 (MQ=255)
                                                             aGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGATTTCACCGTTCTTCGGCTACATGGAGACTGCCCCGCCGGGGATATTTTCTGGCGCGATGGTCCAATGTTTGTTTGTCTGGATGATGCACGGAATGGTCCAGcc                                                                                     >  2:210599/1‑151 (MQ=255)
                                                              gCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGGATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATTGTCCAGcc                                                                                     >  1:333466/1‑150 (MQ=255)
                                                                  cAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGGAGATTTTACCGTTCTGCGGCTACATGGAGACTGCCCCGCCGGGGATATTTTCTGGCGCGGTGGGCCCATGTTTTTTGTTCTGGATGATGCACGTAATGGTCCAGCCGTTc                                                                                 >  2:473678/1‑150 (MQ=255)
                                                                     aTGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTTTGCGGCTACATGGAGACTGCCACGCCGGGAATATTTTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGGTGATGGACGTAATGGTCCAGCCGTTCAgg                                                                              >  1:855189/1‑150 (MQ=255)
                                                                     aTGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTTACCGTTCTGCGGCTACATGGGGACTGCCACGCCGGGGATATTTTCTGGCGCGGTGGTCCCATGTTTGTTGATCTGGGTGATGCCCGTAATGGTCCAGCCGTTCAgg                                                                              >  1:25148/1‑150 (MQ=255)
                                                                           tGATTGCCGCCGTTACAGCACACTGGCGGGGAGATTTTACCGTTCTGCGGCTACATGGAGAATGCCCCGCCGGGGATATTCTCTTGCGCGGGGGGCCCATGTTTGTTGTTCTGGGTGATGCACGTAATTGTCCAGCCGGTCAGGGTTtgt                                                                        >  2:772573/1‑150 (MQ=255)
                                                                             aTTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGGATATTCTCTTGCGCGATGGGCCAAAGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGGTTTTTgg                                                                      >  1:180183/1‑150 (MQ=255)
                                                                              ttGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCCTTCTGCGGCCACATGGAGACTGCCCCGCCGGGGATATTTTCTGGCGCGATGGGCCAATGTTTGTTGGTCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTgg                                                                      >  2:115857/1‑149 (MQ=255)
                                                                                     cgTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCCGCTACATGGAGACTGCCACGCCGGGGATATTTTCTGGCGCGATGGTCCCATGTTTGTTGATCTGGGTGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTc                                                             >  1:590342/1‑151 (MQ=255)
                                                                                                      gggAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCg                                             >  2:812278/1‑150 (MQ=255)
                                                                                                                  ccGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGc                                 >  1:786403/1‑150 (MQ=255)
                                                                                                                         gCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGaa                         <  1:116510/151‑1 (MQ=255)
                                                                                                                            gCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACt                      <  1:97738/151‑1 (MQ=255)
                                                                                                                               aCATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACtatt                   >  1:125972/1‑151 (MQ=255)
                                                                                                                                       aCTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGc           >  2:572088/1‑151 (MQ=255)
                                                                                                                                       aCTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGc           >  1:787939/1‑151 (MQ=255)
                                                                                                                                        cTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCt          >  2:207884/1‑151 (MQ=255)
                                                                                                                                         tGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCtt         >  1:217601/1‑151 (MQ=255)
                                                                                                                                              cGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATga     >  1:595637/1‑150 (MQ=255)
                                                                                                                                              cGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCCAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATga     >  1:196366/1‑150 (MQ=255)
                                                                                                                                                ccGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATga     >  1:401678/1‑148 (MQ=255)
                                                                                                                                                 cGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATgaaga  >  1:461859/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CACGCAAGCTGTTTGAGGACGCTACACTGATACCTTCCGGGTTGAAAGCGGCATTCCTGAAAGCGACAGATGAGCTGATTGCCGCCGTTACAGCACACTGGCGGGAAGATTTCACCGTTCTGCGGCTACATGGAGACTGCCACGCCGGGAATATTCTCTGGCGCGATGGTCCAATGTTTGTTGATCTGGATGATGCACGTAATGGTCCAGCCGTTCAGGATTTGTGGATGTTGCTCAATGGCGATAAAGCCGAGCAGCGGATGCAACTGGAAACTATTATTGAAGCTTATGAAGA  >  NC_000913/4042881‑4043175

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑