Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,187,241 C→T 8.1% Q816* (CAA→TAA)  purL → phosphoribosylformylglycinamidine synthase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,187,2410CT 8.1% 139.3 / 2.4 49Q816* (CAA→TAA) purLphosphoribosylformylglycinamidine synthase
Reads supporting (aligned to +/- strand):  ref base C (20/25);  new base T (2/2);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.96e-01

CATTCCGGTCGGCAAAGACTCGATGTCGATGAAGACCAAGTGGAGCGACGAGGGCGTCGAGAAGAGCGTCACCTCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAACCGCATGGGCGCCTCGATCCTGGCCCAGACCTACGGCAAGATCGCGGCCCAGGCACCGGACGTCGACGACGCCGAAGACCTCAAGGC  >  NC_002947/1187102‑1187389
                                                                                                                                           |                                                                                                                                                    
cATTCCGGTCGGCAAAGACTCGATGTCGATGAAGACCAAGTGGAGCGACGAGGGCGTCGAGAAGAGCGTCACCTCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCa                                                                                                                                             <  3:296261/149‑1 (MQ=255)
 aTTCCGGTCGGCAAAGACTCGATGTCGATGAAGACCAAGTGGAGCGACGAGGGCGTCGAGAAGAGCGTCACCTCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCAt                                                                                                                                            <  6:158056/149‑1 (MQ=255)
  ttCCGGTCGGCAAAGACTCGATGTCGATGAAGACCAAGTGGAGCGACGAGGGCGTCGAGAAGAGCGTCACCTCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGTAACTGCGCATg                                                                                                                                           >  2:120753/1‑149 (MQ=255)
  ttCCGGTCGGCAAAGACTCGATGTCGATGAAGACCAAGTGGAGCGACGAGGGCGTCGAGAAGAGCGTCACCTCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATg                                                                                                                                           <  3:233194/149‑1 (MQ=255)
    ccGGTCGGCAAAGACTCGATGTCGATGAAGACCAAGTGGAGCGACGAGGGCGTCGAGAAGAGCGTCACCTCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATgg                                                                                                                                          <  2:181578/148‑1 (MQ=255)
      ggTCGGCAAAGACTCGATGTCGATGAAGACCAAGTGGAGCGACGAGGGCGTCGAGAAGAGCGTCACCTCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGAc                                                                                                                                        >  3:219058/1‑148 (MQ=255)
        tCGGCAAAGACTCGATGTCGATGAAGACCAAGTGGAGCGACGAGGGCGTCGAGAAGAGCGTCACCTCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAg                                                                                                                                     >  7:22995/1‑149 (MQ=255)
              aaGACTCGATGTCGATGAAGACCAAGTGGAGCGACGAGGGCGTCGAGAAGAGCGTCACCTCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCTCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCgag                                                                                                                               <  4:219058/149‑1 (MQ=255)
              aaGACTCGATGTCGATGAAGACCAAGTGGAGCGACGAGGGCGTCGAGAAGAGCGTCACCTCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCgag                                                                                                                               <  3:147942/149‑1 (MQ=255)
               aGACTCGATGTCGATGAAGACCAAGTGGAGCGACGAGGGCGTCGAGAAGAGCGTCACCTCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCgaga                                                                                                                              <  6:133241/149‑1 (MQ=255)
               aGACTCGATGTCGATGAAGACCAAGTGGAGCGACGAGGGCGTCGAGAAGAGCGTCACCTCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCgaga                                                                                                                              >  3:10524/1‑149 (MQ=255)
                     gatgtcgatgAAGACCAAGTGGAGCGACGAGGGCGTCGAGAAGAGCGTCACCTCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAgaccgac                                                                                                                         <  1:140719/148‑1 (MQ=255)
                     gatgtcgatgAAGACCAAGTGGAGCGACGAGGGCGTCGAGAAGAGCGTCACCTCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAgaccgac                                                                                                                         >  2:140719/1‑148 (MQ=255)
                         tcgatgAAGACCAAGTGGAGCGACGAGGGCGTCGAGAAGAGCGTCACCTCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGAcctga                                                                                                                     <  4:10524/148‑1 (MQ=255)
                           gatgAAGACCAAGTGGAGCGACGAGGGCGTCGAGAAGAGCGTCACCTCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGTAACTgcgc                                                                                                                                              >  6:283609/1‑121 (MQ=255)
                           gatgAAGACCAAGTGGAGCGACGAGGGCGTCGAGAAGAGCGTCACCTCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGTAACTgcgc                                                                                                                                              <  5:283609/121‑1 (MQ=255)
                           gatgAAGACCAAGTGGAGCGACGAGGGCGTCGAGAAGAGCGTCACCTCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGAcctgatc                                                                                                                   <  3:136714/148‑1 (MQ=255)
                                 gACCAAGTGGAGCGACGAGGGCGTCGAGAAGAGCGTCACCTCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCg                                                                                                            <  2:242626/149‑1 (MQ=255)
                                   ccAAGTGGAGCGACGAGGGCGTCGAGAAGAGCGTCACCTCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGAc                                                                                                          >  1:55834/1‑149 (MQ=255)
                                          gAGCGACGAGGGCGTCGAGAAGAGCGTCACCTCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAgaccgacc                                                                                                                        <  5:324251/128‑1 (MQ=255)
                                          gAGCGACGAGGGCGTCGAGAAGAGCGTCACCTCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAgaccgacc                                                                                                                        >  6:324251/1‑128 (MQ=255)
                                             cgacgaGGGCGTCGAGAAGAGCGTCACCTCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCcgcg                                                                                                <  1:312361/149‑1 (MQ=255)
                                                    ggCGTCGAGAAGAGCGTCACCTCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATTGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCaaaaa                                                                                         <  8:22995/149‑1 (MQ=255)
                                                                         tCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAAcc                                                                                       <  6:326138/130‑1 (MQ=255)
                                                                         tCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAAcc                                                                                       >  5:326138/1‑130 (MQ=255)
                                                                         tCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAACCGCATGGGCGCCTCGATCCt                                                                    >  5:114204/1‑149 (MQ=255)
                                                                          cGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAACCGCATGGGCGCCTCGATCCTg                                                                   <  3:68767/149‑1 (MQ=255)
                                                                           gCCGATGTCGCTGATCATCGCCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGTCCTGATCCTGATCGACCTGGGCCGCGGCAAAAACCGCATGGGCGCCTCGATCCTgg                                                                  <  5:55500/149‑1 (MQ=255)
                                                                                  tCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGTAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAACCGCATGGGCGCCTCGATCCTGGCCCAGAc                                                           <  1:120753/149‑1 (MQ=255)
                                                                                    gCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAACCGCATgg                                                                                 <  3:44360/125‑1 (MQ=255)
                                                                                    gCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAACCGCATgg                                                                                 >  4:44360/1‑125 (MQ=255)
                                                                                      tGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAACCGCATGGGCGCCTCGATCCTGGCCCAGACCTAc                                                       >  4:132239/1‑149 (MQ=255)
                                                                                        atcatcACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAACCGCATGGGCGCCTCGATCCTGGCCCAGACCTAc                                                       <  1:254767/147‑1 (MQ=255)
                                                                                           atcaCCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAACCGCATGGGCGCCTCGATCCTGGCCCAGACCTACGGc                                                    >  5:69321/1‑147 (MQ=255)
                                                                                               ccGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAACCGCATGGGCGCCTCGATCCTGGCCCAGACCt                                                         <  8:187033/138‑1 (MQ=255)
                                                                                               ccGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAACCGCATGGGCGCCTCGATCCTGGCCCAGACCt                                                         >  7:187033/1‑138 (MQ=255)
                                                                                               ccGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAACCGCATGGGCGCCTCGATCCTGGCCCAGACCTACGGc                                                    >  5:305777/1‑143 (MQ=255)
                                                                                                    ttCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAACCGCATGGGCGCCTCGATCCTGGCCCAGACCTACGGCAAGATCGCGGc                                         >  2:27120/1‑149 (MQ=255)
                                                                                                       gCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAACCGCATGGGCGCCTCGATCCTGGCCCAGACCTACGGCAAGATCGCGGccc                                       <  3:132239/148‑1 (MQ=255)
                                                                                                         cccGGTCACCGACATCCGCAAGACCCTGAGCCCTCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTTGGCCGCGGCAAAAACCGCa                                                                                    <  4:45408/101‑1 (MQ=255)
                                                                                                         cccGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGAATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAACCGCa                                                                                    >  3:45408/1‑101 (MQ=255)
                                                                                                          ccGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAACCGCATGGGCGCCTCGATCCTGGCCCAGACCTACGGCAAGATCGCGGCCCAgg                                    >  8:124915/1‑148 (MQ=255)
                                                                                                              tCACCGACATCCGCAAGAGCCTGAGCCCGCAACTGCGCATCGACTAGGTCGAGACGGAGGTTTCCGTGTTCGACCTGGGCCGCGGCAAAAACCGCATGGGCGGGTCGATCATGGACCATAACTAGGTCCAGATCGGGGCCCATGCAGCg                               <  6:101492/149‑1 (MQ=255)
                                                                                                                   gACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAACCGCATGGGCGCCTCGATCCTGGCCCAGACCTACGGCAAGATCGCGGCCCAGGCACCGGACGt                          >  4:21664/1‑149 (MQ=255)
                                                                                                                   gACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAACCGCATGGGCGCCTCGATCCTGGCCCAGACCTACGGCAAGATCGCGGCCCAGGCACCGGACGt                          <  3:150884/149‑1 (MQ=255)
                                                                                                                           cAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAACCGCATGGGCGCCTCGATCCTGGCCCAGACCTACGGCAAGATCGCGGCCCAGGCACCGGACGTCGacgacg                  >  3:85781/1‑149 (MQ=255)
                                                                                                                            aaGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAACCGCATGGGCGCCTCGATCCTGGCCCAGACCTACGGCAAGATCGCGGCCCAGGCACCGGACGTCGACGACGc                 <  6:69321/149‑1 (MQ=255)
                                                                                                                                 ccTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAACCGCATGGGCGCCTCGATCCTGGCCCAGACCTACGGCAAGATCGCGGCCCAGGCACCGGACGTCGACGACGCCg               >  2:14055/1‑146 (MQ=255)
                                                                                                                                          gCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAACCGCATGGGCGCCTCGATCCTGGCCCAGACCTACGGCAAGATCGCGGCCCAGGCACCGGACGTCGACGACGCCGAAGACCTCAAgg   >  1:300868/1‑149 (MQ=255)
                                                                                                                                           cAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAACCGCATGGGCGCCTCGATCCTGGCCCAGACCTACGGCAAGATCGCGGCCCAGGCACCGGACGTCGACGACGCCGAAGACCTCAAgcc  >  5:191422/1‑147 (MQ=255)
                                                                                                                                           |                                                                                                                                                    
CATTCCGGTCGGCAAAGACTCGATGTCGATGAAGACCAAGTGGAGCGACGAGGGCGTCGAGAAGAGCGTCACCTCGCCGATGTCGCTGATCATCACCGGCTTCGCCCCGGTCACCGACATCCGCAAGACCCTGACCCCGCAACTGCGCATGGACAAGGGCGAGACCGACCTGATCCTGATCGACCTGGGCCGCGGCAAAAACCGCATGGGCGCCTCGATCCTGGCCCAGACCTACGGCAAGATCGCGGCCCAGGCACCGGACGTCGACGACGCCGAAGACCTCAAGGC  >  NC_002947/1187102‑1187389

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