Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913324,5800CG25.0% 50.7 / 10.7 24intergenic (+127/+116)ykgG/ykgHLutC family protein; putative electron transport chain YkgEFG component/putative inner membrane protein
Reads supporting (aligned to +/- strand):  ref base C (8/10);  new base G (3/3);  total (11/13)
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.73e-01
Rejected as consensus: Frequency below/above cutoff threshold.

TCTGATTATTGAGGATTGTTGAGAGAGATATGAAGTCTGGACATAAACCTGAAGAAGGCATACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAGGCACTCTCTATCAGAATGTTGTTTATTACTGCCCAGAGTTTTGTTG  >  NC_000913/324432‑324716
                                                                                                                                                    |                                                                                                                                        
tCTGATTATTGAGGATTGTTGAGAGAGATATGAAGTCTGGACATAAACCTGAAGAAGGCATACGAGGTTGAAAATAGAGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTCCTCTTCTGTAAATATCCTGGTATTAc                                                                                                                                          >  5:54106/1‑149 (MQ=255)
             gATTGTTGAGAGAGATATGAAGTCTGGACATAAACCTGAAGAAGGCATACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGata                                                                                                                             >  5:120347/1‑149 (MQ=255)
                            tatGAAGTCTGTACATAAACCTGAAGAAGGCATACGAGGGTGAAAATATCGTTTAGGACAACCTGACTTAACCCCGATCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACAAATGCCAGGATATTTACAGTGAAAAAg                                                                                                              >  4:97152/1‑149 (MQ=255)
                                          aTAAACCTGAAGAAGGCATACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGCATTGGTAATACCAGGATATTTAc                                                                                                                        >  2:106963/1‑125 (MQ=38)
                                                     gaagGCATACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGCATTGGTAATACCAGGATATTTa                                                                                                                         <  8:97786/113‑1 (MQ=38)
                                                           aTACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTc                                                                                                <  8:30234/132‑1 (MQ=255)
                                                           aTACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTc                                                                                                >  7:30234/1‑132 (MQ=255)
                                                               gAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGttcccttccc                                                                           <  6:120347/149‑1 (MQ=255)
                                                                           aGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACgg                                                                                      >  8:27118/1‑126 (MQ=255)
                                                                           aGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACgg                                                                                      <  7:27118/126‑1 (MQ=255)
                                                                           aGCGTTTAGGACAACCTGACTTAACCACGTTCCGTTATGGTTATGTTGCTCTTCTGTAAATATCCTTGTATTTGCAATGCCTGTATATTTAGAGTGTAATAGATTGGTATCCATTTGTGAATACGGGTTCCCTTCCCGGTCATTGTct                                                                <  6:54106/148‑2 (MQ=255)
                                                                               tttAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTg                                                                                                                    >  8:51646/1‑92 (MQ=255)
                                                                               tttAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTg                                                                                                                    <  7:51646/92‑1 (MQ=255)
                                                                                  aGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGc                                                                                                        <  2:141932/101‑1 (MQ=255)
                                                                                    gACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGCATTGGTAATACCAGGa                                                                                                                               <  8:98591/76‑1 (MQ=25)
                                                                                    gACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGCATTGGTAATACCAGGa                                                                                                                               >  7:98591/1‑76 (MQ=25)
                                                                                       aaCCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTAAAAACgg                                                                                      <  6:76239/114‑1 (MQ=255)
                                                                                            gACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTAc                                                                                                                                          <  6:60788/57‑1 (MQ=255)
                                                                                            gACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTAc                                                                                                                                          >  5:60788/1‑57 (MQ=255)
                                                                                            gACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAgg                                               >  4:18685/1‑148 (MQ=255)
                                                                                                                   ttATGCTGCTCTTCTGTTAATATCCTGGTATTACCAATTCCAGTATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAGGCACTCTCTATCAGAATGTTGTtta                       <  3:97152/149‑1 (MQ=255)
                                                                                                                           ctctTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAGGCACTCTCTATCAGAATGTTGTTTATTACTGc                <  3:18685/148‑1 (MQ=255)
                                                                                                                                   tAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAGGCACTCTCTATCAGAATGTTGTTTATTACTGCCCAGAGttt       >  6:156912/1‑149 (MQ=255)
                                                                                                                                      atatCCTGGCATTGGTAATACCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAGGCACTCTCTATCAGAATGTTGTTTATTACTGCCCAGAGTTttgt    >  1:20156/1‑149 (MQ=39)
                                                                                                                                         tCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAGGCACTCTCTATCAGAATGTTGTTTATTACTGCCCAGAGTTttgttg  <  5:156912/148‑1 (MQ=255)
                                                                                                                                                   aCCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCAGGTCATTGTCACTTCAGCTTTATACAGGCACTCTCTATCAGAATGTTGTTTATTACTGCCCAGAGTTttgttg  <  5:113950/138‑1 (MQ=255)
                                                                                                                                                   aCCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCAGGTCATTGTCACTTCAGCTTTATACAGGCACTCTCTATCAGAATGTTGTTTATTACTGCCCAGAGTTttgttg  >  6:113950/1‑138 (MQ=255)
                                                                                                                                                    |                                                                                                                                        
TCTGATTATTGAGGATTGTTGAGAGAGATATGAAGTCTGGACATAAACCTGAAGAAGGCATACGAGGTTGAAAATAGCGTTTAGGACAACCTGACTTAACCCCGTTCCGTCATGGTTATGCTGCTCTTCTGTAAATATCCTGGTATTACCAATGCCAGGATATTTACAGTGAAAAAGACAGGCATCCATTCCTGAAAACGGGTTCCCTTCCCGGTCATTGTCACTTCAGCTTTATACAGGCACTCTCTATCAGAATGTTGTTTATTACTGCCCAGAGTTTTGTTG  >  NC_000913/324432‑324716

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