Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,068,1730TG41.7% 27.3 / 10.4 24G31G (GGT→GGG) ymdFconserved protein YmdF
Reads supporting (aligned to +/- strand):  ref base T (3/11);  new base G (10/0);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 1.51e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CGTTAATTCTATGCAATATGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTC  >  NC_000913/1068034‑1068297
                                                                                                                                           |                                                                                                                            
gacagaTTCTATGCAATATGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGt                                                                                                                              <  1:80505/135‑1 (MQ=255)
cgTTAATTCTATGCAATATGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAATACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGt                                                                                                                              <  2:362675/140‑1 (MQ=255)
cgTTAATTCTATGCAATATGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGt                                                                                                                              <  1:245303/140‑1 (MQ=255)
cgTTAATTCTATGCAATATGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGt                                                                                                                              <  1:314693/140‑1 (MQ=255)
 gTTAATTCTATGCAATATGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTa                                                                                                                             <  1:427454/140‑1 (MQ=255)
           tGCAATATGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTTCGGCAATTTTGCAGAAAACCGCGAAAAAGCATCAGAAGCAGGGAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGaa                                                                                                                   >  2:415121/1‑140 (MQ=255)
                tatGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGAcc                                                                                                              >  1:397470/1‑140 (MQ=255)
                tatGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGAcc                                                                                                              >  1:212305/1‑140 (MQ=255)
                tatGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCCGCAATTTTGGCGAAGACCGCGAAAGAGCATCAGAAGCAGGGAAAAAAGGGGGACAGCACACCGGGGGGAATTTCAAGaactgtt                                                                                                              >  2:80505/1‑135 (MQ=255)
                   gatgTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAAAGCATCAGAAGCAGGGAAAAAAGGGGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGc                                                                                                           >  2:138540/1‑140 (MQ=255)
                          ataCCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGGAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCCCAgcgcgc                                                                                                    >  2:53633/1‑140 (MQ=255)
                                 gACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGaa                                                                                             >  2:367673/1‑140 (MQ=255)
                                 gACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGGAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGaa                                                                                             >  1:130531/1‑140 (MQ=255)
                                 gACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGGAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGaa                                                                                             >  2:168093/1‑140 (MQ=255)
                                  aCGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGAAg                                                                                            >  1:18285/1‑140 (MQ=255)
                                           aggaaTGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGGAGTCCGCGAAAGAGCATCAGGAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTaaaa                                                                                   <  1:27374/137‑1 (MQ=255)
                                           aGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTaaaa                                                                                   <  2:324607/140‑1 (MQ=255)
                                                  gCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTTAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAggggg                                                                            >  1:45805/1‑140 (MQ=255)
                                                       ccATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAaga                                                                       >  1:433318/1‑140 (MQ=255)
                                                                  gTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGc                                                            <  1:260180/140‑1 (MQ=255)
                                                                  gTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGGAAAAAAGGGGGGAAGAGCAGTCACGGc                                                            >  1:376104/1‑140 (MQ=255)
                                                                       ggCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGT‑AAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGc                                                      <  2:223943/140‑1 (MQ=255)
                                                                                 cAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAg                                             <  2:397470/140‑1 (MQ=255)
                                                                                 cAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAg                                             <  1:438457/140‑1 (MQ=255)
                                                                                         cgcgAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTa                                     >  1:134522/1‑140 (MQ=255)
                                                                                                  gCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGAc                            <  2:433318/140‑1 (MQ=255)
                                                                                                      cAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAAt                        <  1:149855/140‑1 (MQ=255)
                                                                                                                         tGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTc     <  2:124974/140‑1 (MQ=255)
                                                                                                                            acagcacagcGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTc  <  1:310999/140‑1 (MQ=255)
                                                                                                                                           |                                                                                                                            
CGTTAATTCTATGCAATATGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTC  >  NC_000913/1068034‑1068297

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