Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,485,154 T→C 19.6% intergenic (+60/‑87) purH → / → purD bifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,485,1540TC19.6% 61.0 / 12.3 38intergenic (+60/‑87)purH/purDbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base T (10/19);  new base C (4/3);  total (16/22)
Fisher's exact test for biased strand distribution p-value = 3.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.09e-01

TTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCG  >  NC_002947/5485007‑5485298
                                                                                                                                                   |                                                                                                                                                
ttCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGAGGGAGCGGGCTTGCCCGCGAATACGAtg                                                                                                                                                 <  1:185049/149‑1 (MQ=255)
    aTGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggtgg                                                                                                                                             <  1:171254/149‑1 (MQ=255)
      gCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAAt                                                                                                                                           >  3:29180/1‑149 (MQ=255)
         tGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggt                                                                                                                                               >  3:138078/1‑142 (MQ=255)
         tGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtggt                                                                                                                                               <  4:138078/142‑1 (MQ=255)
            tGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccg                                                                                                                                     >  7:163428/1‑149 (MQ=255)
                 aaGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACcgccgc                                                                                                                                 <  8:163428/148‑1 (MQ=255)
                   gTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCAtt                                                                                                                              >  5:148830/1‑149 (MQ=255)
                      aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGGCCGCGCATACGAAGGTTGCATCACCCCCGCATTctc                                                                                                                           >  8:167259/1‑147 (MQ=255)
                      aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc                                                                                                                           >  1:72998/1‑149 (MQ=37)
                       tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgcg                                                                                                                          >  7:259962/1‑149 (MQ=37)
                         gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCTGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                         <  8:259962/148‑1 (MQ=25)
                         gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCggg                                                                                                                        <  2:72998/149‑1 (MQ=25)
                         gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCggg                                                                                                                        <  6:148830/149‑1 (MQ=255)
                         gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                         <  6:66676/148‑1 (MQ=255)
                          ctgctgCCGACGAGGCCGGCATCTCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCGGGc                                                                                                                       <  4:29180/149‑1 (MQ=25)
                                  gacgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGCATAAGATGGTGGATTCACCGCCGCATTAGCGGGCCAGCCCGc                                                                                                               >  7:72773/1‑149 (MQ=255)
                                          cGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAAATCCTGTGGGAGCGGGCCTGGCCGCGAATAACATGGTGGGTTAAACGCCCCAATAGCGGGCAACCCCGCTACcacaa                                                                                                       >  8:174040/1‑148 (MQ=255)
                                                               tGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCATTCCCGCCAAAAAGAAGGTGGATTAACCCCCGCAATAGCCCGCAACCACGCTACCCCACTGTTAGGCCCCAGCCACGGGa                                                                                  >  4:225540/1‑149 (MQ=255)
                                                                   aTGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCGTGCCCGCGGATAAGATGGGGGGTTGACCGCCGCATTCGCCGGCAAGCCCGCGCCCACAGGGTTCGGCGCAAGCCACGGGATGAc                                                                              >  6:296637/1‑149 (MQ=255)
                                                                     gcccaaCTTCCGCCACAAATTACGCGGATCCCGCCGCATGCGGACACCCTGGGGCTGCGGGCGGGCCCGCTAATACCATGGTCGATCCCCCCCCGCATGCGGGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACtt                                                                            <  8:189955/144‑1 (MQ=255)
                                                                      cgcCACTTCCGCCACTAATTACTCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                         <  2:213332/103‑1 (MQ=255)
                                                                      cgcCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                         >  1:213332/1‑103 (MQ=255)
                                                                              ccGCCACTAATTACGCGTATCCCGAGGCAAGTCTCGATCCTGGGGGAGCGGGCTTGCCCGCGAATACGATGGGGGATTCACCGCCGCATTCCCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAg                                                                    <  5:135027/148‑1 (MQ=38)
                                                                               cGCCACTAATTCCGCGGATCCCGAGGCACGGCGAGAGCCTGTGGGAGCGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCaa                                                                                           <  1:209992/124‑1 (MQ=37)
                                                                                gCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCACGCGAATAAGATGGTGGGTTCACCGCCGCCTTAGCGGGCAAGCCCCCGCCCACACTGGTCGGCGCAAGCCCCCGGATCACGTGAATGGAgggt                                                                 >  4:239799/1‑147 (MQ=255)
                                                                                    cTAATTCCGCGGATCCCGAGGCAAGGCGAGATCCTGTGGGAGCGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGa                                                             <  1:207298/149‑1 (MQ=255)
                                                                                     tAATTCCGCGTATCCCGTGGCATGGCGAGACCCTGTGGTAGGGGGCTGGCCCGCTAATACTATGGTGTATTCCCCGCCCCATGCGGGGCAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGAc                                                            <  4:271087/149‑1 (MQ=255)
                                                                                           cgcgGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGaaa                                                      <  3:225540/149‑1 (MQ=255)
                                                                                                 tCCCGTGGCAAGTCGAGATCGTGTGGGAGCGGGCTTGCCCGCTAATACTATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGtttt                                                 <  8:72773/148‑1 (MQ=255)
                                                                                                   cctaGGAAAGTCGAGATCGTGTGGGAGCGGGCTTCCCGGCTAATACGATGGTGTATTCACCGCCGAATTCGCGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGa                                               <  7:167259/145‑1 (MQ=255)
                                                                                                   ccGAGGCAAGTCGAGATCCTGTGGGAGCCGGCTTGCCCGCTAATACGATGGGGTATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGa                                               <  6:302412/148‑1 (MQ=255)
                                                                                                       ggCAAGTCGATATCCTGGGGCAGCGGGCTTCCCGGACAATACTATGTTGAATTCACCGCCCAATTCGCGGTAATGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatcatc                                          <  3:51429/149‑1 (MQ=255)
                                                                                                                           gcgCGGGGTTGCCGGCTAATACGTTGGTGGATTCACCGCCGCATGCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                         <  3:129838/144‑1 (MQ=255)
                                                                                                                                ggCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTg                                                <  7:19575/118‑1 (MQ=255)
                                                                                                                                ggCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTg                                                >  8:19575/1‑118 (MQ=255)
                                                                                                                                ggCTTGCCCGCCAATACGATGGTGCATTCACCGCCCCATTCCCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGccc                 <  3:22785/149‑1 (MQ=255)
                                                                                                                                   ttGCCCGCGAATACGATGGTGGATTCACAGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctgg              >  7:176724/1‑149 (MQ=255)
                                                                                                                                      cccGCTAATACGATGGTGGATTCACCGCCGCTTTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct           <  5:327752/149‑1 (MQ=255)
                                                                                                                                        cgcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg          >  3:304334/1‑148 (MQ=37)
                                                                                                                                        cgcgAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg          >  8:237297/1‑148 (MQ=255)
                                                                                                                                               ccGATGGGGCATGCCCCGCCGCATCCGCGGCCACCCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCg  <  3:45024/148‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
TTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTCG  >  NC_002947/5485007‑5485298

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