Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,634,0190TC58.8% 30.7 / 37.4 34I125V (ATT→GTT) ‡nohANohA
Reads supporting (aligned to +/- strand):  ref base T (7/7);  new base C (9/11);  total (16/18)
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.44e-02
Rejected as consensus: Frequency below/above cutoff threshold.

AACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGG  >  NC_000913/1633948‑1634111
                                                                       |                                                                                            
aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAg                                                                         >  1:385815/4‑93 (MQ=0)
aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAg                                                                         >  1:4527264/4‑93 (MQ=0)
aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAg                                                                         >  2:149950/4‑93 (MQ=0)
aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAg                                                                         >  2:4121961/4‑93 (MQ=0)
aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAg                                                                         >  2:4893489/4‑93 (MQ=0)
aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAg                                                                         <  1:595422/90‑1 (MQ=0)
aataTGTCGGTTTTCCAGCTCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAg                                                                         >  2:4666847/4‑93 (MQ=0)
                agctCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAAt                                                         >  2:2148797/4‑93 (MQ=0)
                agctCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAAt                                                         >  1:4471134/4‑93 (MQ=0)
                agctCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAAt                                                         <  2:422155/90‑1 (MQ=0)
                agctCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAAt                                                         >  1:5819651/4‑93 (MQ=0)
                agctCCGGGAAGCGCCGCTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAAt                                                         <  1:6597767/90‑1 (MQ=0)
                        gaagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc                                                 <  1:2931428/89‑1 (MQ=0)
                        gaagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTc                                                 <  1:4801541/89‑1 (MQ=0)
                         aagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCCTCCGCGACAGCACGAACGTGCAGAATGCGGTTTcc                                                <  1:3788575/90‑1 (MQ=0)
                         aagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTcc                                                <  1:6012364/90‑1 (MQ=0)
                         aagcgccgcTGAACCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTcc                                                <  2:493818/90‑1 (MQ=0)
                         aaccgccgcTGCACCGACAGGGGGATCCTGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTcc                                                <  1:4713297/90‑1 (MQ=0)
                                  tgaaCCGACAGAGGAACTCCATCAAGAATGCTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCa                                       <  2:1929950/90‑1 (MQ=0)
                                              ggtaTCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTg                           <  1:490187/90‑1 (MQ=1)
                                                            aagaCTGACAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTc             <  1:6111542/90‑1 (MQ=0)
                                                             atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc            <  1:972635/90‑1 (MQ=0)
                                                             atgcTGGCTACTTCTCCGGCCATCCGCGACAGCACGAACGTGCAGAATGCGGTTTCCACCACCTCAGCGGAATCTTTTGCATTTTTTAGTTcc            <  2:5770112/90‑1 (MQ=0)
                                                                   gCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCg      >  2:2216975/1‑93 (MQ=31)
                                                                   gCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCg      >  2:3185245/1‑93 (MQ=31)
                                                                   gCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCg      >  2:3187896/1‑93 (MQ=31)
                                                                     aaTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTc    <  1:5188413/93‑1 (MQ=33)
                                                                     aaTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTACACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTc    <  2:827503/93‑1 (MQ=18)
                                                                     aaTGTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTc    <  1:905378/93‑1 (MQ=18)
                                                                      aTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCg   >  2:874296/1‑93 (MQ=33)
                                                                      aTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCg   >  2:60215/1‑93 (MQ=33)
                                                                      aTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCg   >  2:556691/1‑93 (MQ=33)
                                                                      aTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCg   <  2:4299945/93‑1 (MQ=33)
                                                                      aTTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCg   >  2:3240968/1‑93 (MQ=33)
                                                                      aTTCCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCg   <  2:4673241/93‑1 (MQ=18)
                                                                       tttCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCgg  >  2:5056572/1‑93 (MQ=33)
                                                                       tttCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCgg  >  2:1666031/1‑93 (MQ=33)
                                                                       |                                                                                            
AACATGTCGGTTTTCCAGTTCCGGAAAACGCCGCTGCACCGACAGGGGGATCCCGTCGAGAATACTGGCAATTTCACCTGCGATCCGCGACAGCACGAAAGTACAGAATGCGGTTTCCACCACTTCAGCGGAGTCTCTGGCATTTTTCAGCTCCTGTGCGTCGG  >  NC_000913/1633948‑1634111

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