Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913753,1250TA38.5% 36.1 / 11.5 26intergenic (‑330/+60)ybgD/gltAputative fimbrial protein YbgD/citrate synthase
Reads supporting (aligned to +/- strand):  ref base T (14/2);  new base A (0/10);  total (14/12)
Fisher's exact test for biased strand distribution p-value = 1.24e-05
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

GGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATAC  >  NC_000913/752986‑753261
                                                                                                                                           |                                                                                                                                        
gggTTATAAATGTGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATATTTACAATTTTTTGTGTTTAGGCatat                                                                                                                                          >  2:775753/1‑140 (MQ=255)
           gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAATTCAAcc                                                                                                                               >  2:758915/1‑140 (MQ=255)
           gCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAcc                                                                                                                               >  1:232103/1‑140 (MQ=255)
            cGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAccc                                                                                                                              >  1:884903/1‑140 (MQ=255)
            cGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAAccc                                                                                                                              <  1:824129/140‑1 (MQ=255)
              aCTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGATGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGc                                                                                                                            >  1:630565/1‑140 (MQ=255)
              aCTACCATGAAGTTTTTAATTCAAAGCATTGGGATGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGACTGATAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCACCCCGc                                                                                                                            >  2:222984/1‑140 (MQ=255)
                           ttttAATTAAAATTTTTGGTTTGTTAATATTTTGAGTTTTTTTTTTTTTTTTAAATTTTTATATAGGTTTTTAAATGGATTTTTTTTTTACAATTTTTTGTTTTTAGGCATAAAAAAATCACCCCGCCATATGAAcggcg                                                                                                               <  2:732701/140‑1 (MQ=255)
                            tttAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAAcggcgg                                                                                                              >  2:559900/1‑140 (MQ=255)
                                  taaaaGTTTTGGTTTGCTAATATTTTGAGTTTTTTTTTTTTTTTTAAATATTTATATAGGTTTTTTAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTaaa                                                                                                        <  2:757173/138‑1 (MQ=255)
                                    aaaGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTAGATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAat                                                                                                      >  2:119709/1‑140 (MQ=255)
                                      aGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAatat                                                                                                    >  2:409017/1‑140 (MQ=255)
                                          ttGGGTTGTTGATAATTTGTGTTTTTTTTTTTTTTTTAAATATTTATATAGTTTTTTTAATGGATTTTTTTTTTACATTTTTTTGTTTTTAGAAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTAc                                                                                                <  2:30612/140‑1 (MQ=255)
                                                            gAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGTTTTTTTTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAt                                                                              <  1:938928/140‑1 (MQ=255)
                                                             aGCTTTTTTATTTTTTTTAAATATCTATATAGGTCTTTTAATGGATTTTATTTTTACATTTTTTTGTTTTTAGGCATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCAtt                                                                             <  2:232103/140‑1 (MQ=255)
                                                                 gTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACg                                                                         >  1:367541/1‑140 (MQ=255)
                                                                  ttCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGc                                                                        >  1:873234/1‑140 (MQ=255)
                                                                       tttttttaAAATATCTATATAGTTTTGTTAATGGTTTTTATTTTTACAATTTTTTGTGTTAAGAAAAAAAAAACCCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGat                                                                   <  2:367541/137‑1 (MQ=255)
                                                                            tttAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGc                                                              >  1:435687/1‑140 (MQ=255)
                                                                                          aaGGTTTTTTAATGGTTTTTTTTTTAAAATTTTTTTTTTTTTAGAAATAAAAAAACCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGt                                                <  2:873234/139‑1 (MQ=255)
                                                                                             gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGACATCGCTTTTAAAGTCGCGtttt                                             >  2:302469/1‑140 (MQ=255)
                                                                                             gTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAAAATTTAAAACTTAGCAATCAACCCTTAACGCTTGCTATCGCTTTTAAAGTAGCGtttt                                             >  2:379907/1‑140 (MQ=255)
                                                                                               ttGTTAATGGATTTTTTTTTTAAATTTTTTTTTTTTTAGACATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTc                                           <  1:409017/139‑1 (MQ=255)
                                                                                                 tttAATGTATTTTATTTTTACATTTTTTTGTTTTTAGAAATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                         <  2:610982/139‑1 (MQ=255)
                                                                                                 tttAATGGTTTTTTTTTTTAAATTTTTTTGTTTTTAGGCATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCat                                         <  2:73617/139‑1 (MQ=255)
                                                                                                 gTTAATGGATTTTATTTTTACATTTTTTTGTGTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGACATCGCTTTTAAAGTCGCGTTTTTCat                                         <  1:302469/140‑1 (MQ=255)
                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                        >  1:830381/1‑140 (MQ=255)
                                                                                                  ttAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCata                                        >  1:195048/1‑140 (MQ=255)
                                                                                                                        ttttttGTTTTTAGGCATAAAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGa                  <  1:835176/140‑1 (MQ=255)
                                                                                                                        ttttttGTGTTTAGGAATAAAAAAATCACCCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGa                  <  1:559900/140‑1 (MQ=255)
                                                                                                                                        atatAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATAc  <  1:119709/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                        
GGGTTATAAATGCGACTACCATGAAGTTTTTAATTGAAAGTATTGGGTTGCTGATAATTTGAGCTGTTCTATTCTTTTTAAATATCTATATAGGTCTGTTAATGGATTTTATTTTTACAATTTTTTGTGTTTAGGCATATAAAAATCAACCCGCCATATGAACGGCGGGTTAAAATATTTACAACTTAGCAATCAACCATTAACGCTTGATATCGCTTTTAAAGTCGCGTTTTTCATATCCTGTATACAGCTGACGCGGACGGGCAATCTTCATAC  >  NC_000913/752986‑753261

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 6 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑