Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,756,5630AG20.7% 52.1 / 10.5 29intergenic (‑52/+113)aldB/yiaYaldehyde dehydrogenase B/L‑threonine dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (10/13);  new base G (3/3);  total (13/16)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.51e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CAATAAAGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCAC  >  NC_000913/3756424‑3756702
                                                                                                                                           |                                                                                                                                           
cAATAAAGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCa                                                                                                                                    <  8:77185/149‑1 (MQ=255)
  aTAAAGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGt                                                                                                                                  <  6:5957/149‑1 (MQ=255)
   tAAAGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCtt                                                                                                                                         >  8:83896/1‑141 (MQ=255)
   tAAAGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCtt                                                                                                                                         <  7:83896/141‑1 (MQ=255)
     aaGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCgg                                                                                                                               <  8:111773/149‑1 (MQ=255)
      aGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAGATTTCGCCAGTcacg                                                                                                                              <  8:130149/149‑4 (MQ=255)
              tAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCtt                                                                                                                      >  5:8869/1‑149 (MQ=255)
              tAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGCGACTGGCGAAGATTTCGCCAGTcacgt                                                                                                                             <  8:145203/142‑5 (MQ=255)
              tAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGCGACTGGCGAAGATTTCGCCAGTcacgt                                                                                                                             >  7:145203/1‑138 (MQ=255)
                           aCTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATc                                                                                                         <  1:81623/149‑1 (MQ=255)
                                     gAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTAc                                                                                               >  5:75684/1‑149 (MQ=255)
                                          cATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACa                                                                                                                          >  4:36288/1‑117 (MQ=255)
                                          cATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACa                                                                                                                          <  3:36288/117‑1 (MQ=255)
                                            tACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAg                                                                                        <  6:75684/149‑1 (MQ=255)
                                                          aaTCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTc                                                                                                                                        <  4:122850/87‑1 (MQ=255)
                                                                       gggggATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCgc                                                                                                                               <  4:37562/83‑2 (MQ=255)
                                                                          ggATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAAt                                                                                                   <  1:181829/108‑1 (MQ=255)
                                                                          ggATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAAt                                                                                                   >  2:181829/1‑108 (MQ=255)
                                                                             ttattGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTAc                                                                                  <  7:28700/122‑1 (MQ=255)
                                                                             ttattGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTAc                                                                                  >  8:28700/1‑122 (MQ=255)
                                                                                                                 taACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGAtt                           >  3:84111/1‑141 (MQ=255)
                                                                                                                 taACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGAtt                           <  4:84111/141‑1 (MQ=255)
                                                                                                                  aaCAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAgg                                                                                       >  8:143635/1‑80 (MQ=255)
                                                                                                                   aCAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTg                  <  6:8869/148‑1 (MQ=255)
                                                                                                                               gTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGc     >  5:130979/1‑149 (MQ=255)
                                                                                                                               gTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGc     >  2:114365/1‑149 (MQ=255)
                                                                                                                                  aCTGGCGAAGATTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGAtt                           >  1:127685/1‑124 (MQ=255)
                                                                                                                                  aCTGGCGAAGATTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGAtt                           <  2:127685/124‑1 (MQ=255)
                                                                                                                                   cTGGCGAAGATTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGACAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCcac  >  2:17778/1‑148 (MQ=255)
                                                                                                                                           |                                                                                                                                           
CAATAAAGTTGTCATAGCGGGCTTTTAACTTGAGGGGGAAACCATACTCGCCGGGCTTAATCTGTGCTGAAGGGGGATTATTGGTCATGATCGTCTCCTTGCGGTGTGAGGTATAACAAGGGTAGACGTGACTGGCGAAATCTTCGCCAGTCGGTAACAGCTTTACGACAGCTATCACGAATTTACGGGCAAGGGTTACAGAGATAGTTGAAAAAGCGGCTAACAATTTGCCAGCCGTTGTGGAAATGATGATTACATCGCTGCGCGATAAATCGCCAC  >  NC_000913/3756424‑3756702

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