Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,113,5480TC22.5% 74.7 / 19.0 40intergenic (+142/‑31)opgH/yceKOPG biosynthetic ACP‑dependent transmembrane UDP‑glucose beta‑1,2 glycosyltransferase; nutrient‑dependent cell size regulator, FtsZ assembly antagonist/outer membrane integrity lipoprotein
Reads supporting (aligned to +/- strand):  ref base T (15/16);  new base C (4/5);  total (19/21)
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.32e-01
Rejected as consensus: Frequency below/above cutoff threshold.

AACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGG  >  NC_000913/1113408‑1113696
                                                                                                                                            |                                                                                                                                                    
aaCGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTAt                                                                                                                                              >  1:5369/1‑149 (MQ=255)
aaCGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAt                                                                                                                                              <  4:79686/149‑1 (MQ=255)
                gcgcTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAagag                                                                                                                               <  2:152464/148‑1 (MQ=255)
                              aGGCCTAAATCGTTCGTGCAATTTAATGATTTTGAAAGACTTTGTGTGTCGGATAAGGCGTTAACGCCGCATCCGGAAAAAACAACGAGCCAATAAAAATACCGGCGGTATGCCTGTATTTTTTTACGAAAGAGGTATCAAATGAGtt                                                                                                                 <  8:31842/148‑1 (MQ=38)
                               ggCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAtttttatt                                                                                                                                       <  8:97126/125‑4 (MQ=255)
                               ggCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAtttttatt                                                                                                                                       >  7:97126/1‑122 (MQ=255)
                                   tACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATtg                                                                                                           >  5:80935/1‑149 (MQ=255)
                                         gTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGa                                                                                                     <  2:5369/149‑1 (MQ=255)
                                               gCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCAGGAAAAAACAACGAGCAAATAAAAATACCGGCGTTATGCAGGTATTTTTTTACGAAAGAGGTATCAAATGCTTTTAATTGTGGTCAGcatca                                                                                               >  7:115399/1‑149 (MQ=38)
                                                       ttGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGAc                                                                                        <  6:80935/148‑1 (MQ=255)
                                                                 aaGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATtg                                                                                                           >  1:67200/1‑119 (MQ=255)
                                                                 aaGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATtg                                                                                                           <  2:67200/119‑1 (MQ=255)
                                                                         gtgGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAtttt                                                                                                                                           <  5:40871/79‑1 (MQ=25)
                                                                         gtgGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAtttt                                                                                                                                           >  6:38401/1‑79 (MQ=25)
                                                                         gtgGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAtttt                                                                                                                                           >  6:40871/1‑79 (MQ=25)
                                                                         gtgGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCATAACGCCGGTAtttt                                                                                                                                           <  5:38401/79‑1 (MQ=25)
                                                                           gggTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCtgt                                                                        <  5:34039/144‑1 (MQ=255)
                                                                           gggTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCtgt                                                                        >  6:34039/1‑144 (MQ=255)
                                                                               cGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCat                                                                >  5:42268/1‑148 (MQ=255)
                                                                                 gATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGaaa                                                                                                                                  >  6:8734/1‑80 (MQ=255)
                                                                                  attaGGCTTTCACGCCGCTTCCGGCTAATAGTTCGAGCATATAAAATTTCCGGCGTTATGCCGGTATTTTTTTAGGAAAGATGTTTCTAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGTGTGGCTGTGGCAGCattatt                                                            <  8:115399/146‑1 (MQ=255)
                                                                                     aGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTaa                                                                                                              >  7:28224/1‑96 (MQ=255)
                                                                                     aGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTaa                                                                                                              <  8:28224/96‑1 (MQ=255)
                                                                                                  gcATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCtg                                                                         >  1:94303/1‑120 (MQ=255)
                                                                                                  gcATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCtg                                                                         <  2:94303/120‑1 (MQ=255)
                                                                                                     tCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATtgtg                                                                                                         <  6:136567/85‑1 (MQ=255)
                                                                                                               acaacGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGt                                                                                                                             >  1:80939/1‑55 (MQ=255)
                                                                                                               acaacGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGt                                                                                                                             <  2:80939/55‑1 (MQ=255)
                                                                                                                    gAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATg                                <  6:19796/143‑1 (MQ=255)
                                                                                                                    gAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATg                                >  5:19796/1‑143 (MQ=255)
                                                                                                                        cAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAata                      >  3:103479/1‑149 (MQ=255)
                                                                                                                              aaaTACCGGCATAACGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCa                                                     <  7:133481/112‑1 (MQ=255)
                                                                                                                              aaaTACCGGCATAACGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCa                                                     >  8:133481/1‑112 (MQ=255)
                                                                                                                                   ccGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGgt            <  6:42268/148‑1 (MQ=255)
                                                                                                                                    cGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGgtgt          >  6:134426/1‑149 (MQ=255)
                                                                                                                                         ttATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGt                                                                              <  6:161926/76‑1 (MQ=255)
                                                                                                                                         ttATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGt                                                                              >  5:161926/1‑76 (MQ=255)
                                                                                                                                           aTGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATc                 >  2:92920/1‑135 (MQ=255)
                                                                                                                                           aTGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATc                 <  1:92920/135‑1 (MQ=255)
                                                                                                                                           aTGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATg    <  4:103479/148‑1 (MQ=255)
                                                                                                                                            tGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATggg  <  5:134426/149‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                                    
AACGTAGTTGCCTGATGCGCTACGCTTATCAGGCCTACATCGTTCCTGCAATTTATTGATTTTGCAAGACTTTGTGGGTCGGATAAGGCGTTCACGCCGCATCCGGCAAAAACAACGAGCCAATAAAAATACCGGCGTTATGCCGGTATTTTTTTACGAAAGAGGTATCAAATGCGTTTAATTGTGGTGAGCATCATGGTGACCTTGCTGAGTGGCTGTGGCAGCATTATTAGTCGCACTATACCGGGGCAGGGGCATGGCAACCAATATTATCCCGGTGTGCAATGGG  >  NC_000913/1113408‑1113696

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: ‑