Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,745,5090AT46.9% 21.3 / 18.8 33T236S (ACC→TCC) yiaN2,3‑diketo‑L‑gulonate TRAP transporter large permease protein
Reads supporting (aligned to +/- strand):  ref base A (6/11);  new base T (14/1);  total (21/12)
Fisher's exact test for biased strand distribution p-value = 9.54e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTA‑CCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCCTCGG  >  NC_000913/3745370‑3745656
                                                                                                                                           |                                                                                                                                                    
gtggtggCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTCTCCCGGCCCGGTTACCCCCA‑ccca                                                                                                                                           >  1:88493/1‑148 (MQ=255)
        aGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTA‑CCGAAGCCgggg                                                                                                                                   <  2:392744/150‑1 (MQ=255)
           ccAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTTGCTCCGCCCCGTTTTTCCCTACCCCAACCCGGGGc                                                                                                                                  >  1:354568/1‑149 (MQ=255)
             aGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTGCTCCGGCCTGTTTTCCCCTA‑CCGAAGACGGGGCCGtt                                                                                                                              >  2:647992/1‑150 (MQ=255)
              gTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTA‑CCGAAGCCGGGGCCGTtgc                                                                                                                            <  2:583461/151‑1 (MQ=255)
                    ttAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTGCTCCGGCCTGTTTTCCCCTA‑CCCAAGCCGGGGCCGTTGCTGCttt                                                                                                                      >  2:130116/1‑151 (MQ=255)
                     tAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTC‑CCGCAGCCGGGGCCGGTgctgct                                                                                                                        >  1:37331/1‑148 (MQ=255)
                      aaTTTCCCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTA‑CCGAAGCCGGGGCCGTTGCTGCTTTCt                                                                                                                    <  2:277964/151‑1 (MQ=255)
                       aTTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTTCCCCTA‑CCCAAGCCGGGGCCGTTGCTGCTTTCt                                                                                                                    >  2:429496/1‑150 (MQ=255)
                        tttACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTCGCTCCGGCCTGTTTTCCCCCA‑CCGAAGACGGGGCCGGTGCTGCTTTCt                                                                                                                    >  2:490329/1‑149 (MQ=255)
                              ccGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTA‑CCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTg                                                                                                            <  2:243382/151‑1 (MQ=255)
                                     cAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTTCCCCTT‑CCCAAGGCCGGGCCGTTGTTGCTTTCTTTGCGCTGt                                                                                                           >  2:215196/1‑145 (MQ=255)
                                       aaaaGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCTCTCCCGCCCGTTTACCCCTC‑CCTAAGGCCGGGGCGGTGCTGCTTTCTATGCTGTGTTTGt                                                                                                       >  1:53847/1‑147 (MQ=255)
                                       aaaaGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTA‑CCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGcc                                                                                                   <  1:219828/151‑1 (MQ=255)
                                          aGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTTGCTCCCGCCTGTTTTCCCCTA‑CCCAAACCGGGGGCGGTGTTGCTTTTTTTTGGCTGttttttgc                                                                                                    >  1:302250/1‑142 (MQ=255)
                                                       gAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTA‑CCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTg                                                                                   <  2:37330/151‑1 (MQ=255)
                                                        aaaTCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTGCCCCGGCCCGTTTTCCCCTA‑CCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTTTTTTCCGt                                                                                    >  2:365948/1‑149 (MQ=255)
                                                         aaTCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTA‑CCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGaa                                                                                 <  2:27126/151‑1 (MQ=255)
                                                          aTCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTGCTCCCGCCTGTTTTCCCCTA‑CCGAAGCCAGGGCCGTTGCTGCTTTTTATTCGCTGTTTTTCTCCCCAGTTAttt                                                                                         >  1:554323/1‑142 (MQ=255)
                                                           tCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTGCTTCGGCCCGTTTTCCCCTA‑CCCCAGACCGGGCCGGTGCTGCTTTCTTTGCGCCGTTTTTTGCCACAGTTTTTTAcc                                                                                      >  1:547055/1‑144 (MQ=255)
                                                                          tctcCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTTGCTCCGGCCTGTTTTCCCCTA‑ACCAAGCCGGGGCCGTTGCTGCTTTCTTTTCGCTGTTTTTTGCCACAGTTATTTACCGTGGAAAGACCTTcccccc                                                                   >  2:343390/1‑146 (MQ=255)
                                                                               ggTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGATGTCTACCCTCCCTGTTTCACCCTC‑CCCGATGGGGGGCGGGTTGTGCTTTTTTTTTGCGGTTTTTGGCCCCCTGTTTTTATCCTGAAATGTCCTTCGCCCCCCTTTg                                                             >  1:306036/1‑149 (MQ=255)
                                                                                   tCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCTCTCCGGCCTGTTTTCCCCTA‑CCGAAAACGGGGCCGTTGCTGCTTTTTTTGCGCTGTTTGTCGCCCCAGTTATTTTCCGTGAAATGACCTTTGCCACACCTTGGCATGt                                                       >  2:315480/1‑151 (MQ=255)
                                                                                     tGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGGTTTGCTCCGGCCCGTTTACCCCCA‑CCCAAGCCCGGGCCGGTGCTTCTTTCTTTGCGCTGTTTTTCGCCACCGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATTTa                                                      >  2:52393/1‑150 (MQ=255)
                                                                                      gggCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTA‑CCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGt                                                       <  1:475394/148‑1 (MQ=255)
                                                                                                tcttcCAGTCATCATTATTGGCGGGTTTTGCTCCCGCCTGGTTTCCCCCA‑CCCAAGCCGGGGGCGTTGCTGCTTTTTTTGCGCTGTTTTTCGCCCCAGTTATTTACCGTGAAAAGACCTTCGCCCCACTTTGGCATGTACTCAtttg                                              >  2:586232/1‑145 (MQ=255)
                                                                                                   tcCAGTCATCATTATTGGCGGGTTTTGGCCCGGCCTGTTTTCCCCTA‑CCCAAACCGGGGCCGGTGCTGCTTTCTATTCGTTGTTTGTCGGCCCAGTTATTTTCCGGGAAAAGAGCTTCGCCACACTTTTGCCTGTACTCATTTGTgcgg                                         >  1:295957/1‑149 (MQ=255)
                                                                                                          atcatTATTGACGGGTTTCGCTCCGGCCTGGTTTCCCCTG‑CCGAAGCCGGTGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTCCCGTGAAATGACCTGCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAccca                                <  1:121426/151‑1 (MQ=255)
                                                                                                                   ggCGGGTTTCGCTCCGGCCTGTTTACCCCTA‑CCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTgg                       <  1:143454/151‑1 (MQ=255)
                                                                                                                        gTTTCGCTCCGGCCTGTTTACCCCTA‑CCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATg                  <  1:215195/151‑1 (MQ=255)
                                                                                                                             gCTCCGGCCTGTTTACCCCTA‑CCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCt             >  1:548479/1‑151 (MQ=255)
                                                                                                                                       gTTTACCCCTA‑CCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCCt     <  1:602191/149‑1 (MQ=255)
                                                                                                                                        tttACCCCTA‑CCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCCTCgg  >  1:319402/1‑151 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
GTGGTGGCAGGCCAGTCGCCTTAATTTACCCCGCCAGCAAAAAGCAACGATGCAGGAAATCTGGCACTCTTTTGTCTCCGGTATCTGGGCGCTGTTTCTTCCAGTCATCATTATTGGCGGGTTTCGCTCCGGCCTGTTTACCCCTA‑CCGAAGCCGGGGCCGTTGCTGCTTTCTATGCGCTGTTTGTCGCCACAGTTATTTACCGTGAAATGACCTTCGCCACACTTTGGCATGTACTCATTGGCGCGGCGAAAACCACCTCAGTGGTTATGTTTCTGGTGGCCTCGG  >  NC_000913/3745370‑3745656

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