Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,235,656 T→C I131T (ATC→ACC)  preA → dihydropyrimidine dehydrogenase, NADH‑dependent, subunit C

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,235,6560TC91.2% 174.3 / 6.1 57I131T (ATC→ACC) preAdihydropyrimidine dehydrogenase, NADH‑dependent, subunit C
Reads supporting (aligned to +/- strand):  ref base T (2/3);  new base C (33/19);  total (35/22)
Fisher's exact test for biased strand distribution p-value = 3.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.19e-01

ACCCGGACAAAGTATTGATCGCTTCGATCATGGGGGAAAATGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGATCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGCAAAGCCCGGAGCTGGTAGAAAAATATTGTCGGGCAGT  >  NC_000913/2235560‑2235752
                                                                                                |                                                                                                
aCCCGGACAAAGTATTGATCGCTTCGATCATGGGGGAAAATGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGa                                                                                               >  2:1430983/1‑100 (MQ=255)
aCCCGGACAAAGTATTGATCGCTTCGATCATGGGGGAAAATGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGa                                                                                               >  2:931603/1‑100 (MQ=255)
 cccGGACAAAGTATTGATCGCTTCGATCATGGGGGAAAATGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAg                                                                                              >  1:1076970/1‑100 (MQ=255)
 cccGGACAAAGTATTGATCGCTTCGATCATGGGGGAAAATGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAg                                                                                              >  1:178288/1‑100 (MQ=255)
         aaGTATTGATCGCTTCGATCATGGGGGAAAATGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACtt                                                                                      >  1:508610/1‑100 (MQ=255)
         aaGTATTGATCGCTTCGATCATGGGGGAAAATGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACtt                                                                                      >  1:314628/1‑100 (MQ=255)
           gTATTGATCGCTTCGATCATGGGGGAAAATGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTtct                                                                                    <  2:527851/100‑1 (MQ=255)
             aTTGATCGCTTCGATCATGGGGGAAAATGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTcc                                                                                  <  1:1215015/100‑1 (MQ=255)
             aTTGATCGCTTCGATCATGGGGGAAAATGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTcc                                                                                  <  1:1213721/100‑1 (MQ=255)
              ttGATCGCTTCGATCATGGGGGAAAATGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCt                                                                                 <  1:409646/100‑1 (MQ=255)
                gATCGCTTCGATCATGGGGGAAAATGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGt                                                                               >  2:1687159/1‑100 (MQ=255)
                   cGCTTCGATCATGGGGGAAAATGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTcc                                                                             <  1:755284/99‑1 (MQ=255)
                   cGCTTCGATCATGGGGGAAAATGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCg                                                                            >  1:1582774/1‑100 (MQ=255)
                   cGCTTCGATCATGGGGGAAAATGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCg                                                                            >  2:1647843/1‑100 (MQ=255)
                         gATCATGGGGGAAAATGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATg                                                                      >  1:1744327/1‑100 (MQ=255)
                         gATCATGGGGGAAAATGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATg                                                                      >  1:552393/1‑100 (MQ=255)
                          aTCATGGGGGAAAATGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGATCGAGTGTAACTTCTCCTGTCCGCAAATGa                                                                     <  2:1832086/100‑1 (MQ=255)
                           tCATGGGGGAAAATGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGAc                                                                    <  2:278086/100‑1 (MQ=255)
                           tCATGGGGGAAAATGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGAc                                                                    >  1:782027/1‑100 (MQ=255)
                               gggggAAAATGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTtct                                                                >  2:1194682/1‑100 (MQ=255)
                                ggggAAAATGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGATCGAGTGTAACCTCTCCTGTCCGCAAATGACTtctc                                                               <  1:1257947/100‑1 (MQ=255)
                                 gggAAAATGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCa                                                              >  2:1297363/1‑100 (MQ=255)
                                  ggAAAATGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCAt                                                             <  1:726753/100‑1 (MQ=255)
                                    aaaaTGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGATCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGc                                                           >  1:133097/1‑100 (MQ=255)
                                     aaaTGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCTAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCg                                                          <  1:1194680/100‑1 (MQ=255)
                                     aaaTGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCg                                                          <  2:1071265/100‑1 (MQ=255)
                                       aTGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGAt                                                        >  1:846233/1‑100 (MQ=255)
                                        tGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATg                                                       >  2:908519/1‑100 (MQ=255)
                                         gAGCAGCCATGGGAGGAGCTGGCGCGCCTGGTGCAAGCAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATgg                                                      >  1:920938/1‑100 (MQ=255)
                                           gcagcaATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGt                                                    >  1:11463/1‑100 (MQ=255)
                                              gcaATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGATCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGc                                                 >  2:627922/1‑100 (MQ=255)
                                                     gaggagCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCg                                          >  2:1654754/1‑100 (MQ=255)
                                                        gagCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGc                                       <  2:1065160/100‑1 (MQ=255)
                                                        gagCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGc                                       <  1:1713671/100‑1 (MQ=255)
                                                           cTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGCaaa                                    >  2:92794/1‑100 (MQ=255)
                                                           cTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGCaaa                                    >  1:1430448/1‑100 (MQ=255)
                                                            tGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGCAAAg                                   >  2:640622/1‑100 (MQ=255)
                                                                 cgcCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGCAAAGCCCgg                              <  1:1687149/100‑1 (MQ=255)
                                                                  gcCTGGTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGCAAAGCCCGGa                             >  2:311680/1‑100 (MQ=255)
                                                                      ggTGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGCAAAGCCCGGAGCTg                         >  2:1000349/1‑100 (MQ=255)
                                                                        tGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGCAAAGCCCGGAGCTgg                        >  2:1482754/1‑99 (MQ=255)
                                                                        tGCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGCAAAGCCCGGAGCTGGt                       <  1:497548/100‑1 (MQ=255)
                                                                         gCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGCAAAGCCCGGAGCTGGTa                      <  2:178289/100‑1 (MQ=255)
                                                                         gCAAGAAGCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGCAAAGCCCGGAGCTGGTa                      >  1:999968/1‑100 (MQ=255)
                                                                            agaagCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGCAAAGCCCGGAGCTGGTAGaa                   <  2:551514/100‑1 (MQ=255)
                                                                              aagCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGCAAAGCCCGGAGCTGGTAGaaaa                 >  1:1455069/1‑100 (MQ=255)
                                                                               agCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGCAAAGCCCGGAGCTGGTAGaaaaa                <  2:134787/100‑1 (MQ=255)
                                                                                gCTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGCAAAGCCCGGAGCTGGTAGAAAAat               >  1:1157679/1‑100 (MQ=255)
                                                                                 cTGGCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGCAAAGCCCGGAGCTGGTAGAAAAata              <  2:655651/100‑1 (MQ=255)
                                                                                  tGGCGCGGATATGATCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGCAAAGCCCGGAGCTGGTAGAAAAatat             <  1:1338272/100‑1 (MQ=255)
                                                                                   ggCGCGGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGCAAAGCCCGGAGCTGGTAGAAAAATAtt            <  2:94450/100‑1 (MQ=255)
                                                                                       cgGATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGCAAAGCCCGGAGCTGGTAGAAAAATATTGTCg        >  2:343991/1‑100 (MQ=255)
                                                                                         gATATGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGCAAAGCCCGGAGCTGGTAGAAAAATATTGTCgg       >  2:84315/1‑99 (MQ=255)
                                                                                           tatGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGCAAAGCCCGGAGCTGGTAGAAAAATATTGTCGGGCa    >  2:274482/1‑100 (MQ=255)
                                                                                           tatGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGCAAAGCCCGGAGCTGGTAGAAAAATATTGTCGGGCa    >  2:175140/1‑100 (MQ=255)
                                                                                            atGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGCAAAGCCCGGAGCTGGTAGAAAAATATTGTCGGGCAg   <  2:359997/100‑1 (MQ=255)
                                                                                             tGACCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGCAAAGCCCGGAGCTGGTAGAAAAATATTGTCGGGCAGt  >  2:1040273/1‑100 (MQ=255)
                                                                                                |                                                                                                
ACCCGGACAAAGTATTGATCGCTTCGATCATGGGGGAAAATGAGCAGCAATGGGAGGAGCTGGCGCGCCTGGTGCAAGAAGCTGGCGCGGATATGATCGAGTGTAACTTCTCCTGTCCGCAAATGACTTCTCATGCGATGGGTAGCGATGTCGGGCAAAGCCCGGAGCTGGTAGAAAAATATTGTCGGGCAGT  >  NC_000913/2235560‑2235752

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