Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,068,1730TG31.2% 47.1 / 11.5 32G31G (GGT→GGG) ymdFconserved protein YmdF
Reads supporting (aligned to +/- strand):  ref base T (11/11);  new base G (10/0);  total (21/11)
Fisher's exact test for biased strand distribution p-value = 6.02e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.82e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GCAATATGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCA  >  NC_000913/1068046‑1068306
                                                                                                                               |                                                                                                                                     
gCAATATGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAAt                                                                                                                           >  1:371153/1‑140 (MQ=255)
      tgatgTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCg                                                                                                                     <  1:514108/140‑1 (MQ=255)
         tgTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGGAAAAAAAGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAg                                                                                                                  >  2:514108/1‑140 (MQ=255)
                aCCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGAc                                                                                                                        <  1:352203/127‑1 (MQ=255)
                aCCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGAc                                                                                                                        >  2:352203/1‑127 (MQ=255)
                     gagacaGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGaa                                                                                                      <  1:480708/135‑1 (MQ=255)
                     gACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGa                                                                                                       <  1:230190/139‑1 (MQ=255)
                     gACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGaa                                                                                                      >  2:231755/1‑140 (MQ=255)
                     gACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGa                                                                                                       >  2:230190/1‑139 (MQ=255)
                     gACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGGAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGaa                                                                                                      >  2:402995/1‑140 (MQ=255)
                     gACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGGAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAAAGACCCGCAGCGCGCATCTGaa                                                                                                      >  2:114239/1‑140 (MQ=255)
                        ggAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCa                                                                                                   <  2:491273/140‑1 (MQ=255)
                          aGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGGAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCCGCGCGCATCTGAAGCAgg                                                                                                 >  1:64421/1‑140 (MQ=255)
                                       cAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAggtggt                                                                                    <  1:209554/140‑1 (MQ=255)
                                       cAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAggtggt                                                                                    <  1:585536/140‑1 (MQ=255)
                                       cAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGGAAAAAAGGTGGACAGCACAGCGGGGGGAAATTTAAGAATGACCCGCAGCGCGCATCTGTAGCAGGTaaaaactgtctc                                                                                    >  2:480708/1‑133 (MQ=255)
                                         aaCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCCGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTaa                                                                                  >  1:125763/1‑140 (MQ=255)
                                         aaCCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGGAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGGAAAAAAGGGGGTaa                                                                                  >  2:58341/1‑140 (MQ=255)
                                           ccATCGAGGCGGTTCCGGCAATTTTTCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGGGGACAGCACAGCGGGGGGAATTTCAAGAAAGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGGGGTAaga                                                                                >  2:288672/1‑140 (MQ=255)
                                                      gTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGGAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGc                                                                     >  1:470768/1‑140 (MQ=255)
                                                                          gaCCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCggg                                                 <  2:371153/140‑1 (MQ=255)
                                                                          gaCCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCggg                                                 <  1:231755/140‑1 (MQ=255)
                                                                           aCCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGc                                                >  2:384545/1‑140 (MQ=255)
                                                                             cgcgAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTa                                              >  2:493411/1‑140 (MQ=255)
                                                                               cgAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAAt                                            >  2:350821/1‑140 (MQ=255)
                                                                               cgAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAAt                                            >  1:192149/1‑140 (MQ=255)
                                                                                  aaGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCaa                                         >  1:16226/1‑140 (MQ=255)
                                                                                           aGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATg                                <  1:402995/140‑1 (MQ=255)
                                                                                             aaGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCa                              <  1:281733/140‑1 (MQ=255)
                                                                                                  ggTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCAttttt                         >  2:210438/1‑140 (MQ=255)
                                                                                                  ggTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCAttttt                         >  2:392547/1‑140 (MQ=255)
                                                                                                                     acagcGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCaaaaa      >  2:322393/1‑140 (MQ=255)
                                                                                                                         cGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCa  >  1:545857/1‑140 (MQ=255)
                                                                                                                               |                                                                                                                                     
GCAATATGATGTCTATACCCAGACGGAGGTCAGTAATGGCAAACCATCGAGGCGGTTCCGGCAATTTTGCAGAAGACCGCGAAAGAGCATCAGAAGCAGGTAAAAAAGGTGGACAGCACAGCGGGGGTAATTTCAAGAATGACCCGCAGCGCGCATCTGAAGCAGGTAAAAAAGGTGGTAAGAGCAGTCACGGCAAAAGCGACAACTAGCCGGGCTAATCAATGACGAATGCATTTTTGTCTGTAGCTCGTCAAAAAGCCA  >  NC_000913/1068046‑1068306

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