Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,624,5140AT41.2% 38.1 / 13.3 34I491F (ATT→TTT) ppkpolyphosphate kinase, component of RNA degradosome
Reads supporting (aligned to +/- strand):  ref base A (4/16);  new base T (14/0);  total (18/16)
Fisher's exact test for biased strand distribution p-value = 2.28e-06
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCACGTAAAGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTG  >  NC_000913/2624367‑2624664
                                                                                                                                                   |                                                                                                                                                      
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)
                                                                                                                                              aCTTTATTGAAAGCCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATaa       <  1:620612/151‑1 (MQ=255)
                                                                                                                                                   aTTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCtt   <  2:376528/150‑1 (MQ=255)
                                                                                                                                                   aTTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTg  >  2:343073/1‑151 (MQ=255)
                                                                                                                                                   |                                                                                                                                                      
TCACGTAAAGAAAACGGTGAAGTGGTGCGTTATGCACACATCGGGACCGGGAACTTTAACGAAAAAACCGCGCGTCTTTATACTGACTATTCGTTGCTGACCGCCGATGCGCGCATCACCAACGAAGTACGGCGGGTATTTAACTTTATTGAAAACCCATACCGTCCGGTGACATTTGATTATTTAATGGTATCGCCGCAAAACTCCCGCCGCCTATTGTATGAAATGGTGGACCGCGAGATCGCCAACGCGCAGCAAGGGCTGCCCAGTGGTATCACCCTGAAGCTAAATAACCTTG  >  NC_000913/2624367‑2624664

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