Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,233,609 C→T Q4* (CAG→TAG)  yeiS → DUF2542 domain‑containing protein YeiS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,233,6090CT94.1% 44.1 / NA 17Q4* (CAG→TAG) yeiSDUF2542 domain‑containing protein YeiS
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (8/8);  minor base G (0/1);  total (8/9)
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.93e-01

ATTCGTGAATTTGCCGCCCGTTTCGGTGCGCTGGCTGACCTTTATATTTTTAAACGTGAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCGGTGTATGTCTGGCGAGCAAAAAATCCCGGACTCTTTTTCGCTTATATGGTGGCATATATCGGCTTCGGAATTTTATCTATCGGCATGATTGTTTAT  >  NC_000913/2233415‑2233821
                                                                                                                                                                                                  |                                                                                                                                                                                                                    
aTTCGTGAATTTGCCGCCCGTTTCGGTGCGCTGGCTGACCTTTATATTTTTAAACGTGAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACGTGAATTGCAAAAGAAACAAGGAAAGTATTGATGGATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGaa                                                                                                                                                              <  1:325519/251‑1 (MQ=255)
             ccgccCGTTTCGGTGCGCTGGCTGACCTTTATATTTTTAAACGTGAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTACGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGGATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCa                                                                                                                                                           >  1:292144/1‑241 (MQ=255)
                   gTTTCGGTGCGCTGGCTGACCTTTATATTTTGAAACGGGAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATTCACCAGGTACCGGAAGATGCGCAGGGTTATCCCGCCGGCGCACCCGAACAGTTACTTGAATTACAAAAGAAACAATGAAAGTAGTTATGGATGTAGAGCAGTTTTGTGTCGGTGCCGTTTTTTTCCTTAGTCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGgcgcg                                                                                                                                           <  1:348131/251‑1 (MQ=255)
                                      ccTTTATATTTTTAAACGTGAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTAGGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATTTATAGCAGTTTTTTGTCGTTGCCGTTGTTTTCCTTATCCCGATATTGTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACg                                                                                                                        <  1:255352/251‑1 (MQ=255)
                                         ttATATTTTTAAACGTGACCCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGAGGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATTTATAGCAGTTTTTTGTCGTTGCCTTTTTCTTCCTTATCCCGATATTTTTTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACggg                                                                                                                      <  1:252546/250‑1 (MQ=255)
                                              tttttAAACGTGAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGgcgcn                                                                                                                                           <  1:215998/224‑2 (MQ=255)
                                                                         ggCCGCTGGACCCTATTCCGGCTATGCACCCGGTACCGGAAGATGCGCAGGGGTCTCCAGCCGTCACACCCGAACAGTCACTTGAATTACAAAAGAAACAAGGAATGTAGATATGGATGAATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGATGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGACCCGGTGTATGTCTGGc                                                                                     >  1:145839/1‑251 (MQ=255)
                                                                                        ttCCGGCTATGCCCCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACACGAACAGTTACTTGACTTACAAAAGAAACAAGGAAAGTAGTTATTGATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCGGTGTATGTCTGGCGAGCAAAAAATccc                                                                       <  1:20763/250‑1 (MQ=255)
                                                                                                      cAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGCTATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCGGTGTATGTCTGGCGAGCAAAAAATCCCGGACTCTTTTTCGCt                                                        >  1:58956/1‑251 (MQ=255)
                                                                                                      cAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCGGTGTATGTCTGGCGAGCAAAAAATCCCGGACTCTTTTTCGCt                                                        >  1:58967/1‑251 (MQ=255)
                                                                                                        ggTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTcc                                                                                                                                                                                        >  1:98785/1‑121 (MQ=255)
                                                                                                               gAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCGGTGTATGTCTGGCGAGCAAAAAATCCCGGACTCTTTTTCGCTTATAtggtg                                               <  1:236894/251‑1 (MQ=255)
                                                                                                                                      cnTCNCACCCNAACANTTACTTGAATTACAANAGAAACAAGGAAAGTAGTTNNGNATGTAnagcag                                                                                                                                                                                                                 <  2:134466/66‑1 (MQ=255)
                                                                                                                                               caaaCAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGTATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTg                                                                                                                                                                        <  1:243477/96‑1 (MQ=255)
                                                                                                                                               cGAACAGTNANNNGNNTTACAAAAGAAACAAGGAAAGTAGTTATGGATGTATAGCAGTTTNTTGTCGTTGCCGTTTNTTNCCTTNTCCNGATATTTTg                                                                                                                                                                        >  2:243477/1‑98 (MQ=255)
                                                                                                                                                     gTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGCACCGGTGTATGTCTGGCGAGCAAAAAATcc                                                                        >  1:148376/1‑188 (MQ=255)
                                                                                                                                                            gAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATTTTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCGGTGTATGTCTGGCGAGCAAAAAATCCCGGACTCTTTTTCGCTTATATGGTGGCATATATCGGCTTCGGAATTTTATCTATCGGCATGATTGTTTAt  >  1:194582/1‑251 (MQ=255)
                                                                                                                                                                                    tagtTATGGATGTATAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCGGTGTATGTCTGGCg                                                                                    <  1:194635/145‑1 (MQ=255)
                                                                                                                                                                                                  |                                                                                                                                                                                                                    
ATTCGTGAATTTGCCGCCCGTTTCGGTGCGCTGGCTGACCTTTATATTTTTAAACGTGAACCGCGTTTTTTAGGGCCGCTGGTCCCTATTCCGGCTATGCACCAGGTACCGGAAGATGCGCAGGGGTATCCCGCCGTCACACCCGAACAGTTACTTGAATTACAAAAGAAACAAGGAAAGTAGTTATGGATGTACAGCAGTTTTTTGTCGTTGCCGTTTTTTTCCTTATCCCGATATTTTGTTTCCGCGAAGCATGGAAAGGCTGGCGCGCAGGCGCGATTGATAAACGGGTTAAAAATGCACCGGAACCGGTGTATGTCTGGCGAGCAAAAAATCCCGGACTCTTTTTCGCTTATATGGTGGCATATATCGGCTTCGGAATTTTATCTATCGGCATGATTGTTTAT  >  NC_000913/2233415‑2233821

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