Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,624,5080AT48.5% 21.9 / 18.5 33N489Y (AAC→TAC) ppkpolyphosphate kinase, component of RNA degradosome
Reads supporting (aligned to +/- strand):  ref base A (3/14);  new base T (16/0);  total (19/14)
Fisher's exact test for biased strand distribution p-value = 9.77e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CCTGATTTCACGTAAAGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAAT  >  NC_000913/2624360‑2624657
                                                                                                                                                    |                                                                                                                                                     
ccTGATTTCACGTAAAGAAAACGGGGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCCCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTaa                                                                                                                                                      <  1:168785/150‑1 (MQ=255)
   gATTTCACGTAAAGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTAAttt                                                                                                                                                  >  2:185032/1‑151 (MQ=255)
       tCACGTAAAGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTAtt                                                                                                                                               <  2:577766/150‑1 (MQ=255)
        cACGTAAAGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGCCTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTg                                                                                                                                              <  2:73070/150‑1 (MQ=255)
          cGTAAAGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGTCCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGaaa                                                                                                                                           <  2:602126/151‑1 (MQ=255)
               aGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGTTTTTACTTTTTTTTAAAccc                                                                                                                                       >  1:687960/1‑150 (MQ=255)
               aGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGCATTTTAATTTTTTTAAAAacc                                                                                                                                       >  2:208012/1‑150 (MQ=255)
                    aCGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTTATTTATTGAAAACCCACAcc                                                                                                                                  >  1:676187/1‑150 (MQ=255)
                         gAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTAACATTTTTTCAAACCCaca                                                                                                                                    >  1:339591/1‑141 (MQ=255)
                              ggtgCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACa                                                                                                                       <  2:673443/151‑1 (MQ=255)
                                   gTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTTTTGAAAAACCATAACCTCCgggggaattt                                                                                                                    >  2:543833/1‑141 (MQ=255)
                                     tATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGCAGCACGGCGGACGTTTTACTTTTTTCTAAACACCCACCGACCgggtgcattt                                                                                                                    >  1:155485/1‑139 (MQ=255)
                                       tGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTAATTTTTTGAAAAACCATACCCTCCGGGGACATTTTatt                                                                                                                >  2:538173/1‑149 (MQ=255)
                                          acacaTCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTACTTTTTTTAAAAACCATACCCTCCGGGGGCAttttttttttt                                                                                                            >  2:447739/1‑142 (MQ=255)
                                                 gggACCGGGAACTTTAACGAAAAAACCGCGCGTTTTTATACTGACTATTCGTTGCTTACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTa                                                                                                     <  1:437125/150‑1 (MQ=255)
                                                 gggACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTACTTTTTTGAAAAACCACACCGTCCGGGGGCATTTGATTATTTTATgg                                                                                                       >  2:537092/1‑148 (MQ=255)
                                                     ccGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTAATTTATTTAAAAACCATACCGTCCGGTGACATTTGTTTATTTAATGGTATtgcc                                                                                                >  2:328868/1‑151 (MQ=255)
                                                      cGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTTACTTTTTTTAAAACCCATACCGTCCGGGGGCATTTGATTTTTTTATGGTATcgccc                                                                                               >  1:197220/1‑150 (MQ=255)
                                                       gggAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTAACTTTTTTTAAAAACCATACCCTCCCGTGGCATTTGTTTATTTAATGGTATtgccg                                                                                               >  1:684735/1‑150 (MQ=255)
                                                             tttAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTAATTTTTTTAAAACCCATACCGTCCGGGGACATTTGATTTTTTAATGGTATCGCCGCaaaa                                                                                          >  2:645375/1‑149 (MQ=255)
                                                                    aaaaacCCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCccgcc                                                                                  <  1:344187/150‑1 (MQ=255)
                                                                    aaaaacCCCCGCGTTTTTATCGGGCCTTTTCGTTGTTGCCCGCCGATGCGCGCATCACCAACGAAGGACGGCGGGTATTTAACTTTATTTAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCccgccg                                                                                 <  1:468315/151‑1 (MQ=255)
                                                                             cgcgTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGt                                                                         <  2:576905/150‑1 (MQ=255)
                                                                                     tataCTGACTATTCGTTTCTTACCGCCGATGCGCGCATCACCAACGAAGTACGGCGTGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATTAAAtg                                                                 <  1:528180/150‑1 (MQ=255)
                                                                                              tATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTGTTTTACTTTATTTAAAACCCCTACCCCCCGGGGGCATTTTGTTTTTTTATTGTATTACCCCAAAACTCCCCCCGCCCATTGTATTAAAtggggggacg                                                         >  1:108038/1‑145 (MQ=255)
                                                                                                  cGTTGCTGGCCGCCTATGCGCGCATCACCACCGAAGTACGGCGGGTATTTAACTTAAACGAAAACCCATACCGTCCGGTGACGTTTGATTATTTAATCGGATCGACGCACAACTCCCGCCGCCTATAGTATGAAATGGTGGACCGCgag                                                     <  1:189391/149‑1 (MQ=255)
                                                                                                  cGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACTGCGGTTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGTTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCgaga                                                    <  1:332047/150‑1 (MQ=255)
                                                                                                     tGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTTACTTTTTTGAAAACACAGACCGTCCGGTGACAATTGATTATTTTATAGTATCGCCGCAAAACTCCCCCCGCCTTTTGTATGAAATGGGGGAGCGCGGGATc                                                  >  1:159034/1‑149 (MQ=255)
                                                                                                     tGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGGAGTTTACTTTTTTTAAAAACCATACCGTCCGGTGGCATTTGATTTTTTAAAGGTAATGCCGCACAAACCCCGCCCCCTATTTTTTGAAAATGTGGGCCGCGGGAGCg                                                 >  1:415528/1‑150 (MQ=255)
                                                                                                       cTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTTAATTTATTTAAAAACCATACCCGCCGGGGGCATTTTATTATTTTATGGGATCGCCGCCAAAAACCCCCCGCCCTTTTTTTGAAAAGGGGGGACGCGAGATCGccc                                              >  2:170833/1‑150 (MQ=255)
                                                                                                         gACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCa                                              <  2:336953/149‑1 (MQ=255)
                                                                                                              ccgATGCGCGCATCACCAACGAAGTACGGCGGGGATTTTAATTTTTTGAAAAACCCTACCGTCCGGGGGCATTTTATTATTTTATGGTATCGCCGCAAAACTCCCGCCGCCCCTTGTATGAAATGGTGGACCGCGGGATCGCCCACccgc                                        >  2:390722/1‑150 (MQ=255)
                                                                                                                   gcgcgcATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAg                                  <  2:276092/151‑1 (MQ=255)
                                                                                                                       gcATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCt                              <  1:536825/151‑1 (MQ=255)
                                                                                                                                                    aaCTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAAt  <  2:524528/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CCTGATTTCACGTAAAGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAAT  >  NC_000913/2624360‑2624657

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