Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,206,768 A→G V255V (GTT→GTC)  purH ← bifunctional AICAR transformylase/IMP cyclohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,206,7680AG98.0% 173.7 / NA 51V255V (GTT→GTC) purHbifunctional AICAR transformylase/IMP cyclohydrolase
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (22/28);  minor base C (0/1);  total (22/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 = 1.00e+00

AGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAGCTTCTTAATGAAGTTCAGGTTCAGCGTGCGTGGGAAG  >  NC_000913/4206631‑4206900
                                                                                                                                         |                                                                                                                                    
aGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGAc                                                                                                                                    <  1:1513194/140‑1 (MQ=255)
aGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTAACGGCGATTTTGCTAAAAGAGACGGCTTTACCCTGGAc                                                                                                                                    >  1:650639/1‑140 (MQ=255)
     ggAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTggg                                                                                                                               >  1:842422/1‑140 (MQ=255)
      gAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGt                                                                                                                              <  1:1412550/140‑1 (MQ=255)
            ccGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGt                                                                                                                        >  2:1515921/1‑140 (MQ=255)
            ccGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGt                                                                                                                        >  2:1166546/1‑140 (MQ=255)
                tAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCa                                                                                                                    >  1:876669/1‑140 (MQ=255)
                   ccACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACg                                                                                                                 >  2:874024/1‑140 (MQ=255)
                      cgccgcAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGt                                                                                                                              >  1:892714/1‑124 (MQ=255)
                      cgccgcAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGt                                                                                                                              <  2:892714/124‑1 (MQ=255)
                      cgccgcAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAg                                                                                                              >  2:1463751/1‑140 (MQ=255)
                          gcAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGctt                                                                                                          <  2:842422/140‑1 (MQ=255)
                          gcAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGctt                                                                                                          <  2:128480/140‑1 (MQ=255)
                                  tGGCGTGCTTCACAATCGCACAGGCCGGCTCGGCGGACTGTTTCACGCTCTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTcaca                                                                                                  <  2:650639/140‑1 (MQ=255)
                                  tGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTcaca                                                                                                  <  2:390064/140‑1 (MQ=255)
                                  tGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTcaca                                                                                                  <  2:1918/140‑1 (MQ=255)
                                  tGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTcaca                                                                                                  <  2:1470121/140‑1 (MQ=255)
                                    gCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACAtt                                                                                                <  1:1463751/140‑1 (MQ=255)
                                     cGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACAtt                                                                                                <  1:988043/139‑1 (MQ=255)
                                     cGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACAtt                                                                                                >  2:988043/1‑139 (MQ=255)
                                         cTTCACAATCACACATGCCGGCTCGGCGAACCCTTGCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCtctt                                                                                           <  2:493642/140‑1 (MQ=255)
                                         cTTCACAATCACACACGCCGGCCCGTCGAGCTCTTTCTCGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTCCCCTGGACCTGGGTTGCGGTAGCAACGGGGGCTTCTTTCACATTCtctt                                                                                           <  1:532591/140‑1 (MQ=255)
                                           tCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGCTATCGGCGATGATGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCtcttct                                                                                         >  1:1301602/1‑140 (MQ=255)
                                               aaTCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATTAGAGAGGGGCTTACCCCGGACCTGGGTTACGGGAGCAACGGAGGTATCTTGCACATTCTCTTCtatat                                                                                     >  2:532591/1‑140 (MQ=255)
                                                 tCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACATGGGTTGCAGTAACAAATGAGGCGTGTTTCTAATTCTCTACTATATAg                                                                                   >  2:741015/1‑140 (MQ=255)
                                                       aTGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCa                                                                             >  1:1468588/1‑140 (MQ=255)
                                                         gCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGc                                                                           <  1:1166546/140‑1 (MQ=255)
                                                         gCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGc                                                                           <  1:874024/140‑1 (MQ=255)
                                                         gCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGc                                                                           <  1:317820/140‑1 (MQ=255)
                                                             gCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTCGCAACGGAAGCTCCTTTCACATTCTATTCTATATAGAAGGCAGCctgc                                                                       >  2:521258/1‑140 (MQ=255)
                                                                       ctctTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGt                                                             >  2:209746/1‑140 (MQ=255)
                                                                          tttCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTtct                                                          <  1:257767/140‑1 (MQ=255)
                                                                                gCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGt                                                    <  1:838864/140‑1 (MQ=255)
                                                                                  aCTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGCCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTaa                                                  <  2:1301602/140‑1 (MQ=255)
                                                                                            cgcATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCa                                        <  2:1099691/140‑1 (MQ=255)
                                                                                             gcATCGGTATCGGCGCTGTTGTAATCAGACAGGGCTTGACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAg                                       <  1:435045/140‑1 (MQ=255)
                                                                                             gcATCGGTATCGGCGATGTTGTTATCAGAGAGGGCTTGACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAg                                       <  1:521258/140‑1 (MQ=255)
                                                                                              cATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAGc                                      >  1:1433581/1‑140 (MQ=255)
                                                                                              cATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAGc                                      >  2:869610/1‑140 (MQ=255)
                                                                                              cATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAGc                                      >  1:930783/1‑140 (MQ=255)
                                                                                              cATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGAGGCAGTGCTCGCCGTAACGCATATCCAGc                                      >  1:627783/1‑140 (MQ=255)
                                                                                                  ggTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAGcttct                                  <  2:1468588/140‑1 (MQ=255)
                                                                                                      tCGGCGATGTTGTTATAAGAGAAGGCTTTACGCTGCACCGGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAGCTTCTTAAt                              <  2:627783/140‑1 (MQ=255)
                                                                                                        ggCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAGCTTCTTAATGa                            <  1:209746/140‑1 (MQ=255)
                                                                                                        ggCGATGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAGCTTCTTAATGa                            <  2:930783/140‑1 (MQ=255)
                                                                                                            aTGTTGTTATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTAAATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAGCTTCTTAATGAAgtt                        <  2:876669/140‑1 (MQ=255)
                                                                                                              gttgttATAAGAGAGGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAGCTTCTTAATGAAgttca                      >  1:1467536/1‑140 (MQ=255)
                                                                                                                       agagagGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAGCTTCTTAATGAAGTTCAGGTTCAgcg             <  1:869610/140‑1 (MQ=255)
                                                                                                                       agagagGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAGCTTCTTAATGAAGTTCAGGTTCAgcg             <  2:1433581/140‑1 (MQ=255)
                                                                                                                           agGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAGCTTCTTAATGAAGTTCAGGTTCAgcgtgcg         >  1:1458404/1‑140 (MQ=255)
                                                                                                                           agGGCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAGCTTCTTAATGAAGTTCAGGTTCAgcgtgcg         >  1:1154313/1‑140 (MQ=255)
                                                                                                                            gggCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAGCTTCTTAATGAAGTTCAGGTTCAgcgtgcgt        >  2:1161032/1‑140 (MQ=255)
                                                                                                                            gggCTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAGCTTCTTAATGAAGTTCAGGTTCAgcgtgcgt        >  1:862099/1‑140 (MQ=255)
                                                                                                                               cTTTACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAGCTTCTTAATGAAGTTCAGGTTCAGCGTGCGTggg     >  1:1238398/1‑140 (MQ=255)
                                                                                                                                  tACCCTGGACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAACTTCTTAATGAAGTTCAGGTTGAGCGTGCGTGGGAAg  >  1:608889/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                    
AGAATGGAATTGCCGATAGCCACGCCGCAAGGGTTGGCGTGCTTCACAATCACACATGCCGGCTCGGCGAACTCTTTCACGCACTCCAGCGCCGCATCGGTATCGGCGATGTTGTTATAAGAGAGGGCTTTACCCTGAACCTGGGTTGCGGTAGCAACGGAGGCTTCTTTCACATTCTCTTCTATATAGAAGGCAGCCTGCTGGTGGCTGTTCTCGCCGTAACGCATATCCAGCTTCTTAATGAAGTTCAGGTTCAGCGTGCGTGGGAAG  >  NC_000913/4206631‑4206900

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