Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 772,209 C→A T251N (ACC→AAC)  cydA → cytochrome d terminal oxidase, subunit I

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913772,2090CA98.2% 189.0 / ‑3.3 54T251N (ACC→AAC) cydAcytochrome d terminal oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base A (29/24);  total (29/25)
Fisher's exact test for biased strand distribution p-value = 4.63e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.04e-01

CCTTTGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAACCAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAGAGGAGACGAACAA  >  NC_000913/772113‑772305
                                                                                                |                                                                                                
ccTTTGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACaa                                                                                               <  1:518635/100‑1 (MQ=255)
    tGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTg                                                                                           >  2:1054884/1‑100 (MQ=255)
    tGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTg                                                                                           >  2:456170/1‑100 (MQ=255)
              gCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTg                                                                                 <  2:1303278/100‑1 (MQ=255)
              gCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTg                                                                                 >  2:727936/1‑100 (MQ=255)
               ccAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTg                                                                                 <  2:734604/99‑1 (MQ=255)
                      cGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAg                                                                         <  1:961469/100‑1 (MQ=255)
                        gTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTg                                                                       <  2:438232/100‑1 (MQ=255)
                         tATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTgg                                                                      >  2:1184475/1‑100 (MQ=255)
                         tATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTgg                                                                      <  1:489862/100‑1 (MQ=255)
                            ggCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGaa                                                                   >  2:75102/1‑100 (MQ=255)
                              ctgctgTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAAc                                                                 >  1:296391/1‑100 (MQ=255)
                               tgctgTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACg                                                                <  1:97836/100‑1 (MQ=255)
                                gctgTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGc                                                               <  1:328460/100‑1 (MQ=255)
                                   gTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAAc                                                            >  2:1289734/1‑100 (MQ=255)
                                    ttCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAAcc                                                           <  1:722893/100‑1 (MQ=255)
                                      ctgtctgtTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTg                                                         >  2:358345/1‑100 (MQ=255)
                                       tgtctgtTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTCTTGAAGCCGCGTGGGAAACGCAACCTgc                                                        >  1:354202/1‑100 (MQ=255)
                                       ngtctgtTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGNAAAACAAACTGGCTGCTATTGAAGCCGNGTGGGAAACGCAACCTgc                                                        >  2:554853/2‑100 (MQ=255)
                                        gtctgtTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTgcg                                                       >  1:1178841/1‑100 (MQ=255)
                                        gtctgtTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTgcg                                                       <  1:1151366/100‑1 (MQ=255)
                                          ctgtTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGcc                                                     >  1:756449/1‑100 (MQ=255)
                                            gtTATTGTTCTGGGTGATGAATCCGGCTACGAGATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCctg                                                   >  2:1174257/1‑100 (MQ=255)
                                              tATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCctgct                                                 >  2:875836/1‑100 (MQ=255)
                                              tATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCctgat                                                 >  1:675271/1‑98 (MQ=255)
                                                ttGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCtgctgc                                               <  2:192845/100‑1 (MQ=255)
                                                ctGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCtgctgc                                               <  1:1169863/99‑1 (MQ=255)
                                                  gTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAAAAAACTGGCTGCTATTGAAGCGGAGTGGGAAACTGAACATGGGCctgcgtcg                                              >  2:1169852/1‑95 (MQ=255)
                                                   ttCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCtt                                            >  1:704600/1‑100 (MQ=255)
                                                    tCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAACCAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCttt                                           <  1:1034776/100‑1 (MQ=255)
                                                    tCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCttt                                           <  2:635327/100‑1 (MQ=255)
                                                      tGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTAc                                         <  2:969823/100‑1 (MQ=255)
                                                      tGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTAc                                         <  2:971056/100‑1 (MQ=255)
                                                       gggTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTAct                                        <  1:935318/100‑1 (MQ=255)
                                                         gTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTActct                                      >  1:287035/1‑100 (MQ=255)
                                                         gTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTActct                                      >  2:1250860/1‑100 (MQ=255)
                                                            atgaATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGtt                                   >  2:228305/1‑100 (MQ=255)
                                                             tgaATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTAATGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTc                                  >  1:217350/1‑100 (MQ=255)
                                                               aaTCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCgg                                >  1:448204/1‑100 (MQ=255)
                                                                 tCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCa                              <  1:872796/100‑1 (MQ=255)
                                                                     gCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTcc                          <  1:727941/100‑1 (MQ=255)
                                                                       tACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTg                        <  2:261915/100‑1 (MQ=255)
                                                                       tACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTg                        >  2:1127706/1‑100 (MQ=255)
                                                                         cGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGAt                      >  1:521767/1‑100 (MQ=255)
                                                                         cGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGAt                      >  2:13277/1‑100 (MQ=255)
                                                                            aaTGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAg                   <  1:75103/100‑1 (MQ=255)
                                                                                    aCGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAgagga           <  1:492223/100‑1 (MQ=255)
                                                                                      gTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAGAGgaga         <  1:742848/100‑1 (MQ=255)
                                                                                       tGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAGAGGAGAc        >  2:19003/1‑100 (MQ=255)
                                                                                        gCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAGAGGAGAcn       <  1:1146929/100‑2 (MQ=255)
                                                                                         cAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGACACACAACCTGCGCCTGCTGTCTTTACTCTGTTCGGCATTCCTGATCAGGAAGAGGAGACGa      <  2:711574/100‑1 (MQ=255)
                                                                                         cAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAGAGGAGACGa      >  1:575907/1‑100 (MQ=255)
                                                                                           gAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAGAGGAGACGAAc    >  1:435539/1‑100 (MQ=255)
                                                                                           gAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAGAGGAGACGAAc    >  1:545750/1‑100 (MQ=255)
                                                                                             aaaaCAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAGAGGAGACGAACaa  >  2:81241/1‑100 (MQ=255)
                                                                                             aaaaCAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAGAGGAGACGAACaa  >  1:1048863/1‑100 (MQ=255)
                                                                                                |                                                                                                
CCTTTGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAACCAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAGAGGAGACGAACAA  >  NC_000913/772113‑772305

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