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.8% 240.0 / 0.8 71I31F (ATC→TTC) dgtdeoxyguanosine triphosphate triphosphohydrolase
Reads supporting (aligned to +/- strand):  ref base A (1/2);  new base T (32/36);  total (33/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01

GGACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCC  >  NC_000913/179233‑179419
                                                                                              |                                                                                            
ggACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGttcttc                                                                                         >  2:476147/1‑100 (MQ=255)
 gACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCg                                                                                        <  1:1079794/100‑1 (MQ=255)
  aCATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCga                                                                                       >  2:584363/1‑100 (MQ=255)
   cATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCgag                                                                                      <  2:709727/100‑1 (MQ=255)
     tGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGagag                                                                                    <  1:1230300/100‑1 (MQ=255)
      ggCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGc                                                                                   <  1:951044/100‑1 (MQ=255)
      ggCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGc                                                                                   <  1:1076388/100‑1 (MQ=255)
      ggCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACAGCAGGGCGGTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGc                                                                                   >  2:255521/1‑100 (MQ=255)
        cacaGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGa                                                                                 >  2:559770/1‑100 (MQ=255)
        cacaGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGa                                                                                 >  2:322219/1‑100 (MQ=255)
            gattgattTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATcgc                                                                             <  1:322219/100‑1 (MQ=255)
              ttgattTCCAAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCgg                                                                           <  1:358926/100‑1 (MQ=255)
                gattTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGc                                                                         <  1:435311/100‑1 (MQ=255)
                 attTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAGCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCg                                                                        <  2:1334693/100‑1 (MQ=255)
                 attTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCg                                                                        <  1:1013188/100‑1 (MQ=255)
                  tttCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGt                                                                       >  1:15329/1‑100 (MQ=255)
                    tCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTat                                                                     <  1:1001427/100‑1 (MQ=255)
                    tCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTat                                                                     >  2:1482833/1‑100 (MQ=255)
                    tCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTat                                                                     >  2:1169323/1‑100 (MQ=255)
                      cGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatca                                                                   <  2:28455/100‑1 (MQ=255)
                      cGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatca                                                                   <  2:981407/100‑1 (MQ=255)
                       gAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTatcat                                                                  <  1:556021/100‑1 (MQ=255)
                        aaaaaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTatcat                                                                  <  1:781862/99‑1 (MQ=255)
                         aaaaaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTAtcatca                                                                <  1:283029/100‑1 (MQ=255)
                            aaaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAAct                                                             >  2:46865/1‑100 (MQ=255)
                             aaaTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAActc                                                            >  1:1426380/1‑100 (MQ=255)
                               aTAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACtctc                                                          >  1:1570138/1‑100 (MQ=255)
                                tAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTcc                                                         <  1:584365/100‑1 (MQ=255)
                                 aaaCTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTcc                                                         <  2:1496391/99‑1 (MQ=255)
                                 aaaCTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCg                                                        <  2:931907/100‑1 (MQ=255)
                                  aaCTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCgg                                                       <  2:1263163/100‑1 (MQ=255)
                                  aaCTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCgg                                                       <  1:881820/100‑1 (MQ=255)
                                   aCTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGc                                                      <  1:1311164/100‑1 (MQ=255)
                                   aCTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGc                                                      <  1:35686/100‑1 (MQ=255)
                                         aTCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcg                                                >  1:1179547/1‑100 (MQ=255)
                                          tcgtcgTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgt                                               >  2:634396/1‑100 (MQ=255)
                                          tcgtcgTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgt                                               >  2:1497770/1‑100 (MQ=255)
                                          tcgtcgTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCCTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgt                                               >  1:309951/1‑100 (MQ=255)
                                            gtcgtTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgtcg                                             <  2:117991/100‑1 (MQ=255)
                                            gtcgtTACCGTTCACCGCAGGGCGTGAAAGCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATtcgtcg                                             <  2:1298669/100‑1 (MQ=255)
                                              cgtTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCgtcgtc                                           >  1:567317/1‑100 (MQ=255)
                                               gtTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCt                                          >  1:209725/1‑100 (MQ=255)
                                                ttACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTg                                         <  1:1497779/100‑1 (MQ=255)
                                                    cGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaac                                     >  1:1236722/1‑100 (MQ=255)
                                                     gTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGcaaca                                    >  1:840462/1‑100 (MQ=255)
                                                        cACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACaaaa                                 <  2:15329/100‑1 (MQ=255)
                                                         aCCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGTCGTATCGTCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAg                                >  1:314333/1‑100 (MQ=255)
                                                            gaaGGGCGTTAAAACCGAACATAAGATCCTGCGGTTCTGCGAGAGCGAGCGCGGGCGTATCATCAACGCTCCGGCAATTCGTCGTCGGCAACAAAAGAcc                             >  2:920546/3‑100 (MQ=255)
                                                              aGGGCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCa                           <  1:1313790/100‑1 (MQ=255)
                                                                ggCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAgg                         <  1:1482842/100‑1 (MQ=255)
                                                                 gCGTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGt                        >  2:1551015/1‑100 (MQ=255)
                                                                   gTTAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGttt                      >  1:804377/1‑100 (MQ=255)
                                                                     tAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACATAAGACCCAAGttttt                    >  2:819649/1‑100 (MQ=255)
                                                                     tAAAACCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGttttt                    >  2:902807/1‑100 (MQ=255)
                                                                      aaaaCCGACCATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTCTCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTc                   <  1:294153/100‑1 (MQ=255)
                                                                      aaaaCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTc                   >  2:403400/1‑100 (MQ=255)
                                                                      aaaaCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTc                   <  1:1288270/100‑1 (MQ=255)
                                                                       aaaCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTcc                  <  1:422501/100‑1 (MQ=255)
                                                                       aaaCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTcc                  <  1:295155/100‑1 (MQ=255)
                                                                       aaaCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTcc                  >  2:491431/1‑100 (MQ=255)
                                                                       aaaCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTcc                  <  2:567316/100‑1 (MQ=255)
                                                                         aCCGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCAc                <  1:379494/100‑1 (MQ=255)
                                                                          ccGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACt               >  2:1369728/1‑100 (MQ=255)
                                                                           cGAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTg              >  1:1343264/1‑100 (MQ=255)
                                                                           cGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTg              <  1:287640/100‑1 (MQ=255)
                                                                            gAACATGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTgg             <  1:24400/100‑1 (MQ=255)
                                                                                 tGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAgcgc        >  2:1402321/1‑100 (MQ=255)
                                                                                 tGAGATCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTCTTCCACTGGAgcgc        >  2:1079121/1‑100 (MQ=255)
                                                                                   agaTCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCaa      >  1:1481911/1‑100 (MQ=255)
                                                                                    gaTCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACACAAGACCCAGGTTTTTCCACTGGAGCGCAAt     <  2:667600/100‑1 (MQ=255)
                                                                                    gaTCCTGCGGTTCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAAt     >  1:1016202/1‑100 (MQ=255)
                                                                                       acTGCGGTTGTTCGAGAGAGATCGCGGGCGTATCATCAGCTGTCCGGCGATGAGTAGTCGGCAAAAAACGACAGAGGTGGTGTCACTGGAGCGCAAttac  >  2:379912/2‑97 (MQ=255)
                                                                                              |                                                                                            
GGACATGGCACAGATTGATTTCCGAAAAAAAATAAACTGGCATCGTCGTTACCGTTCACCGCAGGGCGTTAAAACCGAACATGAGATCCTGCGGATCTTCGAGAGCGATCGCGGGCGTATCATCAACTCTCCGGCAATTCGTCGTCTGCAACAAAAGACCCAGGTTTTTCCACTGGAGCGCAATGCC  >  NC_000913/179233‑179419

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