Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913752,0650GC21.2% 73.4 / 12.9 34D35E (GAC→GAG) ybgQputative outer membrane protein
Reads supporting (aligned to +/- strand):  ref base G (13/13);  new base C (3/4);  total (17/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

ATGCATGGAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGC  >  NC_000913/751922‑752206
                                                                                                                                               |                                                                                                                                             
aTGCATGGAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGTTcgc                                                                                                                                           >  1:98344/1‑148 (MQ=255)
aTGCATGGAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGCTcgc                                                                                                                                           <  2:98344/148‑1 (MQ=255)
        aaTGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGCGACCGCATTGAta                                                                                                                                  <  4:76490/149‑2 (MQ=255)
            gTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATc                                                                                                                              >  5:104906/1‑149 (MQ=255)
             ttttGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCg                                                                                                                             <  6:104906/149‑1 (MQ=255)
             ttttGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGCGACCGCATTGAtatttca                                                                                                                             <  3:6663/149‑7 (MQ=255)
                    aCCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTATCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACa                                                                                                                       <  7:94878/148‑1 (MQ=255)
                                                      cTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTc                                                                                                                    <  6:175730/117‑1 (MQ=255)
                                                      cTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTc                                                                                                                    >  5:175730/1‑117 (MQ=255)
                                                      cTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCAt                                                                                    >  6:131456/1‑149 (MQ=255)
                                                          ttttGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGCGACCGCATTGAtatt                                                                                                                                <  3:34041/101‑4 (MQ=255)
                                                          ttttGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGCGACCGCATTGAtatt                                                                                                                                >  4:34041/1‑98 (MQ=255)
                                                                        aCCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGa                                                                  >  4:32723/1‑149 (MQ=255)
                                                                         ccGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAAc                                                                                                                        <  6:131704/94‑1 (MQ=255)
                                                                         ccGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAAc                                                                                                                        >  5:131704/1‑94 (MQ=255)
                                                                               cGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACaaaa                                                           <  3:32723/149‑1 (MQ=255)
                                                                                    aCAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCa                                                                                                 <  3:53335/106‑1 (MQ=255)
                                                                                    aCAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCa                                                                                                 >  4:53335/1‑106 (MQ=255)
                                                                                              cACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAat                                             <  5:131456/148‑1 (MQ=255)
                                                                                               aCCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGcc                                                                                >  8:142934/1‑112 (MQ=255)
                                                                                               aCCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGcc                                                                                <  7:142934/112‑1 (MQ=255)
                                                                                                  gggAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGa                                                                  <  4:118766/123‑1 (MQ=255)
                                                                                                  gggAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGa                                                                  >  3:118766/1‑123 (MQ=255)
                                                                                                        aaTGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAAtt                                                                                                          >  7:156080/1‑77 (MQ=255)
                                                                                                        aaTGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAAtt                                                                                                          <  8:156080/77‑1 (MQ=255)
                                                                                                          tGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTc                                 >  5:5351/1‑148 (MQ=255)
                                                                                                             cTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCAtt                                                                                               <  2:82607/83‑1 (MQ=255)
                                                                                                             cTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCAtt                                                                                               >  1:82607/1‑83 (MQ=255)
                                                                                                                                tGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATAGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACaaaa                                                           <  2:128406/100‑1 (MQ=255)
                                                                                                                                tGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATAGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACaaaa                                                           >  1:128406/1‑100 (MQ=255)
                                                                                                                                cgataaatCAATGCGCGACCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATg           >  1:64229/7‑148 (MQ=255)
                                                                                                                                      atCAATGCGCGACCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAggcg    >  2:75873/1‑149 (MQ=255)
                                                                                                                                         aaTGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAggcggc  >  5:59442/1‑148 (MQ=255)
                                                                                                                                         aaTGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAggcggc  <  6:59442/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                             
ATGCATGGAATGGTTTTGTCACCCTTTTTTTGCCAGTTAATTTTTTGCCCATTACTGATTTTGTTGTTATTCACCGCAACGCTAACAAAATATTCACCGGGAGCAATGACTCCTTTTTCCTTTAATAATGAAATATCAATGCGGTCGCGCATTGATTTATCGAGAACATTCAGATTGAATTCGACCGCCATTGCGCAAGGCATCGCCATGACCAGGCAGGATACAAAAGAGAGTCGATAAATATTCACGGTGTCCATACCTGATAAATATTTTATGAAAGGCGGC  >  NC_000913/751922‑752206

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