Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,961,588 C→G R86T (AGG→ACG)  gppA ← guanosine pentaphosphatase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,961,5880CG99.5% 494.3 / NA 184R86T (AGG→ACG) gppAguanosine pentaphosphatase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base G (100/83);  minor base T (0/1);  total (100/84)
Fisher's exact test for biased strand distribution p-value = 4.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.97e-01

Only 100 of 184 total aligned reads displayed.

CACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACCTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATT  >  NC_000913/3961515‑3961660
                                                                         |                                                                        
cACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTc                                                                        >  2:1401081/1‑76 (MQ=255)
 cccTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCg                                                                       <  1:529416/75‑1 (MQ=255)
 aCCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCg                                                                       >  1:1686039/1‑76 (MQ=255)
    tGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCg                                                                    >  2:1161389/1‑76 (MQ=255)
     ggaCCGGACAGCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGt                                                                   <  2:756577/74‑1 (MQ=255)
      tACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTa                                                                  <  1:4434648/76‑1 (MQ=255)
      tACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTa                                                                  >  2:4812236/1‑76 (MQ=255)
      gaCCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTa                                                                  <  2:3684361/75‑1 (MQ=255)
       cccGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCGCCCGCATTGACGGCAAGGCGTAACGTCGCCGTAg                                                                 <  2:2278607/75‑1 (MQ=255)
       aCCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAg                                                                 <  2:3353332/76‑1 (MQ=255)
       aCCCGACAACCGAGACTTTCCTGCGCTTTGGCAATAAAACCACCCGCATTGACGGCAAGGCGTAACGTCGCCGCAg                                                                 >  2:3062828/1‑76 (MQ=255)
        ccGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGc                                                                >  1:3013805/1‑76 (MQ=255)
         cGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCg                                                               >  1:4577118/1‑76 (MQ=255)
          ggACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGa                                                              >  1:4030050/1‑76 (MQ=255)
           gACAGCCGGGGATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGAc                                                             <  2:668847/76‑1 (MQ=255)
             cAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACaa                                                           <  2:4299744/76‑1 (MQ=255)
             cAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACaa                                                           >  2:5366144/1‑76 (MQ=255)
             cAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACaa                                                           <  1:4158626/76‑1 (MQ=255)
              aaCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAc                                                          <  2:315644/76‑1 (MQ=255)
              aaCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAc                                                          <  2:1588022/76‑1 (MQ=255)
               aCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcg                                                         <  2:1904303/76‑1 (MQ=255)
               aCCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcg                                                         >  1:324977/1‑76 (MQ=255)
                ccGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcgc                                                        >  2:693992/1‑76 (MQ=255)
                ccGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcgc                                                        <  1:1008927/76‑1 (MQ=255)
                ccGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcgc                                                        <  2:796033/76‑1 (MQ=255)
                 cGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcgcg                                                       >  1:576748/1‑76 (MQ=255)
                 cGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAcgcg                                                       >  1:874914/1‑76 (MQ=255)
                  gagaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGa                                                      >  2:5526728/1‑76 (MQ=255)
                  gagaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGa                                                      >  1:4993970/1‑76 (MQ=255)
                   agaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGaa                                                     <  2:1624876/76‑1 (MQ=255)
                    gaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAAt                                                    >  2:1733630/1‑76 (MQ=255)
                    gaATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAAt                                                    >  1:2557901/1‑76 (MQ=255)
                    aaaTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAAt                                                    <  1:447192/75‑1 (MQ=255)
                     gaTTTCCTGCGCTTTGGCAATAAAAGCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAAtt                                                   <  2:5577472/75‑1 (MQ=255)
                     aaTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAAtt                                                   <  2:4422733/76‑1 (MQ=255)
                      aTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAAttt                                                  <  1:1034531/76‑1 (MQ=255)
                      aTTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAAttt                                                  <  1:4378678/76‑1 (MQ=255)
                       tttCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTg                                                 >  1:2412242/1‑76 (MQ=255)
                       tttCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTg                                                 >  1:1811432/1‑76 (MQ=255)
                       tttCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTg                                                 <  2:1232316/76‑1 (MQ=255)
                        ttCGTGCGCTGTGGCAATAAACTCGCCCGCATTGACGGCAAGGCGTAACGGCGCCGTCGCGCCAGCGCGAATTTGc                                                <  2:2961234/76‑1 (MQ=255)
                        ttCCTGCGCTGTGTCAACAAAATCACCTGCCATGACGCCAAGCACAGACGTCGCCGGAGCGACAACGCGGATTTGc                                                >  2:3469003/1‑76 (MQ=255)
                           cTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCgag                                             <  1:5037609/76‑1 (MQ=255)
                           cTGCGCGTTGGCAATAAAAACACCCGCATGGACGGCAAGGCGTAACGGCGCCGTAGCGACAACGCGAATTTGCgag                                             <  2:243935/76‑1 (MQ=255)
                              cgcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCgaggga                                          >  1:5554392/1‑76 (MQ=255)
                               gcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCgagggag                                         >  2:3562167/1‑76 (MQ=255)
                               gcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCgagggag                                         >  1:423441/1‑76 (MQ=255)
                               gcTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCgagggag                                         <  1:693992/76‑1 (MQ=255)
                                cTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACACCGCGAATTTGCgcgggagg                                        >  2:5770808/1‑76 (MQ=255)
                                cTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCgagggcgg                                        >  1:5385367/1‑76 (MQ=255)
                                 tttGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGggagga                                       <  2:4719264/76‑1 (MQ=255)
                                  ttGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGaa                                      >  2:2688536/1‑76 (MQ=255)
                                  ttGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACGGCGCCGTAGCGACATAGAGCATTTGCGCGGGAGGaa                                      >  1:3866136/1‑76 (MQ=255)
                                    ggCAATAAAATCCCCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAata                                    <  2:5041597/76‑1 (MQ=255)
                                      cAATAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGACGTAGCGACAACGCGAATTTGTGAGGGAGGGATATc                                  >  1:5505830/1‑76 (MQ=255)
                                       aaTAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCt                                 <  2:3243128/76‑1 (MQ=255)
                                       aaTAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCt                                 >  1:3725940/1‑76 (MQ=255)
                                        aTAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCtt                                <  1:2684565/76‑1 (MQ=255)
                                        aTAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCtt                                <  2:2063351/76‑1 (MQ=255)
                                        aTAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCtt                                >  2:4946797/1‑76 (MQ=255)
                                        aTAAAATCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCtt                                <  1:1401081/76‑1 (MQ=255)
                                          aaaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGc                              <  2:2948258/76‑1 (MQ=255)
                                           aaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCa                             >  2:5731437/1‑76 (MQ=255)
                                           aaaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCa                             <  1:4334442/76‑1 (MQ=255)
                                            aaTCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAg                            >  2:2310640/1‑76 (MQ=255)
                                              tcccccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAAAGCGAATTTGCGAGGGAGGAATATCTTGCAGAc                          >  1:354487/4‑76 (MQ=255)
                                              tCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGAc                          <  1:4528216/76‑1 (MQ=255)
                                              tCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGAc                          <  1:4710224/76‑1 (MQ=255)
                                              tCACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGAc                          >  2:5226203/1‑76 (MQ=255)
                                               cACCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACg                         <  1:3992910/76‑1 (MQ=255)
                                                aCCCGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTCGCGAGGGAGCAATATCTTCCACACGt                        >  2:52831/1‑76 (MQ=255)
                                                  ccGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTc                      >  1:5154899/1‑76 (MQ=255)
                                                   cGCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCa                     <  2:2864906/76‑1 (MQ=255)
                                                    gCATTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAg                    <  1:52606/76‑1 (MQ=255)
                                                      aTTGGCGGCCAGGGGTAACTTCGCCGTAGCGACACCGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCa                  <  2:3866135/76‑1 (MQ=255)
                                                      aTTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGcg                  >  2:1828988/1‑75 (MQ=255)
                                                      aTTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCa                  >  2:2339581/1‑76 (MQ=255)
                                                      aTTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCa                  <  1:3405849/76‑1 (MQ=255)
                                                      aTTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCa                  >  1:1862225/1‑76 (MQ=255)
                                                      aTTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCa                  >  1:2863354/1‑76 (MQ=255)
                                                      aTTGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCa                  >  1:3801587/1‑76 (MQ=255)
                                                       ttGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCCGACGTTCCGCaa                 >  2:2606339/1‑76 (MQ=255)
                                                       ttGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTAAGCaa                 >  2:5416648/1‑76 (MQ=255)
                                                        tGACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCaac                >  2:1248809/1‑75 (MQ=255)
                                                         gACGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAAc               >  2:3504788/1‑76 (MQ=255)
                                                          aCGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACa              >  1:376333/1‑76 (MQ=255)
                                                          aCGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACa              <  2:4441124/76‑1 (MQ=255)
                                                           cGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAg             <  1:4284068/76‑1 (MQ=255)
                                                           cGGCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAg             <  1:5731439/76‑1 (MQ=255)
                                                            ggCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAgg            >  1:2217834/1‑76 (MQ=255)
                                                            ggCAAGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAgg            <  1:5556242/76‑1 (MQ=255)
                                                               aaGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgcgc         >  1:197574/1‑76 (MQ=255)
                                                               aaGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgcgc         >  1:4653682/1‑76 (MQ=255)
                                                               aaGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGgcgc         >  2:5796389/1‑76 (MQ=255)
                                                               aaGGCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGCTCAGCAAACAGgcgc         >  1:1259109/1‑76 (MQ=255)
                                                                agtcGTCACGTCGCCGCAGCGACACCGCGAATTTGCGAGGGAGGAATATCTTGCAGGCGTTCAGCAAACAGGCGCa        <  1:1589299/73‑1 (MQ=255)
                                                                  gCGTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCaaa      <  1:618480/76‑1 (MQ=255)
                                                                    gTAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACa    >  2:198030/1‑76 (MQ=255)
                                                                     tAACGTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACAt   <  1:276907/76‑1 (MQ=255)
                                                                      aaCGTCGCCGTAGCGACAACGCGAATTTGCGTGGGAGGAATATCTTGCAGACGTTCAGCAAACCGGCGCAAACAtt  <  2:3919252/76‑1 (MQ=255)
                                                                         |                                                                        
CACCTGTACCGGACAACCGAGAATTTCCTGCGCTTTGGCAATAAAATCACCCGCATTGACGGCAAGGCGTAACCTCGCCGTAGCGACAACGCGAATTTGCGAGGGAGGAATATCTTGCAGACGTTCAGCAAACAGGCGCAAACATT  >  NC_000913/3961515‑3961660

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