Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,241,7860CT20.0% 61.0 / 10.1 30intergenic (‑42/+41)ygjV/uxaAImp‑YgjV family inner membrane protein/altronate hydrolase
Reads supporting (aligned to +/- strand):  ref base C (11/13);  new base T (3/3);  total (14/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 = 9.83e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AAGCGACGTTCGTCACGATTGAAAAATGTTGTGATACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGATCGCCAGTTCGCGGAAGTCGTTACGCTCGTTACAGGTTTGCTTACCGTTGGCAAACTCAACGATG  >  NC_000913/3241656‑3241917
                                                                                                                                  |                                                                                                                                   
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                                                                                                                      <  2:23379/113‑1 (MQ=39)
                                 aTACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGTTTCATAATGTGAAGCGCCGt                                                                                                                      >  1:23379/1‑113 (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                                                                                           >  8:136717/1‑73 (MQ=255)
                                                                                                    aGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACttt                                                                                           <  7:136717/73‑1 (MQ=255)
                                                                                                           aTGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGATCGCCAGTTCGCGGAAGTCGTTACGCTc                                       <  4:89375/118‑1 (MQ=255)
                                                                                                               aaaaGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAAAAAGACATTACGACTTTAAAGCGTTACGCCGCTTTTGAAAAACGCCAGTAGGCGGAAGTCGTTACGATCGTTAAAGGTTTGCTTACCGTTGGCAAAAGCAACGa    >  6:162912/1‑149 (MQ=255)
                                                                                                                        cGGCGTTTCATAATGTGAAGCGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTg                                                                       >  6:133923/1‑73 (MQ=25)
                                                                                                                        cGGCGTTTCATAATGTGAAGCGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTg                                                                       <  5:133923/73‑1 (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)
                                                                                                                               tCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGATCGCCAGTTCGCGGAAGTCGTTACGCTCGTTACAGGt                              <  7:90649/107‑1 (MQ=255)
                                                                                                                               tCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGATCGCCAGTTCGCGGAAGTCGTTACGCTCGTTACAGGt                              >  8:90649/1‑107 (MQ=255)
                                                                                                                                  |                                                                                                                                   
AAGCGACGTTCGTCACGATTGAAAAATGTTGTGATACCAATCAGAAAGGCGATGACACCCACGCCCTGGGCCAGCCAATACGCGGTCATGATAAATCCTTAGCAGGTATGGAAAAGCAAACGGCGCTTCACATTATGAAACGCCGTTTTTTATTAACAACTCATTTCGACTTTATAGCGTTACGCCGCTTTTGAAGATCGCCAGTTCGCGGAAGTCGTTACGCTCGTTACAGGTTTGCTTACCGTTGGCAAACTCAACGATG  >  NC_000913/3241656‑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: