Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,287,826 C→T A267T (GCA→ACA)  purU ← formyltetrahydrofolate hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,287,8260CT98.4% 216.7 / ‑5.1 64A267T (GCA→ACA) purUformyltetrahydrofolate hydrolase
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (26/37);  total (27/37)
Fisher's exact test for biased strand distribution p-value = 4.22e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

GCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCC  >  NC_000913/1287690‑1287953
                                                                                                                                        |                                                                                                                               
gccgcTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCa                                                                                                                               <  1:823209/139‑1 (MQ=255)
gccgcTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCa                                                                                                                               <  2:256158/139‑1 (MQ=255)
gccgcTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCa                                                                                                                               <  1:253948/139‑1 (MQ=255)
gccgcTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCa                                                                                                                               <  1:505880/139‑1 (MQ=255)
  cgcTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGt                                                                                                                             <  2:177024/139‑1 (MQ=255)
  cgcTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGt                                                                                                                             <  2:74636/139‑1 (MQ=255)
     tGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACt                                                                                                                          <  2:929491/139‑1 (MQ=255)
     tGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACt                                                                                                                          <  1:596495/139‑1 (MQ=255)
          aGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGta                                                                                                                     <  2:117514/139‑1 (MQ=255)
            tgtTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGtata                                                                                                                   >  2:175685/1‑139 (MQ=255)
                 cGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACt                                                                                                                          >  2:1082092/1‑127 (MQ=255)
                 cGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACt                                                                                                                          <  1:1082092/127‑1 (MQ=255)
                 cGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCa                                                                                                              >  2:339570/1‑139 (MQ=255)
                      agaAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGc                                                                                                               >  2:762234/1‑133 (MQ=255)
                      agaAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGc                                                                                                               <  1:762234/133‑1 (MQ=255)
                      agaAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACt                                                                                                         >  2:1107519/1‑139 (MQ=255)
                             tGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGt                                                                                                                             <  1:196434/112‑1 (MQ=255)
                             tGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGt                                                                                                                             >  2:196434/1‑112 (MQ=255)
                                aCTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCa                                                                                                                               <  2:710513/107‑1 (MQ=255)
                                       gCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGa                                                                                        >  2:695862/1‑139 (MQ=255)
                                                  aTCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGtt                                                                                                <  2:469732/120‑1 (MQ=255)
                                                  aTCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGtt                                                                                                >  1:469732/1‑120 (MQ=255)
                                                      gTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTgcgc                                                                         <  2:766418/139‑1 (MQ=255)
                                                        aaCGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTgcgcgc                                                                       <  1:1093675/139‑1 (MQ=255)
                                                              ttAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCata                                                                 >  2:981013/1‑139 (MQ=255)
                                                                     tAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                          <  1:1107519/139‑1 (MQ=255)
                                                                     tAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                          <  1:585561/139‑1 (MQ=255)
                                                                     tAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                          <  2:322750/139‑1 (MQ=255)
                                                                     tAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                          <  1:807524/139‑1 (MQ=255)
                                                                     tAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGGATAGTGCATGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                          <  1:1049700/139‑1 (MQ=255)
                                                                     tAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCACTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                          <  2:268716/139‑1 (MQ=255)
                                                                       aCCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTa                                                                                                                            >  2:769043/1‑71 (MQ=255)
                                                                       aCCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTa                                                                                                                            <  1:769043/71‑1 (MQ=255)
                                                                           gTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTgtgt                                                    >  2:887856/1‑139 (MQ=255)
                                                                             tATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCa                                                                                                                               >  2:625601/1‑62 (MQ=255)
                                                                             tATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                          <  2:438894/131‑1 (MQ=255)
                                                                             tATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                          >  1:438894/1‑131 (MQ=255)
                                                                             tATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAg                                                  >  1:282617/1‑139 (MQ=255)
                                                                                   tGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGa                                            <  2:316858/139‑1 (MQ=255)
                                                                                   tGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGa                                            <  2:430401/139‑1 (MQ=255)
                                                                                      tGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCTCGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                          <  2:748965/122‑1 (MQ=255)
                                                                                      tGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                          >  1:748965/1‑122 (MQ=255)
                                                                                       ggCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGa                                        >  1:1035560/1‑139 (MQ=255)
                                                                                       ggCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGa                                        >  1:62408/1‑139 (MQ=255)
                                                                                                    taatCGTTCGATTACCGTCAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTc                           <  2:747678/139‑1 (MQ=255)
                                                                                                      atCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACg                                                                                                  >  2:795502/1‑66 (MQ=255)
                                                                                                      atCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACg                                                                                                  <  1:795502/66‑1 (MQ=255)
                                                                                                      atCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                          <  2:827327/106‑1 (MQ=255)
                                                                                                      atCGTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAg                                                          >  1:827327/1‑106 (MQ=255)
                                                                                                         gTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCa                      >  2:637446/1‑139 (MQ=255)
                                                                                                         gTTCGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCa                      <  2:282617/139‑1 (MQ=255)
                                                                                                            cGATTACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCatga                   >  1:134217/1‑139 (MQ=255)
                                                                                                               ttACCGTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCAtgatga                >  1:99456/1‑139 (MQ=255)
                                                                                                                    gTAAACAAAGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTggg           <  2:861219/139‑1 (MQ=255)
                                                                                                                           aaGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATGGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGt    <  1:866928/139‑1 (MQ=255)
                                                                                                                            aGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTgcgc                                                                         >  2:74373/1‑69 (MQ=255)
                                                                                                                            aGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTgcgc                                                                         <  1:74373/69‑1 (MQ=255)
                                                                                                                            aGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCtgt                                                      >  1:1098991/1‑88 (MQ=255)
                                                                                                                            aGACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCtgt                                                      <  2:1098991/88‑1 (MQ=255)
                                                                                                                             gACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTtct                                                                                           <  2:839356/50‑1 (MQ=255)
                                                                                                                             gACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTtct                                                                                           >  1:839356/1‑50 (MQ=255)
                                                                                                                             gACGCGCTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTcc  >  1:455255/1‑139 (MQ=255)
                                                                                                                                 cgcTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGAt                                           >  2:885813/1‑94 (MQ=255)
                                                                                                                                 cgcTGTGTCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGAt                                           <  1:885813/94‑1 (MQ=255)
                                                                                                                                        |                                                                                                                               
GCCGCTGTAAAGTGTTACGTTGAGAAAAATGAACTGGTTGCGTAATTTTCATCCGTAACGGATTAAAGGTAACCAGTTATTTTTGCTGGCGATTAAAGAATAATCGTTCGATTACCGTAAACAAAGACGCGCTGTGCCAGTACTTTGTATAGTGCACGACTTAAGACGTTTTTCTCGACGTCACGACCTGCGCGCATCATATCTTCAGCTGTGTAGGTATGATCGACATGAATAACGTCCTGCATGATGATTGGGCCTTCGTCC  >  NC_000913/1287690‑1287953

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