Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,261,0550AC57.1% 7.0 / 15.7 14V455G (GTG→GGG) tdcEpyruvate formate‑lyase 4/2‑ketobutyrate formate‑lyase
Reads supporting (aligned to +/- strand):  ref base A (6/0);  new base C (0/8);  total (6/8)
Fisher's exact test for biased strand distribution p-value = 3.33e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.25e-03

AGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCT  >  NC_000913/3260923‑3261173
                                                                                                                                    |                                                                                                                      
aGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCg                                                                                                       >  1:234544/1‑150 (MQ=255)
                 caccccccgccctcccATGAAGGGATCGGGGCTGTCCATCCCTTTTTCGTAATCCACCCCGCCTTCCATCGGCGGTGCTTTTTTCGGCCCGTCCAGAATTTTCCGCTCCTCGTCCCCCCCGCCGTTAATTGCGTAGagcag                                                                                               <  1:84527/127‑1 (MQ=255)
                   ccGCCACCCAGCCCATGAAGTGATCGAGTCTTTCCATCATTTTGTCGTAATCCAGCCCTTCGCCCACCAGCGGTGTTTTTTTCGGCCCGACCTGAATCTTCGGTTTCTCTTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCg                                                                                       <  1:254115/147‑1 (MQ=255)
                        agccagTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGtt                                                                                >  1:254202/1‑149 (MQ=255)
                             gtacgtaaaGTGACGGAGGCGTTCCGTCTCTTTGCCGTAATCCCGCCCGAGGTCCATCAGGGGTGGGGTTTTCTCCCGGCCTTGAACTTTATGTCCCTCTCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGcgcg                                                                          <  1:131036/143‑1 (MQ=255)
                              tCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTg                                                                        >  2:352671/1‑151 (MQ=255)
                                                    ccATCCCTTTGTCTTAACCCACCCCGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGCCCTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCAttt                                                      <  2:146546/147‑1 (MQ=255)
                                                           tttGTCGAACCCCAGCCCGTCGCCCCTCGGCGGTGCTGTTTTCGGCCCGCCGTGAATCTTCAGCTTCTCGCCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTAcc                                               <  1:248727/147‑1 (MQ=255)
                                                                   aacccAGCCCGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTg                                    <  2:337542/147‑1 (MQ=255)
                                                                    aTCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTggg                                  >  1:306527/1‑151 (MQ=255)
                                                                                     tCAGCGGGGCTGTTTTCGCCCCGACCGGAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCaa                  <  1:177973/150‑1 (MQ=255)
                                                                                         cGGTGCTGTTTTCGCCCCGCCCTAAATCTTCAGCTTCTCGTCCCCCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATcg               <  2:321087/149‑1 (MQ=255)
                                                                                                   ttCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGc   >  2:220591/1‑151 (MQ=255)
                                                                                                    tCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGCAGAGCAGCGTTTTCGCCAGGTTAGCCCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCt  >  1:38886/1‑151 (MQ=255)
                                                                                                                                    |                                                                                                                      
AGCGCGCTGATGTACTGCACCGCCAGCCAGTCCATGAAGTGATCGAGGCTGTCCATCACTTTGTCGTAATCCAGCACGTCGTCCATCAGCGGTGCTGTTTTCGGCCCGACCTGAATCTTCAGCTTCTCGTCCACCCCGCCGTTAATTGCGTAGAGCAGCGTTTTCGCCAGGTTAGCGCGTGCACCAAAGAACTGCATTTGCTTACCAATCACCATTGGGCTGACGCAGCAGGCAATCGCGTAATCGTCGCT  >  NC_000913/3260923‑3261173

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