Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 989,969 T→C S130G (AGT→GGT)  pncB ← nicotinate phosphoribosyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913989,9690TC97.3% 136.9 / NA 37S130G (AGT→GGT) pncBnicotinate phosphoribosyltransferase
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (15/21);  minor base A (0/1);  total (15/22)
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.00e+00

GATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAA  >  NC_000913/989846‑990095
                                                                                                                           |                                                                                                                              
gATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACCGATAACCGCCAGCaa                                                                                                                 <  1:243309/139‑1 (MQ=255)
     aaGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACCGATAACCGCCAGCAAAGGaa                                                                                                            >  2:127465/1‑139 (MQ=255)
     aaGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACAGATGACCGCCAGCAAAGGaa                                                                                                            <  2:1216028/139‑1 (MQ=255)
      aGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGCCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACCGATAACCGCCAGCAAAGGAAc                                                                                                           <  1:1422115/139‑1 (MQ=255)
                          gcggTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCAc                                                                                       >  1:1266447/1‑139 (MQ=255)
                          gcggTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCAc                                                                                       >  2:470882/1‑139 (MQ=255)
                                aaCGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACg                                                                                 >  1:1374474/1‑139 (MQ=255)
                                       agaAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCt                                                                          <  2:210721/139‑1 (MQ=255)
                                         aaGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTa                                                                        <  2:431889/139‑1 (MQ=255)
                                            tCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAAc                                                                     >  2:786906/1‑139 (MQ=255)
                                                 tAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACCGATAACCGCCAGCAAAGGAAGTTCCCCGAGGATGACTTCACGCCACGGGCCGCTTAAACGAata                                                                <  1:1008853/139‑1 (MQ=255)
                                                         tttCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCtt                                                        <  1:470882/139‑1 (MQ=255)
                                                           tCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGc                                                      <  2:1258961/139‑1 (MQ=255)
                                                                       aGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGcc                                                                                    >  1:1097068/1‑97 (MQ=255)
                                                                       aGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGcc                                                                                    <  2:1097068/97‑1 (MQ=255)
                                                                            ttGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACg                                     <  1:539241/139‑1 (MQ=255)
                                                                            ttGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACg                                     <  1:707285/139‑1 (MQ=255)
                                                                            ttGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACg                                     <  1:786906/139‑1 (MQ=255)
                                                                             tGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCt                                                                          >  2:154274/1‑101 (MQ=255)
                                                                             tGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCt                                                                          <  1:154274/101‑1 (MQ=255)
                                                                              gcgcAACGTCGGCCTGCGGTGAGCGATAGTGATGTACCATTTCACCGATAACCGCCAGCaa                                                                                                                 <  1:429518/61‑1 (MQ=255)
                                                                              gcgcAACGTCGGCCTGCGGTGAGCGATAGTGATGTACCATTTCACCGATAACCGCCAGCaa                                                                                                                 >  2:429518/1‑61 (MQ=255)
                                                                               cgcAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACCGATAACCGc                                                                                                                       >  1:797259/1‑54 (MQ=255)
                                                                               cgcAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACCGATAACCGc                                                                                                                       <  2:797259/54‑1 (MQ=255)
                                                                                              cGGTGAGCGATAGCGATGTACCATTTCACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTa                   >  1:422081/1‑139 (MQ=255)
                                                                                                   agcgatagcgatGTACCATTTCACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCg              >  2:1498417/1‑139 (MQ=255)
                                                                                                           cgatGTACCATTTCACCGATAACCGCCAGCAAAGGaa                                                                                                            <  2:259906/37‑1 (MQ=255)
                                                                                                           cgatGTACCATTTCACCGATAACCGCCAGCAAAGGaa                                                                                                            >  1:259906/1‑37 (MQ=37)
                                                                                                             atGTACCATTTCACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCAt                                                   <  1:1027193/92‑1 (MQ=255)
                                                                                                             atGTACCATTTCACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCAt                                                   >  2:1027193/1‑92 (MQ=255)
                                                                                                               gTACCATTTCACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAc                                <  2:259997/109‑1 (MQ=255)
                                                                                                               gTACCATTTCACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGAc                                >  1:259997/1‑109 (MQ=255)
                                                                                                               gTACCATTTCACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa  <  1:769145/139‑1 (MQ=255)
                                                                                                               gTACCATTTCACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa  <  1:1498417/139‑1 (MQ=255)
                                                                                                               gTACCATTTCACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa  <  1:127465/139‑1 (MQ=255)
                                                                                                                        cACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa  >  1:816372/1‑130 (MQ=255)
                                                                                                                        cACCGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTaa  <  2:816372/130‑1 (MQ=255)
                                                                                                                           |                                                                                                                              
GATGGAAGCGCGACATATCAAGACCGGCGGTTAACGCCGAGAAGTCGACTAATTTGCTTTCCAGCGTGTCGAGGGCTTGCGCAACGTCGGCCTGCGGTGAGCGATAGCGATGTACCATTTCACTGATAACCGCCAGCAAAGGAACTTCCCAGAGGATGACTTCACGCCACGGGCCGCTTAAACGAATATCCAGCTTGCCATTATCGTTGGACACGGTGACTTGTTCCGGGTTAAAGCGGAACTCGCGTAA  >  NC_000913/989846‑990095

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