Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,241,7810CT20.7% 56.2 / 10.8 29intergenic (‑37/+46)ygjV/uxaAImp‑YgjV family inner membrane protein/altronate hydrolase
Reads supporting (aligned to +/- strand):  ref base C (11/12);  new base T (3/3);  total (14/15)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.31e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TAGACCGAAAGCTGCTTTTTGAAGCGACGTTCGTCACGATTGAAAAATGTTGTGATACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGATCGCCAGTTCGCGGAAGTCGTTACGCTCGTTACAGGTTTGCTTACCGTTGGCAAACTCAACGATG  >  NC_000913/3241635‑3241917
                                                                                                                                                  |                                                                                                                                        
taGACCGAAAGCTGCTTTTTGAAGCGACGTTCGTCACGATTGAAAAATGTTGTGATACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGgcgc                                                                                                                                          <  4:42856/147‑1 (MQ=255)
   aCCGAAAGCTGCTTTTTGAAGCGACGTTCGTCACGATTGAAAAATGTTGTGATACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTca                                                                                                                                      >  3:103955/1‑148 (MQ=255)
                     aaGCGACGTTCGTCACGATTGAAAAATGTTGTGATACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGtttt                                                                                                                   <  4:103955/149‑1 (MQ=255)
                       gCGACGTTCGTCACGATTGAAAAATGTTGTGATACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGTTTCATAATGTGAAGCGCCGtttgct                                                                                                                 >  6:152633/1‑146 (MQ=255)
                           cGTTCGTCACGATTGAAAAATGTTGTGATACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGTTTCATAATGTGAAGCGCCGtttgcttttc                                                                                                             <  5:152633/149‑8 (MQ=38)
                             ttCGTCACGATTGAAAAATGTTGTGATACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTaac                                                                                                           <  8:2774/149‑1 (MQ=255)
                                       ttGAAAAATGTTGTGATACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTc                                                                                                 >  1:76210/1‑149 (MQ=255)
                                          aaaaaTGTTGTGATACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCTGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGa                                                                                               <  8:132489/148‑1 (MQ=255)
                                                 ttgtGATACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTtata                                                                                        <  2:76210/148‑1 (MQ=255)
                                                      aTACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGTTTCATAATGTGAAGCGCCGt                                                                                                                      >  1:23379/1‑113 (MQ=39)
                                                      aTACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGTTTCATAATGTGAAGCGCCGt                                                                                                                      <  2:23379/113‑1 (MQ=39)
                                                          cAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTAcgcc                                                                              <  7:95912/149‑1 (MQ=255)
                                                             tCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCt                                                                           >  6:158162/1‑149 (MQ=255)
                                                               aGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCtt                                                                          <  5:158162/148‑1 (MQ=255)
                                                                              cccACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCtt                                                                          >  8:148649/1‑133 (MQ=255)
                                                                                               ccaATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGAt                                                                  <  3:131734/124‑1 (MQ=255)
                                                                                               ccaATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGAt                                                                  >  4:131734/1‑124 (MQ=255)
                                                                                                 aaTACGCGCTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTa                                                                                                                                 <  1:57619/59‑1 (MQ=255)
                                                                                                 aaTACGCGCTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTa                                                                                                                                 >  2:57619/1‑59 (MQ=255)
                                                                                                   tACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGATCGCCAGTTCGCGGAAGTCGTTACGCTCGt                                     >  6:72961/1‑149 (MQ=255)
                                                                                                                         aGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACttt                                                                                           <  7:136717/73‑1 (MQ=255)
                                                                                                                         aGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACttt                                                                                           >  8:136717/1‑73 (MQ=255)
                                                                                                                                aTGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGATCGCCAGTTCGCGGAAGTCGTTACGCTc                                       <  4:89375/118‑1 (MQ=255)
                                                                                                                                    aaaaGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAAAAAGACATTACGACTTTAAAGCGTTACGCCGCTTTTGAAAAACGCCAGTAGGCGGAAGTCGTTACGATCGTTAAAGGTTTGCTTACCGTTGGCAAAAGCAACGa    >  6:162912/1‑149 (MQ=255)
                                                                                                                                             cGGCGTTTCATAATGTGAAGCGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTg                                                                       <  5:133923/73‑1 (MQ=25)
                                                                                                                                             cGGCGTTTCATAATGTGAAGCGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTg                                                                       >  6:133923/1‑73 (MQ=25)
                                                                                                                                              ggCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGATCGCCAGTTCGCGGAAGTCGTTACGCTCGTTACAGGTTTGCTTa                       <  5:31055/120‑1 (MQ=255)
                                                                                                                                              ggCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGATCGCCAGTTCGCGGAAGTCGTTACGCTCGTTACAGGTTTGCTTa                       >  6:31055/1‑120 (MQ=255)
                                                                                                                                              ggCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGATCGCCAGTTCGCGGAAGTCGTTACGCTCGTTACAGGTTTGCTTACCGTTGGCAAACTCAACGATg  >  4:32466/1‑141 (MQ=255)
                                                                                                                                              ggCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGATCGCCAGTTCGCGGAAGTCGTTACGCTCGTTACAGGTTTGCTTACCGTTGGCAAACTCAACGATg  <  3:32466/141‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                        
TAGACCGAAAGCTGCTTTTTGAAGCGACGTTCGTCACGATTGAAAAATGTTGTGATACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGATCGCCAGTTCGCGGAAGTCGTTACGCTCGTTACAGGTTTGCTTACCGTTGGCAAACTCAACGATG  >  NC_000913/3241635‑3241917

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: