Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,931,976 (G)5→6 coding (1096/1179 nt) purT → phosphoribosylglycinamide formyltransferase 2

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,931,9711.G95.6% 143.7 / ‑2.4 45L364R (CTG→CGG) purTphosphoribosylglycinamide formyltransferase 2
Reads supporting (aligned to +/- strand):  ref base . (0/0);  major base G (19/24);  minor base A (1/1);  total (20/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.22e-01

TCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCC  >  NC_000913/1931837‑1932088
                                                                                                                                       |                                                                                                                     
tCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTgggg                                                                                                                   >  2:1037719/1‑139 (MQ=255)
     cAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTggggggtgg                                                                                                              >  1:203473/1‑135 (MQ=255)
        cTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCAc                                                                                                            >  1:362555/1‑139 (MQ=255)
          tCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTg                                                                                                          >  1:373398/1‑139 (MQ=255)
          tCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTggg                                                                                                                    >  1:107469/1‑128 (MQ=255)
          tCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTggg                                                                                                                    <  2:107469/128‑1 (MQ=255)
           cTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTgg                                                                                                         <  1:534175/139‑1 (MQ=255)
             gCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCt                                                                                                       <  2:362555/139‑1 (MQ=255)
                   aTTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCa                                                                                                 <  1:131/139‑1 (MQ=255)
                              aCTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGtgttgt                                                                                      <  1:266489/139‑1 (MQ=255)
                                   ccAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACt                                                                                                           <  2:665434/113‑1 (MQ=255)
                                   ccAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACt                                                                                                           >  1:665434/1‑113 (MQ=255)
                                   ccAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGc                                                                                                        <  2:407117/116‑1 (MQ=255)
                                   ccAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGc                                                                                                        >  1:407117/1‑116 (MQ=255)
                                   ccAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACg                                                                                 >  2:577526/1‑139 (MQ=255)
                                      gTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCa                                                                                                             <  1:269334/108‑1 (MQ=255)
                                      gTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCa                                                                                                             >  2:269334/1‑108 (MQ=255)
                                            aTGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGa                                                                                   <  2:557812/128‑1 (MQ=255)
                                            aTGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGa                                                                                   >  1:557812/1‑128 (MQ=255)
                                                      tGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCa                                                              <  2:203473/139‑1 (MQ=255)
                                                              tgCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGccgccg                                                      >  1:545223/1‑139 (MQ=255)
                                                              tgCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTgggggg                                                                                                                 <  2:306450/79‑2 (MQ=255)
                                                              tgCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTgggggg                                                                                                                 >  1:306450/1‑78 (MQ=255)
                                                               gCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGttgttg                                                                                    <  2:934297/108‑1 (MQ=255)
                                                               gCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGttgttg                                                                                    >  1:934297/1‑108 (MQ=255)
                                                               gCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCgg                                                     <  1:1037719/139‑1 (MQ=255)
                                                                cAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCa                                                                                                 >  1:625166/1‑94 (MQ=255)
                                                                cAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCa                                                                                                 <  2:625166/94‑1 (MQ=255)
                                                                cAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGa                                                    <  2:373398/139‑1 (MQ=255)
                                                                 aGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGc                                                                                                  >  2:99658/1‑92 (MQ=255)
                                                                 aGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGc                                                                                                  <  1:99658/92‑1 (MQ=255)
                                                                  gAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACa                                                  >  1:220566/1‑139 (MQ=255)
                                                                   aaTGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAg                                                 >  2:198637/1‑139 (MQ=255)
                                                                      gCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTgg                                                                                                         <  1:344931/80‑1 (MQ=255)
                                                                      gCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTgg                                                                                                         >  2:344931/1‑80 (MQ=255)
                                                                        cGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTa                                              <  2:808557/137‑1 (MQ=255)
                                                                        cGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTa                                              >  1:808557/1‑137 (MQ=255)
                                                                          tAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAg                                          >  2:283122/1‑139 (MQ=255)
                                                                            ggCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTa                                        <  2:1211664/139‑1 (MQ=255)
                                                                            ggCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTa                                        <  2:220566/139‑1 (MQ=255)
                                                                            ggCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTa                                        <  1:577526/139‑1 (MQ=255)
                                                                             gcgcAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTAc                                       <  2:545223/139‑1 (MQ=255)
                                                                                                       ggTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCTGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCt             >  2:1430026/1‑139 (MQ=255)
                                                                                                       ggTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACa                                      >  2:267321/1‑114 (MQ=255)
                                                                                                       ggTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACa                                      <  1:267321/114‑1 (MQ=255)
                                                                                                       ggTAAGCCGGAAATTGATGGCAGCCGTCGTCTgggggg                                                                                                                 <  2:211145/38‑2 (MQ=255)
                                                                                                       ggTAAGCCGGAAATTGATGGCAGCCGTCGTCTgggggg                                                                                                                 >  1:211145/1‑37 (MQ=255)
                                                                                                         tAAGCCGGAAATTGATGGCAGCCGTCGTCTAGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACa                                                  >  1:834519/1‑100 (MQ=255)
                                                                                                         tAAGCCGGAAATTGATGGCAGCCGTCGTCTAGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACa                                                  <  2:834519/100‑1 (MQ=255)
                                                                                                                  aaTTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCTGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGcc  <  1:1430026/139‑1 (MQ=255)
                                                                                                                  aaTTGATGGCAGCCGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGcc  <  1:198637/139‑1 (MQ=255)
                                                                                                                              ccGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAAcgcg                                                                     <  2:1214465/60‑1 (MQ=255)
                                                                                                                              ccGTCGTCTGGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAAcgcg                                                                     >  1:1214465/1‑60 (MQ=255)
                                                                                                                                       |                                                                                                                     
TCCTGCAGCTTCTGCCGTTATTCTGCCACAACTGACCAGTCAGAATGTCACGTTTGATAATGTGCAGAATGCCGTAGGCGCAGATTTGCAGATTCGTTTATTTGGTAAGCCGGAAATTGATGGCAGCCGTCGTCTGGGGGTGGCACTGGCTACTGCAGAGAGTGTTGTTGACGCCATTGAACGCGCGAAGCACGCCGCCGGACAGGTAAAAGTACAGGGTTAAACCCGGCAAAAAAACGCTACAAAAATGCC  >  NC_000913/1931837‑1932088

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