Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,095,695 T→C C191R (TGC→CGC)  hisA → 1‑(5‑phosphoribosyl)‑5‑[(5‑ phosphoribosylamino)methylideneamino]imidazole‑4‑ carboxamide isomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,095,6950TC96.9% 101.3 / ‑4.1 33C191R (TGC→CGC) hisA1‑(5‑phosphoribosyl)‑5‑[(5‑ phosphoribosylamino)methylideneamino]imidazole‑4‑ carboxamide isomerase
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (17/14);  total (18/15)
Fisher's exact test for biased strand distribution p-value = 4.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.35e-01

GGGCGTTTCACTGGAACAACTGGTGGAAACCTATCTGCCCGTCGGCCTGAAACATGTGCTGTGTACCGATATCTCGCGCGACGGCACGCTGGCAGGCTCTAACGTCTCTTTATATGAAGAAGTGTGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCGCGGCGTAATAGTTGGTCGGGCATTACTGGAAGGTAAATTCACC  >  NC_000913/2095571‑2095826
                                                                                                                            |                                                                                                                                   
gggCGTTTCACTGGAACAACTGGTGGAAACCTATCTGCCCGTCGGCCTGAAACATGTGCTGTGTACCGATATCTCGCGCGACGGCACGCTGGCAGGCTCTTACGTCTCTTTATATGAGGAAGGGCGCGCCAGATATCCGc                                                                                                                      >  2:388496/1‑140 (MQ=255)
   cGTTTCACTGGAACAACTGGTGGAAACCTATCTGCCCGTCGGCCTGAAACATGTGCTGTGTACCGATATCTCGCACGACGGCACGCTGGCAGGCTCTAACGTATCTTTATGTGAAGAAGTGCGCCCCAGATATCCGCAgg                                                                                                                   >  1:418173/1‑140 (MQ=255)
               acaacTGGTGGAAACCTATCTGCCCGTCGGCCTGAAACATGTGCTGTGTACCGATATCTCGCGCGACGGCACGCTGGCAGGCTCTAACGTCTCTTTATATGAAGAAGTGCGCGCCAGATATCCGCAGGTGGCATTTCAGt                                                                                                       >  1:75130/1‑140 (MQ=255)
                caacTGGTGGAAACCTATCTGCCCGTCGGCCTGAAACATGTGCTGTGTACCGATATCTCGCGCGACGGCACGCTGGCAGGCTCTAACGTCTCTCTATATGAAGAAGTGGGCGCCAGAAATCCCCAGGTGGCATTTCAGtc                                                                                                      >  2:255585/1‑140 (MQ=255)
                        ggAAACCTATCTGCCCGTCGGCCTGAAACATGTGCTGTGTACCGATATCTCGCGCGACGGCACGCTGGCAGGCTCTAACGTCTCTTTATATGAAGAAGTGCGTGCCAGATTTCCTCAGGTGGCATTTCAGTCCTCcggcg                                                                                              >  1:372241/1‑140 (MQ=255)
                            aCCTATCTGCCCGTCGGCCTGAAACATGTGCTGTGTACCGATATCTCGCGCGACGGCACGCTGGCAGGCTCTAACGTCTCTTTATATGAAGAAGTGCGCGCCGGATATCCCCAGGTGGCATTTCAGTCCTCCGGCGGTAt                                                                                          >  1:316622/1‑140 (MQ=255)
                               tATCTGCCCGTCGGCTTGAAACGTGCGCTCTGTACCGCTATCTCGCGCGACGGCACGCTGGCAGGCTCTAACGTCTCTCTATATGAAGAAGTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTgg                                                                                       <  1:255585/140‑1 (MQ=255)
                                     cccGTCGGCCTGAAACATGTGCTGTGTACCGATATCTCGCGCGACGGCACGCTGGCAGGCTCTAACGTCTCTTTATATGAAGAAGTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACAt                                                                                 <  2:812770/140‑1 (MQ=255)
                                                          cTGTGTACCGATATCTCGCGCGACGGCACGCTGGCAGGCTCTAACGTCTCTTTATATGAAGAAGTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCg                                                            <  1:775928/140‑1 (MQ=255)
                                                          cTGTGTACCGATATCTCCCGCGCCGGCACGCAGGCAGGCTCTAACGTCTCTTTATATGAAGAAGTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCg                                                            <  2:316622/140‑1 (MQ=255)
                                                               tACCGATATCTCGCGCGACGGCACGCTGGCAGGCTCTAACGTCTCTTTATATGAAGAAGTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCa                                                       <  2:44790/140‑1 (MQ=255)
                                                               tACCGATATCTCGCGCGACGGCACGCTGGCAGGCTCTAACGTCTCTTTATATGAAGAAGTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCa                                                       >  1:674666/1‑140 (MQ=255)
                                                               tACCGATATCTCGCGCGACGGCACGCTGGCAGGCTCTAACGTCTCTTTATATGAAGAAGTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCa                                                       >  1:539123/1‑140 (MQ=255)
                                                                       tctcGCGCGACGGCACGCTGGCAGGCTCTAACGTCTCTTTATATGAAGAAGTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGTGTGGCACTGGtgtg                                               >  1:211130/1‑140 (MQ=255)
                                                                       tctcGCGCGACGGCACGCTGGCAGGCTCTAACGTCTCTTTATATGAAGAAGTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGtgtg                                               >  2:83028/1‑140 (MQ=255)
                                                                         tcGCGCGACGGCACGCTGGCAGGCTTTAAAGTCTCTTGATATGAAGAAGTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTgcg                                             <  1:388496/140‑1 (MQ=255)
                                                                            cgcgACGGCACGCTGGCAGGCTCTAACGTCTCTTTATATGAAGAAGTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcgg                                          >  1:752959/1‑140 (MQ=255)
                                                                              cgACGGCACGCTGGCAGGCTCTAACGTCTCTTTATATGAAGAAGTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCgcggcg                                        >  1:131659/1‑140 (MQ=255)
                                                                               gACGGCACGCTGGCAGGCTCTAACGTCTCTTTATATGAAGAAGTGTGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCGCGGCGt                                       <  2:751928/140‑1 (MQ=255)
                                                                                 cGGCACGCTGGCAGGCTCTAACGTCTCTTTATATGAAGAAGTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCGCGGCGTaa                                     >  1:142764/1‑140 (MQ=255)
                                                                                        cTGGCAGGCTCTAACGTCTCTTTATATGAAGGAGTCCGCCCCAGATATGCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCGCGGCGTAATAGTTgg                              <  2:303964/140‑1 (MQ=255)
                                                                                        cTGGCAGGCTCTAACGTCTCTTTATATGAAGAAGTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGTGTGGCACTGGTGTGCGCGGCGTAATAGTTgg                              <  2:211130/140‑1 (MQ=255)
                                                                                           gCAGGCTCTAACGTCTCTTTATATGAAGAAGTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCGCGGCGTAATAGTTGGTCg                           <  1:126329/140‑1 (MQ=255)
                                                                                             aGGCTCTAACGTCTCTTTATATGAAGAAGTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCGCGGCGTAATAGTTGGTCggg                         <  2:142764/140‑1 (MQ=255)
                                                                                               gCTCTAACGTCTCTTTATATGAAGAAGTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCGCGGCGTAATAGTTGGTCGGGCa                       <  1:715900/140‑1 (MQ=255)
                                                                                                     aCGTCTCTTTATATGAAGAAGTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCGCGGCGTAATAGTTGGTCGGGCATTACTg                 >  2:211918/1‑140 (MQ=255)
                                                                                                      cGTCTCTTTATATGAAGAAGTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCGCGGCGTAATAGTTGGTCGGGCATTACTgg                <  2:131659/140‑1 (MQ=255)
                                                                                                         ctctTTATATGAAGAAGTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCGCGGCGTAATAGTTGGTCGGGCATTACTGGAAg             >  2:87845/1‑140 (MQ=255)
                                                                                                          tctTTATATGAAGAAGTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCGCGGCGTAATAGTTGGTCGGGCATTACTGGAAgg            >  2:434221/1‑140 (MQ=255)
                                                                                                           ctTTATATGAAGAAGTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCGCGGCGTAATAGTTGGTCGGGCATTACTGGAAGGt           >  1:692198/1‑140 (MQ=255)
                                                                                                                 atGAAGAAGTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGGGGCGGCCCTGCGTGGCACTGGTGTGCGCGGCGTAATAGTTGGTCGGGCATTACTGGAAGGTAAATTc     >  1:503330/1‑140 (MQ=255)
                                                                                                                   gaagaaGTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCGCGGCGTAATAGTTGGTCGGGCATTACTGGAAGGTAAATTCAc   <  2:75130/140‑1 (MQ=255)
                                                                                                                    aagaagTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCGCGGCGTAATAGTTGGTCGGGCATTACTGGAAGGTAAATTCAcc  <  2:539123/140‑1 (MQ=255)
                                                                                                                    aagaagTGCGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCGCGGCGTAATAGTTGGTCGGGCATTACTGGAAGGTAAATTCAcc  <  2:674666/140‑1 (MQ=255)
                                                                                                                            |                                                                                                                                   
GGGCGTTTCACTGGAACAACTGGTGGAAACCTATCTGCCCGTCGGCCTGAAACATGTGCTGTGTACCGATATCTCGCGCGACGGCACGCTGGCAGGCTCTAACGTCTCTTTATATGAAGAAGTGTGCGCCAGATATCCGCAGGTGGCATTTCAGTCCTCCGGCGGTATTGGCGACATTGATGATGTGGCGGCCCTGCGTGGCACTGGTGTGCGCGGCGTAATAGTTGGTCGGGCATTACTGGAAGGTAAATTCACC  >  NC_000913/2095571‑2095826

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