Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,485,150 A→G 19.5% intergenic (+56/‑91) 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,1500AG19.5% 76.8 / 15.2 41intergenic (+56/‑91)purH/purDbifunctional phosphoribosylaminoimidazolecarboxamide formyltransferase/IMP cyclohydrolase/phosphoribosylamine‑‑glycine ligase
Reads supporting (aligned to +/- strand):  ref base A (10/23);  new base G (2/6);  total (12/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-03

GGTGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTC  >  NC_002947/5485002‑5485297
                                                                                                                                                         |                                                                                                                                                   
ggtggtTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATa                                                                                                                                                     <  8:248429/149‑1 (MQ=255)
     ttCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtg                                                                                                                                                <  3:5090/149‑1 (MQ=255)
     ttCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGAtg                                                                                                                                                <  7:298673/149‑1 (MQ=255)
     ttCGATGCGTGATGCTGAAGTCATCGCTGCTACCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCggcg                                                                                                                                                >  5:96456/1‑146 (MQ=255)
             gTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGAtt                                                                                                                                         >  5:331213/1‑148 (MQ=255)
                   cTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCAccgcc                                                                                                                                  <  4:323935/149‑1 (MQ=255)
                   cTGAAGTCATCGCTGCTACCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCAccatc                                                                                                                                  <  6:96456/149‑4 (MQ=38)
                      aaGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGAt                                                                                                                                          <  4:219963/138‑1 (MQ=255)
                      aaGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGAt                                                                                                                                          >  3:219963/1‑138 (MQ=255)
                           aTCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc                                                                                                                          <  1:75130/149‑1 (MQ=37)
                            tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCTTAAAAAGATGCTCGAAGTTCCTGTGGGAGCGGGCTTTCCCGCGAATAAGATGGTGGATTAACcgccgcct                                                                                                                              >  1:311442/1‑146 (MQ=255)
                            tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGTGGCTTGCGCCGAACACTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgc                                                                                                                          <  3:11084/149‑1 (MQ=14)
                            tCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTcgcg                                                                                                                         <  6:83360/149‑1 (MQ=37)
                              gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCTCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                        <  2:180772/148‑1 (MQ=255)
                              gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                        <  7:301339/148‑1 (MQ=37)
                              gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATGCGGCGGTGAATCCACCATCGTATTCGCgg                                                                                                                        <  1:349425/148‑1 (MQ=37)
                              gctgctGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCgg                                                                                                                        <  7:53711/148‑1 (MQ=255)
                                   tgcCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCACCGCATTCGCGGGCAAGc                                                                                                                  <  6:319115/149‑1 (MQ=255)
                                      cgacgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCg                                                                                                               <  8:307241/149‑1 (MQ=255)
                                         cgaGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTc                                                                                                            <  6:331213/149‑1 (MQ=255)
                                                         aTGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGGTTGCCCGCCCCATTCGCGGGCCAGCCCGCTCCCCCAGTGTTCGgcgc                                                                                            >  7:183718/1‑149 (MQ=255)
                                                               ttCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTTCCCCCGGATACGATGGGGGGTTGACCGCCGCATTTGCGGGGAAGCCCGCGGCCACCGGGTTCGGCgccagcca                                                                                      >  3:130080/1‑149 (MQ=255)
                                                                  agtGGCTTGCGCCATTTCCGCCACTAATTCGGCGTATCCCGTGGAATGTCTAGATCTTGGGGGCGGGGGTTTGCCCGCTAATACTTTGTTGTATTCCCCGCCGCATTGGGGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACgg                                                                                   <  8:33656/147‑1 (MQ=255)
                                                                  aggggCATGCGCAACTTCCGCCATAAATTCCGCGTACCCCCTGCAACGTCGATATCGTGGGGGCGGGGGCTTGCCGGCTAATACTTGGGGGTATTCACCGCCGCTTTCGGGGGAATGCCCGCTCCCACAGTGTTCGGCGCAAGCCACg                                                                                    <  2:21589/145‑1 (MQ=255)
                                                                            gcCACTTCCGCCACTAATTACGCGTATCCCTTGGCACGTCGAGATCCTGTGGGAGCGGGGTGGCCGGCTAATACGATGGTGGATTCACCGCCGCATTGGCGGGCAAGCCCGCTCCcaca                                                                                                       <  6:311975/119‑1 (MQ=37)
                                                                            gcCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTTCCCGCGCATAAGCTGGTGGGTTGACCGCCGCCTTAGCGGGCCAGCCCGCGCCccca                                                                                                       >  5:311975/1‑119 (MQ=37)
                                                                                   ccGCCACTAATTCCGCGGATCCCGCGGCAAGTCGAGATCCTGGGGGAGCGGGCTTGCCCGCTAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAg                                                                   <  2:311442/148‑1 (MQ=255)
                                                                                      ccACTAATTACGCGGATCCCGAGCCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACt                                                                            <  7:86744/136‑1 (MQ=255)
                                                                                                            gtcAAGTCGATATCCTGTGGGAGGGGGCTTCCGCGCAAATACTATTGTGGATTCCCCGCCTCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGatcatc                                         <  3:241221/147‑1 (MQ=255)
                                                                                                                     tCGAGATCCTGTGGGAGCGGGCTGGCCCGCGAATACGATGGTGGATTCACCGCCACATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGc                                   <  8:183718/149‑1 (MQ=255)
                                                                                                                                  tgGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGTTCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGa                        <  4:284040/149‑1 (MQ=255)
                                                                                                                                     gAGCGGGATGGCCCGCTAATACGATGGTGGTTTCACCGCCGCATTGGCGGGAACGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCa                     <  7:212335/149‑1 (MQ=255)
                                                                                                                                         gggCTTGCCCGCTAAAACGATGGGGGATTCCCCGCCGCATTCGCGGGAAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGcc                 <  7:156053/149‑1 (MQ=255)
                                                                                                                                             ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGa                                              >  8:348054/1‑116 (MQ=255)
                                                                                                                                             ttGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGa                                              <  7:348054/116‑1 (MQ=255)
                                                                                                                                                cgacagAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcct          <  1:324562/144‑1 (MQ=255)
                                                                                                                                                 ccGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg         >  4:174702/1‑149 (MQ=255)
                                                                                                                                                  cgcgAATGCGGCGGTGAATCCACCATCGTATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCcctggcctg         >  7:323718/1‑148 (MQ=37)
                                                                                                                                                      aaTACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAg    <  3:174702/149‑1 (MQ=255)
                                                                                                                                                        tACGATGGTTGATTCACCGCCTCTTTCGCGGGCAAGCCCGCTCCCACATTGTTCGGCGCAAGCCACGGGATCACTTGAATCTAGGTTTTGACATTAAAGTTTTGATCATCg                                        <  4:16913/111‑1 (MQ=255)
                                                                                                                                                        tACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCg                                        >  3:16913/1‑111 (MQ=255)
                                                                                                                                                        tACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTc  >  7:29377/1‑149 (MQ=255)
                                                                                                                                                         |                                                                                                                                                   
GGTGGTTCGATGCGTGATGCTGAAGTCATCGCTGCTGCCGACGAGGCCGGCATCGCGATGGTCTTCACTGGCATGCGCCACTTCCGCCACTAATTACGCGGATCCCGAGGCAAGTCGAGATCCTGTGGGAGCGGGCTTGCCCGCGAATACGATGGTGGATTCACCGCCGCATTCGCGGGCAAGCCCGCTCCCACAGTGTTCGGCGCAAGCCACGGGATCACTTGAATCGAGGTTTTGACATGAAAGTTTTGATCATCGGCAGCGGTGGCCGTGAGCACGCCCTGGCCTGGAAAGTC  >  NC_002947/5485002‑5485297

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