Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913757,1790GT40.0% 8.2 / 13.7 15E425* (GAG→TAG) sdhAsuccinate dehydrogenase, flavoprotein subunit
Reads supporting (aligned to +/- strand):  ref base G (4/5);  new base T (4/2);  total (8/7)
Fisher's exact test for biased strand distribution p-value = 6.08e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: E-value score below prediction cutoff.
Rejected as consensus: Frequency below/above cutoff threshold.

GGGACTGTTTGCCGTTGGTGAAATCGCTTGTGTATCGGTACACGGCGCTAACCGTCTGGGCGGCAACTCGCTGCTGGACCTGGTGGTCTTTGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAA  >  NC_000913/757049‑757315
                                                                                                                                  |                                                                                                                                        
gggACTGTTTGCCGTTGGTGAAATCGCTTGTGTATCGGTACACGGCGCTAACCGTCTGGGCGGCAACTCGCTGCTGGACCTGGTGGTCTTTGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTgcgc                                                                                                                         >  8:66446/1‑148 (MQ=255)
                   gAAATCGCTTGTGTATCGGTACACGGCGCTAACCGTCTGGGCGGCAACTCGCTGCTGGACCTGGTGGTCTTTGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGt                                                                                                     >  6:97034/1‑149 (MQ=255)
                    aaaTCGCTTGTGTATCGGTACACGGCGCTAACCGTCTGGGCGGCAACTCGCTGCTGGACCTGGTGGTCTTTGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCTAGCAGGgcg                                                                                                                                 <  1:94937/120‑1 (MQ=255)
                    aaaTCGCTTGTGTATCGGTACACGGCGCTAACCGTCTGGGCGGCAACTCGCTGCTGGACCTGGTGGTCTTTGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCTAGCAGGgcg                                                                                                                                 >  2:94937/1‑120 (MQ=255)
                    aaaTCGCTTGTGTATCGGTACACGGCGCTAACCGTCTGGGCGGCAACTCGCTGCTGGACCTGGTGGTCTTTGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGtt                                                                                                    <  7:66446/149‑1 (MQ=255)
                                             cgcTAACCGTCTGGGCGGCAACTCGCTGCTGGACCTGGTGGTCTTTGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAAcc                                                                           <  5:97034/149‑1 (MQ=255)
                                                         gggCGGCAACTCGCTGCTGGACCTGGTGGTCTTTGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACaat                                                                >  2:149955/1‑148 (MQ=255)
                                                                               cTGGTGGTCTTTGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATg                                                                                                      >  5:146757/1‑88 (MQ=255)
                                                                               cTGGTGGTCTTTGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATg                                                                                                      <  6:146757/88‑1 (MQ=255)
                                                                                          tGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATcc                                 <  1:149955/146‑1 (MQ=255)
                                                                                               gcgcgGCAGGTCTGCATCTGCAAGAGTCTATCGCCTAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCACGTTGCTATCCGTaaaaa                          >  4:103125/1‑146 (MQ=255)
                                                                                                              atctgcaAGAGTCTATCGCCTAGCAGGGCGGACTGCGAGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGATCAATAAGCGTAACGGTGAAGATGCGGTGGCGATCCGTAAAGCGATGCAAGAATGTATg          >  3:75412/1‑149 (MQ=255)
                                                                                                                    aaGAGTCTATCGCCTAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACaa                                                                 >  1:126370/1‑88 (MQ=255)
                                                                                                                    aaGAGTCTATCGCCTAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACaa                                                                 <  2:126370/88‑1 (MQ=255)
                                                                                                                      gagTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATaa  <  5:108043/149‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                        
GGGACTGTTTGCCGTTGGTGAAATCGCTTGTGTATCGGTACACGGCGCTAACCGTCTGGGCGGCAACTCGCTGCTGGACCTGGTGGTCTTTGGTCGCGCGGCAGGTCTGCATCTGCAAGAGTCTATCGCCGAGCAGGGCGCACTGCGCGATGCCAGCGAGTCTGATGTTGAAGCGTCTCTGGATCGCCTGAACCGCTGGAACAATAATCGTAACGGTGAAGATCCGGTGGCGATCCGTAAAGCGCTGCAAGAATGTATGCAGCATAA  >  NC_000913/757049‑757315

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