Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,290,7160GA47.1% 37.9 / 13.9 34E159K (GAA→AAA) rssBPcnB‑degradosome interaction factor; response regulator
Reads supporting (aligned to +/- strand):  ref base G (16/2);  new base A (0/16);  total (16/18)
Fisher's exact test for biased strand distribution p-value = 6.99e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGAT  >  NC_000913/1290567‑1290853
                                                                                                                                                      |                                                                                                                                         
cAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGAAATTATTACAgg                                                                                                                                           >  2:609022/1‑150 (MQ=255)
   ttAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGAAATTATTACAGGAACt                                                                                                                                       >  1:211012/1‑151 (MQ=255)
        gATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGAAATTATTACAGGAACTACAAc                                                                                                                                  >  2:87075/1‑151 (MQ=255)
                 cGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTc                                                                                                                         >  1:681631/1‑151 (MQ=255)
                         cgAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTg                                                                                                                 >  1:99288/1‑151 (MQ=255)
                                       ccTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGGTTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCGGGTTCAGCAGGTGATTTCCCATTgt                                                                                                     >  2:377349/1‑148 (MQ=255)
                                        cTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGAGATTTCCCATTGCCgg                                                                                                  >  2:431174/1‑151 (MQ=255)
                                                           tttAATTCGCGCGTGGAGGAAGAGAAAAGTTTTTTTCGCGACGGGGATGAAATGTTTAAAACCCCTGCCGCGGCGGCAAAATTATTAAAGAAACTAAAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAAtt                                                                                 <  1:39804/149‑1 (MQ=255)
                                                               aTTCGCGCGTTGAGGAAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGc                                                                           >  1:152912/1‑151 (MQ=255)
                                                                   gcgcgTGGAGAAAGGGAAAGGGTTTTTTCGCGACGGGGATGAATTGGTTAATACCCCTCCCGCAGCGGAAAAATTATTACAGAAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCcgc                                                                        <  1:547587/150‑1 (MQ=255)
                                                                     ggtttGAGGAAGAGAAATGTTTTTTCCGGGGCGGGGATAAAATGGTTAATACCCCTGCCGCGGCGGAAAAATTATTCGAGAAACTAAAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTTTCGTCAATTGGTTGCcgc                                                                        <  2:363431/145‑1 (MQ=255)
                                                                     gcgTGGAGGAAGAGGAATGGTTTTTCGGCGGGTGGGATGCAATGTTTGATACCCCGCCCGCAGCGGCGAAATTATTACAGAAACTACAGCCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCcgc                                                                        <  1:430997/148‑1 (MQ=255)
                                                                            ggaagagaaaGGTTTTTTCCGCGGGGGGGATGAAATGGTTGAAACCCCTGCCGCGGGGGCAAAATTATCACAGAAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAAccc                                                              <  1:483535/151‑1 (MQ=255)
                                                                                    aaggtttttttCGCGACGGGGATGAAATGTTTAATACCCCTGCCGCAGCGGAAAATTTATTACAGAAACTAAAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGt                                                      <  1:609091/145‑1 (MQ=255)
                                                                                    aaggttttttcCGGGCCTGGGATGAATTGTTTAAAACCCCTCCCGCGGCGGAAAATTTATTACAGAAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTg                                                        <  1:741784/143‑1 (MQ=255)
                                                                                    aaagtttttttCGCGAGGGGGAGGAAATGGTTGAAACCCCTCCCGCAGCGGAAAAATTATTACAGAAACAACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTg                                                        <  1:848123/143‑1 (MQ=255)
                                                                                    aaaGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGt                                                      >  1:763515/1‑151 (MQ=255)
                                                                                    aaaGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAACTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGt                                                      >  1:21463/1‑151 (MQ=255)
                                                                                        gCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCtt                                                  >  1:410983/1‑151 (MQ=255)
                                                                                         cTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTg                                                 >  2:685432/1‑151 (MQ=255)
                                                                                            tttCGCGTTTGGGATGCATTGTTTGATACCCCTCCCGCGGCGGAAAAATTATTACAGGACCTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGat                                               <  1:862967/150‑1 (MQ=255)
                                                                                            tttCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGat                                               >  2:134969/1‑150 (MQ=255)
                                                                                                cgcgACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTg                                           >  1:787878/1‑151 (MQ=255)
                                                                                                             aattGGTTAAAACCCCTGCCGCGGGGGAAAAATTATTACAGAAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCg                                <  2:621719/146‑1 (MQ=255)
                                                                                                                 gtttaatcccccgcccGCGGCGGAAAATTTATAACAGAAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGaa                             <  2:734927/134‑1 (MQ=255)
                                                                                                                 ggttaacaccccTGCCGCGGCGGGAAATTTATTAAAGAAACTAAAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACttt                                  <  2:629434/134‑1 (MQ=255)
                                                                                                                  gttaaacccccTGCCGCAGCGGAAAAATTTTAACAGAAACTAAAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAAc                          <  2:667959/142‑1 (MQ=255)
                                                                                                                      aacccccTGCCGCGGGGGAAAAATTATTACAGAAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCt                     <  2:270163/147‑1 (MQ=255)
                                                                                                                      aTAACCCTGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCt                     >  2:411839/1‑151 (MQ=255)
                                                                                                                                cGCAGCGGCGGAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTAt           >  2:171156/1‑151 (MQ=255)
                                                                                                                                    gcggaaaatttattaAAGGACCCAAAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTg         <  2:195313/140‑1 (MQ=255)
                                                                                                                                        cgAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTg    >  2:286791/1‑150 (MQ=255)
                                                                                                                                          aaattttgaggggAACAACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGc        <  1:559853/133‑1 (MQ=255)
                                                                                                                                          aaaTTTTTACAGAAACTAAAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGAt  <  1:604376/150‑1 (MQ=255)
                                                                                                                                          aaaTTATTACAGAACCTCCAGCCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGa   <  2:433468/149‑1 (MQ=255)
                                                                                                                                          aaaTTATTACAGAAACTAAAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGAt  <  1:166764/150‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                         
CAGTTAAAGATCTGAATCGCTTGCGCGAGATGGTTTTTGCCTGTCTCTATCCCAGCATGTTTAATTCGCGCGTTGAGGAAGAGGAAAGGCTTTTTCGCGACTGGGATGCAATGGTTGATAACCCTGCCGCAGCGGCGAAATTATTACAGGAACTACAACCGCCGGTTCAGCAGGTGATTTCCCATTGCCGGGTTAATTATCGTCAATTGGTTGCCGCGGACAAACCCGGCCTGGTGCTTGATATTGCCGCACTTTCGGAAAACGATCTGGCATTTTATTGCCTTGAT  >  NC_000913/1290567‑1290853

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