Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 448,014 G→A S11S (TCC→TCT)  cyoD ← cytochrome bo3 subunit 4

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913448,0140GA96.2% 69.8 / NA 26S11S (TCC→TCT) cyoDcytochrome bo3 subunit 4
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (15/10);  minor base C (1/0);  total (16/10)
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.52e-01

GCACGGTGAAGACAAACGCCGTCATGTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCCAAACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACCG  >  NC_000913/447791‑448223
                                                                                                                                                                                                                               |                                                                                                                                                                                                                 
gCACGGTGAAGACAAACGCCGTCATGTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGAGACGCGCCGCTGTGATCGGTAGAATg                                                                                                                                                                                         >  1:28966/1‑250 (MQ=255)
                     tCATGTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACAAGAACCTGTACCAATCCCATTGCAAGGATTGTTCCCACAATTACGGCCGGAGAGGCAGCTCCTGTCATCCCCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGAGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATc                                                                                                                                                                    >  1:226294/1‑250 (MQ=255)
                                               gTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCGCAGAATTACGGCCGGAGAGGCAGCGCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGGCAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCTATGAGACGCGCCGCGTTGAGCGTTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATc                                                                                                                                         <  1:81792/251‑1 (MQ=255)
                                                                                ngAACCTGAAACACTGCCATTGCCAGGATTGTTCCCAAACTTACGGCCGGAGACAAAGCTCCTGTCATCACCATCCAGACCGGAATCACCGTCCGAATGATCGACAGGATAAAGCCTGACATGTAGGTTTTTACGCTGCCATGAGACGCGCCGCTGTGTTCGGTAGa                                                                                                                                                                                            >  1:90695/2‑167 (MQ=255)
                                                                                   aCCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGAGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGccc                                                                                                                                                                         >  1:261795/1‑183 (MQ=255)
                                                                                                       cAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTGTACGCTGCCATGAGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATcc                                                                                                                                        <  1:348350/196‑1 (MQ=255)
                                                                                                           aTTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGAATCCCCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGAGACGCGCCGCTGTGATCGGGAGAATGCCTCAGTACATCGCCCCCATCa                                                                                                                                                                   <  1:156323/165‑1 (MQ=255)
                                                                                                                  ccAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGAGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCCAAACCACATCCAGGAAGTGc                                                                                                                    <  1:4556/205‑1 (MQ=255)
                                                                                                                                     ggCAGCTCCTGTCATCACCATCCAGAACGGACTCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGAGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAAAAGTGAACACACAGATCCAAACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTg                                                   >  1:167177/1‑251 (MQ=255)
                                                                                                                                     ggCAGCTCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGCTAACGCCTGTCATGTAGGTTTTTACGCTGCCATGAGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCCAAACCACATCCAGGAAGTGCCAGAAc                                                                                                              >  1:366656/1‑192 (MQ=255)
                                                                                                                                            ncTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGAAGGTTTTTACGCTGCCATGAGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAg                                                                                                                                                                  >  1:308791/2‑133 (MQ=255)
                                                                                                                                            ccTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGAGACGCGCCGCTGTGATCGGTAGAATGTCTCATTACATCGCCCCCATCAg                                                                                                                                                                  <  2:308791/133‑1 (MQ=255)
                                                                                                                                                     aCCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGAGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCACCCATCAGATAAACAACAGTGAACACACAGATCCAAACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGAtt                                                     >  1:294400/1‑233 (MQ=255)
                                                                                                                                                     aCCATCCAGAACGGAATCACCGACAGAATGATCGACAGGATAAAGCCTGTCATGTAAGTTTTTACGCTGCCATGCGACGCGCCGCTGTGCTCGGTAGCATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCCAAACCACATCCAGGAAGTGCCAGAACAGGCTCCGGCACATGATGCggg                                                                                        >  1:467800/1‑198 (MQ=255)
                                                                                                                                                           cAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGAGACGCGCAGCTGTGATCGGTAGAATGACTCa                                                                                                                                                                                    >  1:180954/1‑100 (MQ=255)
                                                                                                                                                                         cGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTT‑ACGCTGCCAGGAGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCCAAACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCcgc                                                                <  1:230462/201‑1 (MQ=255)
                                                                                                                                                                            cAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGAGACGCGCCGCTGTGATCGGTAGAATGACTCATTCCATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCCAAACCACATCCAGGAAGTGCCTGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCa                                           >  1:497728/1‑220 (MQ=255)
                                                                                                                                                                              gAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGATGCCATGAGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAAcacaca                                                                                                                                             >  1:212739/1‑120 (MQ=255)
                                                                                                                                                                                  gNNCGACAGNNTAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGAGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATaa                                                                                                                                                              <  2:28966/99‑1 (MQ=255)
                                                                                                                                                                                     cGNCNGNATAAAGCCTGTCATGTAGGTTTTTACGCTNCCATGAGACGCGCCGCTGTGATCNNTANNNTCNNTCAtt                                                                                                                                                                                  <  2:365382/76‑1 (MQ=255)
                                                                                                                                                                                     cGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGAGACGCGCCGCTGTGATCGGTAGAATGaataatt                                                                                                                                                                                  >  1:365382/1‑70 (MQ=255)
                                                                                                                                                                                      gACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGAGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCCAAACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACcg  >  1:170971/1‑251 (MQ=255)
                                                                                                                                                                                                aGCCGGTCATGTAGGTTTTTACGCGGCCATGAGACGCGCCGCTGTGAGCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCCAAACCACATCCAGGAAGTGCCAGAACAg                                                                                                            <  1:71212/135‑1 (MQ=255)
                                                                                                                                                                                                                       cTGCCATGAGACGCGCCTTTGTGATTGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCCAAACCACATCCATGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCtn                                                                          <  1:60075/146‑2 (MQ=255)
                                                                                                                                                                                                                         gCCATGAGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCACCATAAGATAAACAACCGTGAACACAAAGATCCAACCCACATCCAGGACGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGcc                                                                  >  1:321890/1‑152 (MQ=255)
                                                                                                                                                                                                                          ccATGAGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCCAAACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCAc                                                >  1:338164/1‑169 (MQ=255)
                                                                                                                                                                                                                               |                                                                                                                                                                                                                 
GCACGGTGAAGACAAACGCCGTCATGTTCCAGCCTTCATCTGATTTGGTATTCATGTGCAGGAAGCACACCAGATGCACCAGAACCTGTACCACTGCCATTGCCAGGATTGTTCCCAGAATTACGGCCGGAGAGGCAGCTCCTGTCATCACCATCCAGAACGGAATCACCGTCAGAATGATCGACAGGATAAAGCCTGTCATGTAGGTTTTTACGCTGCCATGGGACGCGCCGCTGTGATCGGTAGAATGACTCATTACATCGCCCCCATCAGATAAACAACAGTGAACACACAGATCCAAACCACATCCAGGAAGTGCCAGAACAGGCTCAGGCACATGATGCGGGTACGGTTAGTGCTGGTCAGGCCGCGACGGGCGATTTGCACCATCAGCACCGCCATCCAGATAAGACCAGAAGTGACGTGCAGACCG  >  NC_000913/447791‑448223

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