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,2090CA93.4% 247.8 / 9.0 76T251N (ACC→AAC) cydAcytochrome d terminal oxidase, subunit I
Reads supporting (aligned to +/- strand):  ref base C (2/3);  new base A (33/38);  total (35/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.31e-01

CCTTTGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAACCAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAGAGGAGACGAACA  >  NC_000913/772113‑772304
                                                                                                |                                                                                               
ccTTTGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACaa                                                                                              >  1:1238269/1‑100 (MQ=255)
ccTTTGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACaa                                                                                              <  2:731155/100‑1 (MQ=255)
ccTGTGCTAGCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGGCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACa                                                                                               <  2:1093571/99‑1 (MQ=255)
 cTTTGCTATAGCTGCAAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGCGAAGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACaaa                                                                                             >  2:1438452/1‑100 (MQ=255)
   ttGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGTAATGGGCGACGTGCAGAAAAACAAACt                                                                                           >  2:433257/1‑100 (MQ=255)
   ttGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACt                                                                                           <  1:1712180/100‑1 (MQ=255)
    tGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTg                                                                                          >  2:466224/1‑100 (MQ=255)
     gCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGATTTGGGCGACGTGCAGAAAAACAAACTgg                                                                                         >  1:1709993/1‑100 (MQ=255)
        aTCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGgctg                                                                                      >  2:780919/1‑100 (MQ=255)
        aTCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGgctg                                                                                      >  1:1560335/1‑100 (MQ=255)
           gCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTa                                                                                   >  1:628740/1‑100 (MQ=255)
             tGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAACCAAACTGGCTGCTAtt                                                                                 >  1:562853/1‑100 (MQ=255)
              gCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTg                                                                                >  1:925423/1‑100 (MQ=255)
              gCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTg                                                                                <  2:1238263/100‑1 (MQ=255)
              gCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTg                                                                                <  2:142295/100‑1 (MQ=255)
               ccAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAACCAAACTGGCTGCTATTGa                                                                               <  2:562851/100‑1 (MQ=255)
               ccAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGa                                                                               <  1:1402195/100‑1 (MQ=255)
                   cTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGcc                                                                           <  2:1780731/100‑1 (MQ=255)
                   cTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGcc                                                                           <  1:1880667/100‑1 (MQ=255)
                        gTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTg                                                                      >  1:940803/1‑100 (MQ=255)
                        gTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTg                                                                      <  1:984897/100‑1 (MQ=255)
                         tATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTgg                                                                     >  1:267404/1‑100 (MQ=255)
                           tGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAAGCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGa                                                                   <  2:1244885/100‑1 (MQ=255)
                             gctgctGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGaaa                                                                 >  1:721981/1‑100 (MQ=255)
                              ctgctgTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGaatc                                                                >  2:1056128/1‑98 (MQ=255)
                              ctgctgTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAAc                                                                <  1:357972/100‑1 (MQ=255)
                              ctgctgTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAAc                                                                <  1:922635/100‑1 (MQ=255)
                               tgctgTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACg                                                               <  1:1779087/100‑1 (MQ=255)
                                     tctgtctgTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCt                                                         <  2:706513/100‑1 (MQ=255)
                                      ctgtctgtTATTGTTCTGGGTGAGGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTg                                                        <  2:628738/100‑1 (MQ=255)
                                        gtctgtTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTgcg                                                      <  2:1650646/100‑1 (MQ=255)
                                          tccgatctGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGcc                                                    <  2:384730/93‑1 (MQ=255)
                                          ctgtTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGcc                                                    >  2:1592068/1‑100 (MQ=255)
                                          ctgtTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGcc                                                    <  2:1560325/100‑1 (MQ=255)
                                           tgtTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAACCAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCct                                                   >  1:1619425/1‑100 (MQ=255)
                                           tgtTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCct                                                   >  2:1044884/1‑100 (MQ=255)
                                           tgtTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCct                                                   >  1:664849/1‑100 (MQ=255)
                                           ngtTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCct                                                   >  2:1687817/2‑100 (MQ=255)
                                            gtTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCctg                                                  <  2:653490/100‑1 (MQ=255)
                                               aTTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCctgctg                                               <  1:1864827/100‑1 (MQ=255)
                                                  gTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGccagatcgga                                            >  1:384731/1‑92 (MQ=255)
                                                  gTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCt                                            >  2:1857664/1‑100 (MQ=255)
                                                  gTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCt                                            >  2:1961085/1‑100 (MQ=255)
                                                  gTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCt                                            <  2:1042666/100‑1 (MQ=255)
                                                  gTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCt                                            <  2:66193/100‑1 (MQ=255)
                                                     cTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTa                                         >  2:888873/1‑100 (MQ=255)
                                                      tGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTAc                                        <  2:657797/100‑1 (MQ=255)
                                                      tGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTAc                                        >  2:953641/1‑100 (MQ=255)
                                                      tGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTAc                                        <  1:862165/100‑1 (MQ=255)
                                                      tGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTAc                                        >  1:792946/1‑100 (MQ=255)
                                                       gggTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTAct                                       >  2:460229/1‑100 (MQ=255)
                                                       gggTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTAct                                       <  2:400282/100‑1 (MQ=255)
                                                        ggTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTActc                                      <  1:1156665/100‑1 (MQ=255)
                                                        ggTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTActc                                      <  1:210062/100‑1 (MQ=255)
                                                        ggTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTActc                                      <  1:466225/100‑1 (MQ=255)
                                                            atgaATCCGGCTACGAAATGGGCGACGTGCAGAAAACCAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGtt                                  <  2:1649443/100‑1 (MQ=255)
                                                             tgaATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTc                                 >  1:559578/1‑100 (MQ=255)
                                                              gagtccgtgtACGAAATGGGCTACGTGCAGAAAAACAAACTGGCTGCTAGTTAATCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTc                                 <  1:680934/90‑1 (MQ=255)
                                                              gaATCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCg                                <  1:480630/100‑1 (MQ=255)
                                                                aTCCGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGc                              <  1:67198/100‑1 (MQ=255)
                                                                  ccGGCTACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTTCGCCTGCTGCCTTTACTCTGTTCGGCa                             <  2:881884/99‑1 (MQ=255)
                                                                    ggCTACGAAATGGGCGACGTGCAGAAACACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTc                          >  1:1270453/1‑100 (MQ=255)
                                                                       tACGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTg                       <  2:1102919/100‑1 (MQ=255)
                                                                        aCGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGa                      >  1:727876/1‑100 (MQ=255)
                                                                         cGAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGAt                     <  1:1044889/100‑1 (MQ=255)
                                                                          gAAATGGGCGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATc                    >  2:434974/1‑100 (MQ=255)
                                                                                  cGACGTGCAGAAAACCAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAagag            <  1:526267/100‑1 (MQ=255)
                                                                                  cGACGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAagag            <  1:1592080/100‑1 (MQ=255)
                                                                                    aCGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAgagga          >  2:2018621/1‑100 (MQ=255)
                                                                                     ngTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAgaggag         >  2:316613/2‑100 (MQ=255)
                                                                                     cGTGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAgaggag         <  2:940799/100‑1 (MQ=255)
                                                                                       tGCAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAgaggagcc       >  2:1828419/1‑98 (MQ=255)
                                                                                         cAGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAGAGGAGACGa     >  2:1405842/1‑100 (MQ=255)
                                                                                          aGAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAGAGGAGACGaa    <  1:460230/100‑1 (MQ=255)
                                                                                           gAAAAACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTTCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAGAGGAGACGAAc   >  1:404321/1‑100 (MQ=255)
                                                                                            aaaaACAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAGAGGAGACGAACa  <  1:953645/100‑1 (MQ=255)
                                                                                                |                                                                                               
CCTTTGCTATCGCTGCCAGCTTCGGTATGGCTGCTGTTCTGTCTGTTATTGTTCTGGGTGATGAATCCGGCTACGAAATGGGCGACGTGCAGAAAACCAAACTGGCTGCTATTGAAGCCGAGTGGGAAACGCAACCTGCGCCTGCTGCCTTTACTCTGTTCGGCATTCCTGATCAGGAAGAGGAGACGAACA  >  NC_000913/772113‑772304

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: