Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,852,271 C→T 100% R129C (CGG→TGC)  pncA → nicotinamidase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,852,2710CT91.6% 58.4 / ‑1.2 24R129C (CGG→TGC) pncAnicotinamidase
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base T (12/10);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.82e-01
Rejected as polymorphism: E-value score below prediction cutoff.
Rejected as polymorphism: Variant not supported by required number of reads on each strand.

TCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGAT  >  NC_000913/1852141‑1852407
                                                                                                                                  |                                                                                                                                        
tCACCGTGTCAAGAACGGTGCAGGCGCACAATTACATCCTCTACTGCACCAAAAAGCGATCGCAGCGGTGCTCGATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACct                                                                                                                              <  1:65424/143‑1 (MQ=255)
 cACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACctc                                                                                                                             >  1:23580/1‑143 (MQ=255)
    tgtgtgCAGAACAGTGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGAGCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGTCCTTTGGGAGACAACctcgct                                                                                                                          >  2:263462/1‑143 (MQ=255)
     gtgtgCAGAACAGTGAAGGCGCACAATTACATCCGTTAGTGCACCAAAAAGCGATCGCAGCGGTGTTTCATAAAGTCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCtctctc                                                                                                                         <  2:537428/143‑1 (MQ=255)
                  cgAAGTCGCACAATTACTTCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTAc                                                                                                            <  2:479734/142‑1 (MQ=255)
                   gtaGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTAcg                                                                                                           <  2:41805/141‑1 (MQ=255)
                               ttACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAAt                                                                                               <  1:566051/143‑1 (MQ=255)
                                       gTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCGTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGCCTGGTTACGCGATCATGAAATCGATGAAt                                                                                       >  1:283510/1‑143 (MQ=255)
                                            tGCCCCGAAAAGCGCTCGCAGCGGTGTGCCATAACGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGTTAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATc                                                                                  <  2:435263/143‑1 (MQ=255)
                                            tGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCATTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGCATTGATc                                                                                  >  1:171712/1‑143 (MQ=255)
                                                 cAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTAt                                                                             >  1:332436/1‑143 (MQ=255)
                                                   aaaaGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATgg                                                                           <  1:400883/143‑1 (MQ=255)
                                                           tCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGTTGACTGGTTACGCGATCATGAAATCGCTGCATTGCGCGTTATGGGACTGGCt                                                                   >  2:310790/1‑143 (MQ=255)
                                                                                gCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTCCGTGAACttt                                              >  1:530612/1‑143 (MQ=255)
                                                                                     caTCATTTAGTTGACAGTCACAGATCCTTCTTTGCTAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCGACTGACTATTGCGTGAAGTTTACCGt                                         <  1:222612/142‑1 (MQ=255)
                                                                                     agtCCTTTAGTTGACAGTGACAGTGCCTTTTTTGATACCGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGt                                         <  1:182972/141‑1 (MQ=255)
                                                                                                     gTTACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAg                         >  1:261289/1‑143 (MQ=255)
                                                                                                      ttACAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGCTCGTTATGGGCCTGGCTACTGACTATTGCGTGACGTTTACCGTGCTGGACGCGTTACAGt                        >  1:90195/1‑143 (MQ=255)
                                                                                                         cAGTGCCTTTTTTGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGGGCTGGACGCGGTACAGTGAg                     >  2:318315/1‑143 (MQ=255)
                                                                                                                ttttttGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATaa              >  1:274639/1‑143 (MQ=255)
                                                                                                                   tttGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGGTGCTGACGAGACCTTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGt           >  1:266254/1‑143 (MQ=255)
                                                                                                                   tttGATAACGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTACTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGt           >  1:503878/1‑143 (MQ=255)
                                                                                                                           cGGCCGTTGCCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGa   <  2:332436/143‑1 (MQ=255)
                                                                                                                            ggCCGTTGCCAGCCAACCTCTCTCGATGACTGGGTACGCGATCACGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGAt  <  1:79873/143‑1 (MQ=255)
                                                                                                                                  |                                                                                                                                        
TCACTGTGTGCAGAACAGTGAAGGCGCACAATTACATCCGTTACTGCACCAAAAAGCGATCGCAGCGGTGTTCCATAAAGGCGAAAATCCTTTAGTTGACAGTTACAGTGCCTTTTTTGATAACGGCCGTCGGCAGAAAACCTCTCTCGATGACTGGTTACGCGATCATGAAATCGATGAATTGATCGTTATGGGCCTGGCTACTGACTATTGCGTGAAGTTTACCGTGCTGGACGCGTTACAGTTAGGTTATAAGGTAAACGTGAT  >  NC_000913/1852141‑1852407

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