Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,391,417 C→G 7.5% W30S (TGG→TCG)  nuoN ← NADH:quinone oxidoreductase subunit N

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,391,4170CG 7.5% 125.1 / 2.3 47W30S (TGG→TCG) nuoNNADH:quinone oxidoreductase subunit N
Reads supporting (aligned to +/- strand):  ref base C (18/25);  new base G (2/2);  total (20/27)
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.53e-01

TAAAGCATGGCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCACTGCTGGATAT  >  NC_000913/2391279‑2391558
                                                                                                                                          |                                                                                                                                             
tAAAGCATGGCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGc                                                                                                                                           <  2:44289/143‑1 (MQ=255)
 aaaGCATGGCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCACACCGGCGAGGAAATAACTGCGTCGACGCGCa                                                                                                                                          >  2:279348/1‑143 (MQ=255)
    gCATGGCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTGGCGACGCAAtc                                                                                                                                       <  1:812208/143‑2 (MQ=255)
      aTGGCGAAACCATCAGCGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGga                                                                                                                                     <  1:857745/143‑1 (MQ=255)
              gccgtCAACGCGCAGCAGCGGCGGGACGTCCATAGCGCCCGCCGGGCCAACAAACCAGAGCGAAACCAGCGGCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATca                                                                                                                             <  2:815524/139‑1 (MQ=255)
                cATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATcaca                                                                                                                           <  1:36566/143‑1 (MQ=255)
                  tCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCCCCTGGCCAACAAACCAGAGCCAAACCAGCGCCGCGTTAACCCCAACAACCGAGACCGCAGCGTTGAGGAAAGGATCGCGTCCCCACGCAATCGAGCGCACCACAac                                                                                                                         >  2:484872/1‑143 (MQ=255)
                     aCGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAaccac                                                                                                                      <  1:581465/143‑1 (MQ=255)
                            tCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAGTAACCGGGAGCGTCGCGTTGAGGAGATGCATGCCTCGCCCCGCAGTGGAGAGGCTCGCAAGCACCGCCGTc                                                                                                               >  2:722049/1‑143 (MQ=255)
                              agtggcgTAACGTCCATAGCGCCTGCCTGGCCAACAAACCAGAGCGAAACCAGCGTCGCGTTAAGCCCAATCATCCACAGCGTAGCGTTGATGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCaa                                                                                                             <  2:404962/140‑1 (MQ=255)
                              aGCGGCGTAACGCCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCaa                                                                                                             >  1:630744/1‑143 (MQ=255)
                                      aaCGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCcgacga                                                                                                     <  2:519438/143‑1 (MQ=255)
                                          tCCATAGCGCGCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAg                                                                                                 <  2:708746/143‑1 (MQ=255)
                                            cATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGGCGagcagcc                                                                                               >  1:242503/1‑142 (MQ=255)
                                              tAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGACGGCGTAGCGTGGCGGAAATGATTGAGTCGCCACGCAATGGAGAGCATCACAACCACCACCGCCAAGCCGACGAGCacaaac                                                                                             >  2:18144/1‑138 (MQ=255)
                                               cgcgcCCGCATCGCCAGCACGCCAGAGCGACACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCTCCACCGTCAATCCGACGATCAGCAACg                                                                                            <  2:769166/142‑1 (MQ=255)
                                                 cgcccgccTGGCCAACAACCCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATACCCGAGTGCGTAGCGGTGAGGAAATGATCGCGGCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGt                                                                                          <  2:23377/143‑1 (MQ=255)
                                                  gcccgccTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTa                                                                                         <  1:721097/143‑1 (MQ=255)
                                                  gcccgccTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGAGTAGCGTTGAGGTAATGATTGCGGCGCCACGCCATGGAGAGCCTCACAACCCCCACCGTCAAGCCGACCGTCAGCAACGGTa                                                                                         >  1:335565/1‑143 (MQ=255)
                                                  gaccgcaTGGCCAACAAACCAGAGCTAAACTAGCGCCGCGTGAAGCCCAATAACCGAGAGCGTAGCGTGGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTa                                                                                         <  2:695287/141‑1 (MQ=255)
                                                    ccgccTGGCCAACAAACCAGAGCGAAACCAGCGCCCAGTTAAGCCCAATAACCAAGAGCGGAGCGTCGCGGTAATGATGGCGTCGCCACGCAATGGAGAAGATCACAACCAACACCGACAGACCGACGATCAGCAACgcaaac                                                                                       >  2:570460/1‑138 (MQ=255)
                                                                       gagCTAAGCCAGCGCCGGGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCCGGTTTtgtg                                                                    <  2:494287/143‑1 (MQ=255)
                                                                            aaaCCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGACCATCACAACCACCCCCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGCGGAgct                                                               >  1:873636/1‑141 (MQ=255)
                                                                                    tccgcGTTCAACCCAATAACCGAGAGCGTAGCGTTCAGGAAATCATTGCGTCGCCACGCAATGGATAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCAt                                                       <  2:88272/142‑1 (MQ=255)
                                                                                    gccgcTTTAAGACCAATCACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGCGAGCATCACAACCACCACCGTCAAGCCGACGATCATCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGtttt                                                       <  2:199187/143‑4 (MQ=255)
                                                                                     ccgcTTTAAGCCCAATAATCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATg                                                      <  1:675869/143‑1 (MQ=255)
                                                                                      cgcgTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATgg                                                     >  1:876799/1‑143 (MQ=255)
                                                                                             aGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTa                                              >  1:674909/1‑143 (MQ=255)
                                                                                              gCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGAAGCCACGCAATGGAGAGCATAACAACCACCACCGTCAAGCCGACGAACAGCAGCGGTAGAAGTGCGAGCAGGCTTTGTGGAGTTAATGTCATGGCGATTTaa                                             >  2:810598/1‑142 (MQ=255)
                                                                                               cccAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACg                                            >  1:878154/1‑143 (MQ=255)
                                                                                                  aaTAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCACTGGAGAGCATCACAACCACCACCGGCAAGCCGACGATCAGCAACGGTAGCAGCGAGATCCGGTTGTGTGGAGTTATAGTCATGGCGATTTACGGca                                         >  2:594600/1‑142 (MQ=255)
                                                                                                       ccGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTgta                                    <  1:622783/143‑1 (MQ=255)
                                                                                                        cgggagCGTGGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTgtag                                   <  2:879493/140‑1 (MQ=255)
                                                                                                          agagCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTgtagta                                 <  1:631039/143‑1 (MQ=255)
                                                                                                          agagCGTAGCGTTGAGGAAAGGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTgtagta                                 <  1:574464/143‑1 (MQ=255)
                                                                                                          agagCGGAGCGTTGGGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGGCGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTgtagta                                 <  2:894180/143‑1 (MQ=255)
                                                                                                                    gtTGAGGAAATGGTGGCGGCGCGACGCAATGGGGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAAt                       <  2:4105/143‑1 (MQ=255)
                                                                                                                    gtTGAGGAAATGATTGCGTCGCCACGCAATCGAGAGCATCACAACCACCACCGTCAAGCCGAAGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAAt                       <  1:450847/143‑1 (MQ=255)
                                                                                                                     ttGAGGAAGTGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAAtt                      <  1:103703/143‑1 (MQ=255)
                                                                                                                          gAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTGCGGCCGTGTAGTAGTAACGGCAGTAACaa                 >  1:216728/1‑143 (MQ=255)
                                                                                                                            ctccaTTGCGTGGCGACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAAc               >  1:136775/5‑143 (MQ=255)
                                                                                                                             tccaTTGCGTGGCGACGCAATGGAGAGCATCACAACCACCACCGTAAAGCCGACGATCAGCAACGGTAGCAGTGCAATCAGGATTTGTGGAGGTATTGTCAGGACGATTGACGGCCTGGTAGTAGGAACGGAAATAACAAAcc              >  2:46897/4‑143 (MQ=255)
                                                                                                                              tGATTGCGTCGCCACGCAATGGAGAGCATCAAAACAAACACCGTAAAGACAACGATAAGAAACTGTAAAAAAACAATCAAGTTTTATAGAATTATTGTAATAAAAATTAACGAAATTGTAATAATAAAGGAAAAAACaaaaaa             >  2:519260/1‑140 (MQ=255)
                                                                                                                               gATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGCTCAGGGTTTGCGGAGATATTGTCATGGCGCTTTACGGCCTTGTAGTAGTAACGGAATGAACAAACCAc            >  2:285085/1‑143 (MQ=255)
                                                                                                                                 ttGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTCTTGTGGAGTTATTGTCATGGCGATTGACCGCCTTGTAGTAGTAACGGAATTAACAAACCActg          >  2:207817/1‑143 (MQ=255)
                                                                                                                                  tgccccGCCACGCAACGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCActgc         <  2:111528/138‑1 (MQ=255)
                                                                                                                                  ggcgacaccACGAAATGGAGAGCATCACAACCACCGCCGGCGCGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTGTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCActgc         <  2:53762/136‑1 (MQ=255)
                                                                                                                                      tCGCCACGCAATGGAGAGCATCACAGCCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCACTGCTGGa     <  2:242503/143‑1 (MQ=255)
                                                                                                                                         ccACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATAGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCACTGCTGGatat  >  2:32463/1‑143 (MQ=255)
                                                                                                                                          |                                                                                                                                             
TAAAGCATGGCGAAACCATCAACGCGCATCAGCGGCGTAACGTCCATAGCGCCCGCCTGGCCAACAAACCAGAGCGAAACCAGCGCCGCGTTAAGCCCAATAACCGAGAGCGTAGCGTTGAGGAAATGATTGCGTCGCCACGCAATGGAGAGCATCACAACCACCACCGTCAAGCCGACGATCAGCAACGGTAGCAGTGCGATCAGGTTTTGTGGAGTTATTGTCATGGCGATTTACGGCCTTGTAGTAGTAACGGAATTAACAAACCACTGCTGGATAT  >  NC_000913/2391279‑2391558

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