Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 774,064 C→A 8.4% T341T (ACC→ACA)  cydB → cytochrome d terminal oxidase, subunit II

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913774,0640CA 8.4% 166.0 / 3.8 55T341T (ACC→ACA) cydBcytochrome d terminal oxidase, subunit II
Reads supporting (aligned to +/- strand):  ref base C (27/23);  new base A (3/2);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.45e-01

TGTTCTCCTCCCTGACGCTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGCCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGTTCGGTCGTATCACCAAAGAAGATATTGAACGTAACACCCACTCTCTGTACTAAGTAAGGAGCTAAAAATGTGGTATTTCG  >  NC_000913/773919‑774208
                                                                                                                                                 |                                                                                                                                                
tgtTCTCCTCCCTGACGCTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGCCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTgg                                                                                                                                               <  6:153623/149‑1 (MQ=255)
     tcctccCTGACGCTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGCCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGc                                                                                                                                          >  8:110771/1‑149 (MQ=255)
                gCTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGCCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTgg                                                                                                                               >  8:111026/1‑149 (MQ=255)
                gCTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGCCGGCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTgg                                                                                                                               <  4:8193/149‑1 (MQ=255)
                    gCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGCCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTa                                                                                                                             <  7:110771/147‑1 (MQ=255)
                    gCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGCCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTAc                                                                                                                            <  3:172529/148‑1 (MQ=255)
                       tGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGCCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACATGGGTTGCGGTGGTTCTGGTACCGAt                                                                                                                        >  5:52991/1‑149 (MQ=255)
                       tGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGCCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACATGGGTTGCGGTGGTTCTGGTACCGAt                                                                                                                        >  5:187883/1‑149 (MQ=255)
                         catcatCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGCCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACATGGGTTGCGGTGGTTCTGGTACCGATc                                                                                                                       <  6:187883/148‑1 (MQ=255)
                         catcatCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGCCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACATGGGTTGCGGTGGTTCTGGTACCGATc                                                                                                                       <  6:52991/148‑1 (MQ=255)
                             atcCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGCCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCt                                                                                                                  >  6:57027/1‑149 (MQ=255)
                              tcCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGCCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAAAGACATGACATGGGTTGCGGTGGTACTGGTACAGATCATTCTg                                                                                                                 >  6:97590/1‑149 (MQ=255)
                                              cGCAATGTTCCCGTTTGTGATGCCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGgt                                                                                                 >  6:174518/1‑149 (MQ=255)
                                               gCAATGTTCCCGTTTGTGATGCCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGgtg                                                                                                <  8:7608/149‑1 (MQ=255)
                                               gCAATGTTCCCGTTTGTGATGCCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGgt                                                                                                 <  6:81638/148‑1 (MQ=255)
                                                 aaTGTTCCCGTTTGTGATGCCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGgtgt                                                                                               >  2:150655/1‑148 (MQ=255)
                                                       cccGTTTGTGATGCCGTCCAGCACCATGATGAACGAAATTCTGACAATGTGGGATGCAACTTACAGCCAGCTGACGCTTAACGTCATGACCTGTGTTGCTGTGGTTATGGTACAGATCATTCTGCTCTACACCTACTGGTGTTACTgg                                                                                         >  1:34559/1‑148 (MQ=255)
                                                        ccGTTTGTGATGCCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGaa                                                                                       <  3:91484/149‑1 (MQ=255)
                                                        ccGTTTGTGATGCCGTCCAGCACCATGATGAACGCAAGTCTGACAAGGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACCCCGCCTGGGGTTACTGGa                                                                                        >  5:56775/1‑148 (MQ=255)
                                                          gTTTGTGATGCCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGaaaa                                                                                     >  4:78363/1‑149 (MQ=255)
                                                              gtgATGCCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGtt                                                                                 >  1:184821/1‑149 (MQ=255)
                                                               tgATGCCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGtt                                                                                 <  3:202528/148‑1 (MQ=255)
                                                                             aCCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGtt                                                                                 <  4:60925/134‑1 (MQ=255)
                                                                             aCCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGtt                                                                                 >  3:60925/1‑134 (MQ=255)
                                                                               cATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCGTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGAGCATTCTGCTCTACACCGCCTGGTGTTAGTGGAAAATGTTCGGTCGTATCACCa                                                                   >  2:186615/1‑146 (MQ=255)
                                                                                   atgaACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGTTCg                                                                               <  7:187989/130‑1 (MQ=255)
                                                                                   atgaACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGTTCg                                                                               >  8:187989/1‑130 (MQ=255)
                                                                                   atgaACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGTTCGGTCGTATCACCAaag                                                                >  5:218416/1‑145 (MQ=255)
                                                                                     gaACGCAAGTCTGACAATGTGGGATGCAACTACCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTATGGTACCGATCATTCTGCTCTACACCGCCTCGTGTTACTGGAAAATGTTCGGTCGTATCACCAAAGAAGat                                                           >  6:571/1‑148 (MQ=255)
                                                                                      aaCGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTGACGGCATGGCCGGGGTTGCGGTGGTTCTGGTACCGCTCCTTCTGCTCTACCCCGCCTGGTGTTACTGGAAAATGTTCGGTCCGATCACGAGAGAAGataa                                                         >  6:195923/1‑148 (MQ=255)
                                                                                           aaGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGTTCGGTCGTATCACCAAAGAAGATATTGAAc                                                    >  6:186829/1‑149 (MQ=255)
                                                                                                     aTGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGTTCGGTCGTATCACCAAAGAAGATATTGAACGTAACACCCa                                          <  7:111026/149‑1 (MQ=255)
                                                                                                      tgtgGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGTTCGGTCGTATCACCAAAGAAGATATTGAACGTAACACCCAc                                         >  7:246127/1‑149 (MQ=255)
                                                                                                         gggATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGTTCGGTCGTATCACCAAAGAAGATATTGAACGTAACACCCActct                                      >  3:73617/1‑149 (MQ=255)
                                                                                                          ggATGCAACTTCCAGCCAGCTGACTCTTAACGTCATGACCTGGTTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGTTCGGTCGTATCACCAAAGAAGATATTGAACGTAACACCCActct                                      <  5:97590/148‑1 (MQ=255)
                                                                                                                     ccagccagCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGTTCGGTCGTATCACCAAAGAAGATATTGAACGTAACACCCACTCTCTGTACtaagt                           <  8:246127/148‑1 (MQ=255)
                                                                                                                      cagccagcTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGTTCGGTCGTATCACCAAAGAAGATATTGAACGTAACACCCACTCTCTGTACtaagtaa                         <  1:150655/149‑1 (MQ=255)
                                                                                                                      cagcaagcTGACGCTTAACGTCATGAGCTGTCTTGCGGTGGTTCTGGTACCGATCATTCTTCTCTACACCGCCTGGTGTTACTGGAATATGTTCGGTCGTATCACCTAAGAAGATTTTGAACGTAACACCCACTCTCTGTACtaagtaa                         <  7:111569/149‑1 (MQ=255)
                                                                                                                        gccagcTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGTTCGGTCGTATCACCAAAGAAGATATTGAACGTAACACCCACTCTCTGTACtaagtaag                        >  5:54535/1‑148 (MQ=255)
                                                                                                                            gcTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGTTCGGTCGTATCACCAAAGAAGATATTGAACGTAACACCCACTCTCTGTACTAAGTAAGGAGc                    >  1:80982/1‑148 (MQ=255)
                                                                                                                              tGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATg                                                                                   >  7:96443/1‑83 (MQ=255)
                                                                                                                              tGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATg                                                                                   <  8:96443/83‑1 (MQ=255)
                                                                                                                              tGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATg                                                                                   <  6:83498/83‑1 (MQ=255)
                                                                                                                              tGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATg                                                                                   >  5:83498/1‑83 (MQ=255)
                                                                                                                                aCGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGTTCGGTCGTATCACCAAAGAAGATATTGAACGTAACACACACTCTCTGTACTAAGTAAGGAGCTaaa                >  1:31165/1‑148 (MQ=255)
                                                                                                                                  gCTTAAAGTCATGACCTGGGTTGCGGTTGTTCTTGTACCGATCATTCTGCTCTACACCGACTTGTGTTACTGGAAAATGTTCTGTCGTATCACCAAAGAAGATATTGAACGTAACACCCACTCTCTGTACTAAGTAAGGAGCTAAAAAt             <  2:34559/149‑1 (MQ=255)
                                                                                                                                      aaCGTCATTACCTGGGTTGCGGTGGTTCTGGTACCTATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGTTCGGTCGTATCACCAAAGAAGATATTGAACGTAACACCCACTCTCTGTACTAAGTAAGGAGCTAAAAAtgtg          <  5:57027/148‑1 (MQ=255)
                                                                                                                                      aaCGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGTTCGGTCGTATCACCAAAGAAGATATTGAACGTAACACCCACTCTCTGTACTAAGTAAGGAGCTAAAAAtgtg          <  2:184821/148‑1 (MQ=255)
                                                                                                                                      aaCGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGTTCGGTCGTATCACCAAAGAAGATATTGAACGTAACACCCACTCTCTGTACTAAGTAAGGAGCTAAAAAtgtg          <  1:2883/148‑1 (MQ=255)
                                                                                                                                       aCGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGTTCGGTCGTATCACCAAAGAAGATATTGAACGTAACACCCACTCTCTGTACTAAGTAAGGAGCTAAAAATGTGGt        <  5:195923/149‑1 (MQ=255)
                                                                                                                                          tCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGTTCGGTCGTATCACCAAAGAAGATATTGAACGTAACACCCACTCTCTGTACTAAGTAAGGAGCTAAAAATGTGGTAtt     >  8:55903/1‑149 (MQ=255)
                                                                                                                                          tCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCAGCTCTACACCGCCTCGTGTTACTGGAAAATGTTCGGTCGTATCACCAAAGAAGATATTGAACGTAACACCCACTCTCTGTACTAAGTAAGGAGCTAAAAATGTGGTAtt     >  8:132744/1‑149 (MQ=255)
                                                                                                                                             tGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGTTCGGTCGTATCACCAAAGAAGATATTGAACGTAACACCCACTCTCTGTACTAAGTAAGGAGCTAAAAATGTGGTATTTCg  <  5:174518/149‑1 (MQ=255)
                                                                                                                                              gACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGTTCGGTCGTATCACCAAAGAAGATATTGAACGTAACACCCACTCTCTGTACTAAGTAAGGAGCTAAAAATGTGGTATTTCg  >  5:218392/1‑148 (MQ=255)
                                                                                                                                              gACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGTTCGGTCGTATCACCAAAGAAGATATTGAACGTAACACCCACTCTCTGTACTAAGTAAGGAGCTAAAAATGTGGTATTTCg  >  8:127593/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
TGTTCTCCTCCCTGACGCTGGCCTGCATCATCCTGACAGCCGGTATCGCAATGTTCCCGTTTGTGATGCCGTCCAGCACCATGATGAACGCAAGTCTGACAATGTGGGATGCAACTTCCAGCCAGCTGACGCTTAACGTCATGACCTGGGTTGCGGTGGTTCTGGTACCGATCATTCTGCTCTACACCGCCTGGTGTTACTGGAAAATGTTCGGTCGTATCACCAAAGAAGATATTGAACGTAACACCCACTCTCTGTACTAAGTAAGGAGCTAAAAATGTGGTATTTCG  >  NC_000913/773919‑774208

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