Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 179,327 A→T I31F (ATC→TTC)  dgt → deoxyguanosine triphosphate triphosphohydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913179,3270AT95.4% 219.9 / 1.4 65I31F (ATC→TTC) dgtdeoxyguanosine triphosphate triphosphohydrolase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base T (34/28);  total (35/30)
Fisher's exact test for biased strand distribution p-value = 5.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01

AGGACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCG  >  NC_000913/179232‑179423
                                                                                               |                                                                                                
agGACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCtt                                                                                              <  2:605606/100‑1 (MQ=255)
 ggACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGttcttc                                                                                             >  1:96637/1‑100 (MQ=255)
  gACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCg                                                                                            >  2:1557905/1‑100 (MQ=255)
    cATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCgag                                                                                          <  2:1669891/100‑1 (MQ=255)
    cATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCgag                                                                                          >  1:1600134/1‑100 (MQ=255)
        gCACAGATTGATTTCCGAAAAACAATAAACTGGCCTCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGCGATCCTGCGGttcttccaca                                                                                         >  1:1693087/1‑93 (MQ=255)
         cacaGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGa                                                                                     >  1:927824/1‑100 (MQ=255)
           caGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATc                                                                                   >  1:896512/1‑100 (MQ=255)
            aGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATcg                                                                                  >  1:614703/1‑100 (MQ=255)
                  attTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCg                                                                            <  1:1884647/100‑1 (MQ=255)
                      ccGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatc                                                                        >  1:526569/1‑100 (MQ=255)
                      ccGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatc                                                                        >  1:175697/1‑100 (MQ=255)
                      ccGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatc                                                                        >  2:1150807/1‑100 (MQ=255)
                      ccGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatc                                                                        >  2:295691/1‑100 (MQ=255)
                      ccGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatc                                                                        >  1:1031874/1‑100 (MQ=255)
                       cGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatca                                                                       >  1:406402/1‑100 (MQ=255)
                         aaaaaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatcatc                                                                     >  2:271463/1‑100 (MQ=255)
                         aaaaaaaaTAAACTGGCATCGACGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatcatc                                                                     >  1:1086371/1‑100 (MQ=255)
                          aaaaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTAtcatca                                                                    <  2:1127981/100‑1 (MQ=255)
                          aaaaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTAtcatca                                                                    <  2:1139154/100‑1 (MQ=255)
                          aaaaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTAtcatca                                                                    >  2:974821/1‑100 (MQ=255)
                            aaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAAc                                                                  <  2:284730/100‑1 (MQ=255)
                             aaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAAct                                                                 >  2:935528/1‑100 (MQ=255)
                             aaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAAct                                                                 >  2:750504/1‑100 (MQ=255)
                             aaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAAct                                                                 >  1:1179459/1‑100 (MQ=255)
                             aaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAAc                                                                  <  1:1381056/99‑1 (MQ=255)
                              aaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAActc                                                                >  2:166635/1‑100 (MQ=255)
                                 tAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTcc                                                             <  2:1158611/100‑1 (MQ=255)
                                   aaCTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCgg                                                           >  1:1320202/1‑100 (MQ=255)
                                      tGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCaa                                                        >  2:1273801/1‑100 (MQ=255)
                                           tcgtcgTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgt                                                   <  2:1856607/100‑1 (MQ=255)
                                             gtcgtTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgtcg                                                 <  2:121603/100‑1 (MQ=255)
                                             gtcgtTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgtcg                                                 >  1:948647/1‑100 (MQ=255)
                                                gtTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCt                                              >  2:493044/1‑100 (MQ=255)
                                                 ttACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTg                                             <  2:948648/100‑1 (MQ=255)
                                                  tACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGc                                            >  2:41773/1‑100 (MQ=255)
                                                    ccGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaa                                          <  2:450083/100‑1 (MQ=255)
                                                    ccGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaa                                          >  1:884750/1‑100 (MQ=255)
                                                     cGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaac                                         <  1:1350407/100‑1 (MQ=255)
                                                      gTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaaca                                        >  2:1625654/1‑100 (MQ=255)
                                                      gTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaaca                                        >  2:153145/1‑100 (MQ=255)
                                                           ccGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGa                                   >  2:1450281/1‑100 (MQ=255)
                                                            cGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGAc                                  <  1:1150805/100‑1 (MQ=255)
                                                                gggCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAg                              <  1:1598549/100‑1 (MQ=255)
                                                                   cGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGtt                           <  2:538973/100‑1 (MQ=255)
                                                                      tAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGttttt                        >  2:310109/1‑100 (MQ=255)
                                                                       aaaaCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAGGACCCAGGTTTTTc                       <  2:884751/100‑1 (MQ=255)
                                                                       aaaaCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTc                       <  1:935527/100‑1 (MQ=255)
                                                                        aaaCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTcc                      <  2:1356138/100‑1 (MQ=255)
                                                                        aaaCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTcc                      <  2:105738/100‑1 (MQ=255)
                                                                           ccGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACt                   <  1:1571342/100‑1 (MQ=255)
                                                                           ccGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACt                   >  1:751542/1‑100 (MQ=255)
                                                                            cGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTg                  <  2:654959/100‑1 (MQ=255)
                                                                            cGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTg                  <  2:191799/100‑1 (MQ=255)
                                                                             gAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTgg                 >  1:1550704/1‑100 (MQ=255)
                                                                             gAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTgg                 <  2:927825/100‑1 (MQ=255)
                                                                               aCATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAg               <  1:551308/100‑1 (MQ=255)
                                                                               aCATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAg               <  1:1344291/100‑1 (MQ=255)
                                                                                  tGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgcgc            <  2:1531161/100‑1 (MQ=255)
                                                                                   gagaTCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCa           <  2:1332507/100‑1 (MQ=255)
                                                                                   gagaTCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCa           >  1:1690363/1‑100 (MQ=255)
                                                                                    agaTCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCaa          >  1:1825409/1‑100 (MQ=255)
                                                                                      aTCCTGCGGTTCTCCGAGAGCGATCGCGGGCGTATCATAAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATg        >  1:1200372/1‑100 (MQ=255)
                                                                                        ccTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATgcc      <  1:91131/100‑1 (MQ=255)
                                                                                          tGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATgccgc    <  1:1599818/100‑1 (MQ=255)
                                                                                           gCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATgccgcc   >  2:403604/1‑100 (MQ=255)
                                                                                            cGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGccgccg  >  2:1833017/1‑100 (MQ=255)
                                                                                               |                                                                                                
AGGACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCCGCCG  >  NC_000913/179232‑179423

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