Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,403,813 C→T 10.5% W11* (TGG→TGA)  nuoC ← NADH:quinone oxidoreductase subunit CD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,403,8130CT10.5% 114.5 / 3.7 48W11* (TGG→TGA) nuoCNADH:quinone oxidoreductase subunit CD
Reads supporting (aligned to +/- strand):  ref base C (32/10);  new base T (2/3);  total (34/14)
Fisher's exact test for biased strand distribution p-value = 1.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

TCGCCAACTTCCAGTAATTGTTCACGCTTGATCCACACAACGGGAACCCCGGTGCGAGTCGCCTGAACAGTAAAGGCATCCGGCCCAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCAAATGGTCCGTGGTGACTGTCGGCGCAGTAAGCTTCGCGAAATAGATTTAGGAATAAGCGAAGAAAAATCCCTTTGCCGACAGGCGCAAATTAAA  >  NC_000913/2403673‑2403955
                                                                                                                                            |                                                                                                                                              
tCGCCAACTTCCAGTAATTGTTCACGCTTGAGCCACACAACGGGACCCCCGGGGCGAGTCGCCTGAACAGTAAAGGCATCCGGCCCAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCa                                                                                                                                              <  2:643221/143‑1 (MQ=255)
   ccAACTTGCAGTAATTGTTCAGGCTTAATCCACACAACGGGGACCCCGGTGCGAGTCGCCGGAACAGTAAAAACATACGGCCCAAAACGGTTGCGCAGTTCGCAAAGCACCGGATCATCAAGATGATCACGGGTCTGCCACGc                                                                                                                                           >  2:840377/1‑143 (MQ=255)
   ccAACTTCCAGTAATTGTTCACGCTTGATCCACACAACGGGAACCCCGGTGCGAGTCGCCTGAACAGTAAAGGCATCCGGCCCAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGGCCTGCCAGGc                                                                                                                                           >  2:592427/1‑143 (MQ=255)
           cAGTAATTGTTCACGCTTGATCCACACAACGGGAACCCCGGTGTGAGTCGCCTGAACAGTAAAGGCATCCGGCCCAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTCCtt                                                                                                                                   >  2:746778/1‑143 (MQ=255)
            ggTAGTTGTTCCCGGTTGATCGCCAGCACGGTGGCCCGGGGGCGAGTCGCCTGGAGAGTAAAGGCATCCGGCCGAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTg                                                                                                                                  <  2:253267/142‑1 (MQ=255)
              tAATTGTTCACGCTTGATCCACACAACGGGAACCCCGGTGCGAGTCGCCTGAACAGTAAAGGCGTCCCGCCCAAAACGGTTGCGCAGTTCGCCAATGACCGGATGATCAAGATGATCGCGGGGCGGCCAGGCGGATTcgggcg                                                                                                                                >  2:533065/1‑138 (MQ=255)
                aTTGTTCACGCGTGACCCCCACAACGGGATCCCCGGTGCGAGTCGCCTGAACAGTAAAGGCATCCGGCCCAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGcgcg                                                                                                                              <  1:837279/143‑1 (MQ=255)
                       aCGCTTGATCCACACAACGGGAACCCCGGTGCGAGTCGCCTGAACAGTAAAGGCATCCGGCCCAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGt                                                                                                                       >  1:138510/1‑143 (MQ=255)
                          cTTGATCCACACAACGGGAACCCCGGTGCGAGTCGCCTGAACAGTAAAGGCATCCGGCCCAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCgg                                                                                                                    >  1:544447/1‑143 (MQ=255)
                            tGAGCCACACAACGGGAGCCCCGGTGCGAGTCGCCTGAACAGTAAAGGCATCCGGCCCAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTc                                                                                                                  <  1:756338/143‑1 (MQ=255)
                               tCCACACAACGGGAACCCCGGTGCGAGTCGCCTGAACAGTAAAGGCATCCGGCCCAAAACGGTTGCGCAGTTCGCCAATCACCGGATCGTCTAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTAata                                                                                                               >  1:534241/1‑143 (MQ=255)
                               tCCACACAACGGGAACCCCGGTGCGAGTCGCCTGAACAGTAAAGGCATCCGGCCCAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGGCCGTGCGCGGTTAAGTCGGTCata                                                                                                               >  2:703267/1‑143 (MQ=255)
                                  acacaACGGGAACCCCGGTGATAGTCGCCTGAACAGTAAAGGCATCCGGCCCAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGGTCTGACAAACCGGGTCTTGCGCGGTTAAGTCGGTCATATTg                                                                                                            <  1:695728/143‑1 (MQ=255)
                                      aaCGGGAACCCCGGTGCGAGTCGCCTGAACAGTAAAGGCATCCGGCCCAAAACGGTTGCGCAGTTCGCCAACCACCGGCGCATCAAGATGATCGCGGGTCTGTCAGGCGGGCTCTTGCGCGGTTACGCCGGGCATATCGTTCa                                                                                                        >  2:234958/1‑143 (MQ=255)
                                        cGGGAGCCCCGTTGCGAGTCGCATGAACAGTAAAGTCATCCGGCCCAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGTTCTGCCAGTCGGGTTCTTGCGCGGTTAAGTCTGTCATATTGTTCAcc                                                                                                      <  2:354228/143‑1 (MQ=255)
                                        cGGGAACCCCGGTGCGAGTCGCCTGAACAGTAAAGGCATCCGGCCCAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCAGATTGTTCAcc                                                                                                      >  1:823511/1‑143 (MQ=255)
                                         gggAACCCCGGTGCGAGTCGCCTGAACAGTAAAGGCATCCGGCCCAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTGGGTCATATAGTGCACCa                                                                                                     >  1:631917/1‑143 (MQ=255)
                                            aaCCCCGGTGCGAGTCGCCTGAACAGTAAAGGCATCCGGCCCAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGGTCACCATTg                                                                                                  >  1:496949/1‑143 (MQ=255)
                                                        aGTCGCCTGAACAGTAAAGGCATCCGGCCCAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTGTTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCAAATGGTCCgt                                                                                      >  1:134703/1‑143 (MQ=255)
                                                             ccTGAACAGTAAAGGCATCCGGCCCACAACGGTTGCGCAGTTCGCCAATCACCGGATCATCATGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCAAATGGTCCGTGGTGa                                                                                 <  1:940560/143‑1 (MQ=255)
                                                              cTGAACAGTAAAGGCATCCGGCCCAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATCAATATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCATGGCAAATGGTACGTGGTGAc                                                                                >  1:296792/1‑143 (MQ=255)
                                                               tGAACAGTAAAGGCATCCGGCCCAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCACATGGTCCGTGGTGACt                                                                               >  1:678429/1‑143 (MQ=255)
                                                                gAACAGTAAAGGCATCCGGCCCAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCAAATGGTCCGTGGTGACTg                                                                              >  1:882417/1‑143 (MQ=255)
                                                                         aGGCATCCGGCCCAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCGTTGCAAATGGTCCGTGGTGACTGTCGGCGCAg                                                                     >  1:866/1‑143 (MQ=255)
                                                                          ggCATCCGGCCCAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCAAATGGTCCGTGGTGACTGTCCGCGCAGt                                                                    >  2:881226/1‑143 (MQ=255)
                                                                             aTCCGGGCCAAAACGGTGGCGCAGTTCGCCAATCACCGGATAATCAAGATGATCGCGGGTCTGTCAGGCGGGTTCTTGCGCGGTTAAGTCGGGCATATTGTGCACCATTGCAAATGGTCCGTGGTGACTGTCGGAGCAGTAAg                                                                 <  1:190247/143‑1 (MQ=255)
                                                                             aTCCGGCCCAAAGCGGCTGCGCAGGTCGCCAATCACCGGATCAGCAAGATGATCGCGGGTCTGTCAGGCGGGTCCGTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCAAATGGTCCGTGGTGACTGTCGGCGCAGTAAg                                                                 <  1:70404/143‑1 (MQ=255)
                                                                              tCCGGCCCAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATCTTGTTCACCATTGCCAATGGTCCGTGGTGACTGCCGGCGCAGTAAGc                                                                >  2:891501/1‑143 (MQ=255)
                                                                                ctgtccAAAACGGTTGCGCAGTTCGCCAATCACCGGGTCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCGTGCGCGGTTAAGTCGGTCATATTGTTCCCCATTGCAAATGGTCCGTGGTGACTGTCGGCGCAGTAAGCtt                                                              <  1:342290/139‑1 (MQ=255)
                                                                                 ggCCCAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATCGAGATGATCGCGGGTCTGCCAGGCGGGGTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCAAATGGTCCGTGGTGACTGTCGGCGCAGTAAGcgtg                                                             >  1:901130/1‑140 (MQ=255)
                                                                                 ggCCCAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCAAATGGTCCGTGGTGACTGTCGGCGCAGTAAGCTTc                                                             >  1:396490/1‑143 (MQ=255)
                                                                                   cccAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATTTAGAGGATCGCGGGTCTGTCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCAAATGGTCCGTGGTGACTGTCGGCGCAGTAAGCTTcgc                                                           <  1:731636/143‑1 (MQ=255)
                                                                                            gTTGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCAAATGGTCCGTGGTGACTGTCGGCGCAGTAAGCTTCGCGAAATAGAt                                                  >  1:940169/1‑143 (MQ=255)
                                                                                            gTTGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCAAATGGTCCGTGGTGACTGTCGGCGCAGTAAGCTTCGCGAAATAGAt                                                  >  1:560128/1‑143 (MQ=255)
                                                                                             ttGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCAAATGGTCCGTGGTGACTGTCGGCGCAGTAAGCTTCGCGAAATAGAtt                                                 >  1:854548/1‑143 (MQ=255)
                                                                                                cgcAGTTCGCCGAGCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCAAATGGTCCGTGGTGACTGTCGGCGCAGTAAGCTTCGCGAAATAGATTTAg                                              <  2:788204/143‑1 (MQ=255)
                                                                                                     ttCGCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCAAATGGTCCGTGGTGACTGTCGGCGCAGTAAGCTTCGCGAAATAGATTTAGGAATa                                         >  1:112776/1‑143 (MQ=255)
                                                                                                        gCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCAACTGGTCCGTTGTGACTGTCGGCGCAGTCAGCCCCGCGAAATACAGTTATCAATAAGc                                      >  2:658151/1‑143 (MQ=255)
                                                                                                        gCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCAAATGGTCCGTGGTGACTGTTGGCGCAGTAAGCTTCGCGAAATAGTTTTAGGAATAAGc                                      >  1:650958/1‑143 (MQ=255)
                                                                                                        gCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCAAATGGTCCGTGGTGACTGTCGGCGCAGGAAGCTTCGCGAAATAGATTTAGGAATAGGc                                      >  1:236512/1‑143 (MQ=255)
                                                                                                           aaTCACCGGATCATCAAGATGATCGCGGGTCTGTCAGGCGGGTTCTTGCGCGATTGAGTCGGTCATAGTGTTCACCATTGCAAATGGTCCGGGGAGACCGTCGGCGCAGCAGGCGTCGCGAAATCGATTTGGGGATAGGccag                                   >  2:767730/1‑140 (MQ=255)
                                                                                                           aaTCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCAAATGGTCCGTGGTGACTGTCGGCGCAGGAAGCTCCGCGACATCAATCTAGGAATAAGCgaa                                   >  2:227921/1‑143 (MQ=255)
                                                                                                                ccGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCAAATGGTCCGTGGTGACTGTCGGCGCAGTAAGCTTCGCGAAATAGATTTAGGAATAAGCGAAGaaaa                              >  1:150557/1‑143 (MQ=255)
                                                                                                                         caAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCAAATGGTCCGTGGTGACTGTCGGCGCAGTAAGCTTCGCGAAATAGATGTAGGACTAAGCCAAACGACATCCCTTTg                     >  2:444178/1‑143 (MQ=255)
                                                                                                                          gaGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCAAATGGTCCGTGGTGACTGTCGGCGCAGTAAGCTTCGCGAAATAGATTTAGGAATAAGCGAAGAAAAATCCCTTTGc                    <  1:488255/142‑1 (MQ=255)
                                                                                                                                   gcgGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCAAATGGTCCGTGGTGACTGTCGGCGCAGTAAGCTTCGCGAAATAGATTTAGGAATAAGCGAAGAAAAGTCGCTTTGCCGCCCGgcg           >  1:75243/1‑143 (MQ=255)
                                                                                                                                          tGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCAAATGGTCCGTGGTGACTGTCGGCGCAGTAAGCTTCGCGAAATAGATTTAGGAATAAGCGAAGAAAAATCCCTTTGCCGACAGGCGCAAATTa    <  1:225773/143‑1 (MQ=255)
                                                                                                                                           gCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCAAATGGTCCGTGGTGACTGTCGGCGCAGTAAGCTTCGCGAAATAGATTTAGGAATAAGCGAAGAAAAATCCCTTTGCCGACAGGCGCAAATTaa   >  2:560639/1‑143 (MQ=255)
                                                                                                                                            ccATGCCGGTCCTTGCCCGCTTAAGTCTGTCCTCTTGTTCACCACTGCAAAGGGTCCCTGGTGACTGTCGGCGCAGTAAGCTTCGCGAAATAGATTTAGGAATAAGCGAAGAAAAATCCCTTTGCCGACAGGCGCAAATTaaa  <  1:373214/143‑1 (MQ=255)
                                                                                                                                            |                                                                                                                                              
TCGCCAACTTCCAGTAATTGTTCACGCTTGATCCACACAACGGGAACCCCGGTGCGAGTCGCCTGAACAGTAAAGGCATCCGGCCCAAAACGGTTGCGCAGTTCGCCAATCACCGGATCATCAAGATGATCGCGGGTCTGCCAGGCGGGTTCTTGCGCGGTTAAGTCGGTCATATTGTTCACCATTGCAAATGGTCCGTGGTGACTGTCGGCGCAGTAAGCTTCGCGAAATAGATTTAGGAATAAGCGAAGAAAAATCCCTTTGCCGACAGGCGCAAATTAAA  >  NC_000913/2403673‑2403955

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