Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,002,2171.G25.0% 44.3 / 23.2 24coding (411/1071 nt)elfGputative fimbrial‑like adhesin protein
Reads supporting (aligned to +/- strand):  ref base . (10/8);  new base G (4/2);  total (14/10)
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.

GGGCATTCGCCGCGGATTCACCGTTAGTCTATTTAGGCGACAACTGGTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTCCCTGAAAATATTGCAGCCGTTCGTTGGCTCGGTCGTCATTCCTCCTA  >  NC_000913/1002081‑1002350
                                                                                                                                         |                                                                                                                                     
gggCATTCGCCGCGGATTCACCGTTAGTCTATTTAGGCGACAACTGGTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAgg                                                                                                                                     <  1:539367/139‑1 (MQ=255)
             ggATTCACCGTTAGTCTATTTAGGCGACAACTGGTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGATACCCGa                                                                                                                        <  2:586430/139‑1 (MQ=255)
                       ttAGTCTATTTAGGCGACAACTGGTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGATACCCGATGGCATAttt                                                                                                              >  2:355658/1‑139 (MQ=255)
                        tAGTCTATTTAGGCGACAACTGGTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGATACCCGATGGCATATTTg                                                                                                             <  1:355658/139‑1 (MQ=255)
                                    gCGACAACTGGTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCg                                                                                                  >  2:724427/1‑139 (MQ=255)
                                    gCGACAACTGGTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCgg                                                                                                 >  2:86067/1‑139 (MQ=255)
                                     cGACAACTGGTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGATACCCGATg                                                                                                                      <  1:312650/117‑1 (MQ=255)
                                     cGACAACTGGTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGATACCCGATg                                                                                                                      >  2:312650/1‑117 (MQ=255)
                                     cGACAACTGGTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGt                                                                                                <  1:86067/139‑1 (MQ=255)
                                              gTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGttt                                                                                        >  2:395858/1‑139 (MQ=255)
                                                                                 tATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCtt                                                    >  2:48281/1‑139 (MQ=255)
                                                                                   ttGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCttt                                                   >  2:370589/1‑139 (MQ=255)
                                                                                        ggTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTCCCTg                                             >  2:712713/1‑139 (MQ=255)
                                                                                                   gTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGt                                                                                                <  2:751463/77‑1 (MQ=255)
                                                                                                   gTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGt                                                                                                >  1:751463/1‑77 (MQ=255)
                                                                                                              gCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTCCCTGAAAATATTGCAGCCGTTCGTTg                       >  2:102069/1‑139 (MQ=255)
                                                                                                              gCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTCCCTGAAAATATTGCAGCCGTTCGTTg                       >  1:159165/1‑139 (MQ=255)
                                                                                                               cGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCa                                                                        <  2:143113/90‑1 (MQ=255)
                                                                                                               cGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCa                                                                        >  1:143113/1‑90 (MQ=255)
                                                                                                                 gTTCCTTTCGAAAACGTGGGCACAGGGGGGGATAcc                                                                                                                           >  1:719141/1‑36 (MQ=255)
                                                                                                                 gTTCCTTTCGAAAACGTGGGCACAGGGGGGGATAcc                                                                                                                           <  2:719141/36‑1 (MQ=255)
                                                                                                                     cTTTCGAAAACGTGGGCACAGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTCCCTGAAAATATTGCAGCCGTTCGTTGGCTCGgt                <  1:754155/139‑1 (MQ=255)
                                                                                                                       ttCGAAAACGTGGGCACAGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTCCCTGAAAATATTGCAGCCGTTCGTTGGCTCGgtcg              <  1:48281/139‑1 (MQ=255)
                                                                                                                         cGAAAACGTGGGCACAGGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTCCCTGAAAATATTGCAGCCGTTCGTTGGCTCGgtcgt             <  1:724427/139‑1 (MQ=255)
                                                                                                                                   ggCACAGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTCCCTGAAAATATTGCAGCCGTTCGTTGGCTCGGTCGTCATTCCTCCTa  >  1:239891/1‑139 (MQ=255)
                                                                                                                                         |                                                                                                                                     
GGGCATTCGCCGCGGATTCACCGTTAGTCTATTTAGGCGACAACTGGTACAAAATTAATGACTATCTTGCCGCCAAAGTTTTATTGCAGGTTAAAGGCAGTTCTCCTACTGCGGTTCCTTTCGAAAACGTGGGCACAGGGGGGGATACCCGATGGCATATTTGCGACCCTGGCGGTCAACGTTTAGGTGGGCAGGGGGCAAGCGGTAATAGCGGTAGCTTTTCCCTGAAAATATTGCAGCCGTTCGTTGGCTCGGTCGTCATTCCTCCTA  >  NC_000913/1002081‑1002350

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 35 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: