Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,291,674 G→A A246V (GCC→GTC)  ccmH ← holocytochrome c synthase CcmH component

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,291,6740GA83.3% 40.0 / ‑3.6 18A246V (GCC→GTC) ccmHholocytochrome c synthase CcmH component
Reads supporting (aligned to +/- strand):  ref base G (3/0);  new base A (5/10);  total (8/10)
Fisher's exact test for biased strand distribution p-value = 6.86e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGTTGGCAGGTAAGAGTTTCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGG‑TTGTCGTTGGGATCAGA‑‑TGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGCCAACGCGGCCCAACATAATCCAGCCTTCTATATCTCCCGGATTTTTTTGCAGTTGAGTACGCATCCCCAGCGCAAGACGCGACATCTCTTCTTCGTTGAGCGGATCGGCTTTCGGATCCAGCGCCCTGTCCAGCAACGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAG  >  NC_000913/2291430‑2291916
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
aTCGTGGGCAGTTAAGAGTTTCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCTCTGCTGCGCAAAGGCATTAAACGCATACATGCTTAGCACACGTATATTGCTATGGTCCGTTCCCACCAGCGGACGTAGCAGTTCACCGCCGAGGCGG‑TTGTCGTTGGGATCAGA‑‑TGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGACGTAt                                                                                                                                                                                                                                               <  1:161861/250‑1 (MQ=255)
                               tcCCACGAGGCAACGGCTTCGCCAAATCGATGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGG‑TTGTCGTTGGGGTCCGA‑‑TGAACGAGTCAACGCTTCCGCGTATCCCAATGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGGGGCGTATTCATCGGTGGCGATACTGGCGTTACCCAgcgc                                                                                                                                                                                                               >  1:73464/1‑251 (MQ=255)
                                     gctgCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTTTGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGG‑TTGTCGTTGGGATCAGA‑‑TGAACGAGTCAACGCTTCAGCGTATCCCAGTTCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGACGTATGCAGCGGTGGCGATACTGTCGTTTCCCAGCTCCATgac                                                                                                                                                                                                         <  2:102694/248‑3 (MQ=255)
                                                      aaaTCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGG‑TTGTCGTTGGGATCAAA‑‑AGAACGAGTCAACGATTCAGCGTATCCCCGTGCAGAATCACTGTTTTTCGGATCGGAGCGATACGAAGTGGCGTATGCATCTGTGGCGATAACGGCGTTGCCCATCGCCATGCCAACGCCGCCCCACATa                                                                                                                                                                                         >  1:119896/1‑250 (MQ=255)
                                                        aTCGATGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGGGTAGT‑GTTGGGAGCAGC‑‑TGAACGAGTTGACGCGTAACCGTAACCCAGTGCAGCATCACAGTTTTTCGGCTCGAGGCGGTACGTCGGGACGTATGCCGCGGTGGGGATTATTGCGTGACCCAGCGCCaaagcaagtcggcgccacattatc                                                                                                                                                                                      >  1:132136/1‑228 (MQ=255)
                                                               ctgctCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGG‑TTGTCGTTTGGATCCAA‑‑AGAACGAGTCAACGCTTCAGCGTATCCCCGTGCAGCCCTACTGGTTTTCGGATCGAGGCGATACGAAGTGAAGTCTGCATCGGTGGGGCTACTGGCCGTACCCAGCGCCATGCCAACGCGGCCCAACCTAATcccacctt                                                                                                                                                                               >  2:128673/1‑245 (MQ=255)
                                                                    cAAAGGCATTAAACGCATACATTCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCGGTTCACCGCCGAGGCGG‑TTGTCGTTGGGATCAGA‑‑TGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGACGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGCCAACGCGGCCCAACATAATCCAGCCTTCtata                                                                                                                                                                          <  2:65748/251‑1 (MQ=255)
                                                                                               gCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGG‑TTGTCGTTGGGGTCAGA‑‑TGAACGAGTCAACGCCTCAGCGGATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGACGGATGCCTCGGGGGCGATACTGGCGTTAACCAGCGCCATGCCAAACCGGCCCAACATAATACCGCCTTCTATATCTCCCGGGttttttttccgtggagt                                                                                                                                                >  2:15150/1‑240 (MQ=255)
                                                                                                         atTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGG‑TTGTCGTTGGGATCAGA‑‑TGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGACGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGCCAACGCGGCCCAACATAATCCAGCCTTCTATATCTCCCGGATTTTTTTGCAGTTGAGTACGCATcccc                                                                                                                                      <  1:207730/250‑1 (MQ=255)
                                                                                                            gCTAGGGTCCGTTCTCACCAGCTGAAGTAGCAGTTCACCGCCGAGGCGG‑TTGTCGTTGGGATCAGA‑‑TGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTGTTTCGGATCGAGTCGATACGCAGTGACGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGCCAACGCGGCCCAACATAATCCAGCCTTCTATATCTCCCGGATTTTTTTGCAGTTGAGTACGCATCCCCAgcg                                                                                                                                  <  2:98317/251‑1 (MQ=255)
                                                                                                                  gTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGG‑TTGTCCTTGGGATCAGA‑‑TGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGACGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGCCAACGCGGCCCAACATAATCCAGCCTTCTATATCTCCCGGATTTTTTTGCAGTTGAGTACGCATCCCCAGCGCAAGAc                                                                                                                            <  2:210619/251‑1 (MQ=255)
                                                                                                                        tctcACCAGCTGACGTAGCAGTTCACCGCCGAGGCGG‑TTGGCGGTGGG‑‑CAGAACTAAGCGAATCAACACTTCAGCGTCTCCCAGGGCAGCGTCACTGTTTTTTTGGTATAGGCGATGCACAGTGGCCGATTCATCGGGGGCGATACTGGACTTACCCAGCGCCATGCCAACGAGGCCCCACATAATCCACCCTTCTTTATATCCCGGATTTTTTTTCAGGTGTATACGGATCCCCCGCGCAAGAcgccac                                                                                                                       >  1:115842/1‑247 (MQ=255)
                                                                                                                                                                   gTTGGGATCAGA‑‑TGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGACGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGCCAACGCGGCCCAACATAATCCAGCCTTCTATATCTCCCGGATTTTTTTGCAGTTGAGTACGCATCCCCAGCGCAAGACGCGACATCTCTTCTTCGTTGAGCGGATCGGCTTTCGGATCCAGCGccc                                                                            >  2:28429/1‑251 (MQ=255)
                                                                                                                                                                                                      gTATCCCAGTGCAGCATCACTGTTTTTCGGAGCGAGGCGATACGCAGTGACGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGCCAACGCGGCCCAACATAATCCAGCCTTCTATATCTCCCGGATTTTTTTGCAGTTGAGTACGCATCCCCAGCGCAAGACGCGACATCTCTTCTTCGTTGAGCGGATCGGCTTTCGGATCCAGCGCCCTGTCCAGCAACGCCGGAGCCTGTGCCGTTGCct                                           <  2:195725/251‑1 (MQ=255)
                                                                                                                                                                                                        aTCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGACGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGCCAACGCGGCCCAACATAATCCAGCCTTCTATATCTCCCGGAGTTTTTTGCAGTTGATTACGCATCCCCAGCGCAAGACGCGACATCTCTTCTTCGTTGAGCGGATCGGCTTTCGGATCCAGCGCCCTGTCCAGCAACGCCGGAGCCTGTGCCGTGGCctg                                          <  1:242017/250‑1 (MQ=255)
                                                                                                                                                                                                                                               aCGCAGTGACGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGCCAACGCGGCCCAACATAATCCAGCCTTCTATATCTCCCGGATTTTTTTGCAGTTGAGTACGCATCCCCAGCGCAAGACGCGACATCTCTTCTTCGTTGAGCGGATCGGCTTTCGGATCCAGCGCCCTGTCCAGCAACGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTTCCAGTCTGgtag  <  2:209986/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 gCAGTGACGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGCCAACGCGGCCCAACATAATCCAGCCTTCTATATCTCCCGGATTTTTTTGCAGTTGAGTACGCATCCCCAGCGCAAGACGCGACATCTCTTCTTCGTTGAGCGGATCGGCTTTCGGATCCAGCGCCCTGTCCAGCAACGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGcc            <  1:78485/239‑1 (MQ=255)
                                                                                                                                                                                                                                                 gCAGTGACGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGCCAACGCGGCCCAACATAATCCAGCCTTCTATATCTCCCGGATTTTTTTGCAGTTGAGTACGCATCCCCAGCGCAAGACGCGACATCTCTTCTTCGTTGAGCGGATCGGCTTTCGGATCCAGCGCCCTGTCCAGCAACGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGcc            >  2:78485/1‑239 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                  
ATCGTTGGCAGGTAAGAGTTTCAACATCATCTCCCACGCGGCAACGGCTTCGCCAAATCGCTGCTGCTCAAAGGCATTAAACGCATACATGCTTAGCACACGGATATTGCTATGGTCCGTTCTCACCAGCTGACGTAGCAGTTCACCGCCGAGGCGG‑TTGTCGTTGGGATCAGA‑‑TGAACGAGTCAACGCTTCAGCGTATCCCAGTGCAGCATCACTGTTTTTCGGATCGAGGCGATACGCAGTGGCGTATGCATCGGTGGCGATACTGGCGTTACCCAGCGCCATGCCAACGCGGCCCAACATAATCCAGCCTTCTATATCTCCCGGATTTTTTTGCAGTTGAGTACGCATCCCCAGCGCAAGACGCGACATCTCTTCTTCGTTGAGCGGATCGGCTTTCGGATCCAGCGCCCTGTCCAGCAACGCCGGAGCCTGTGCCGTGGCCTGCTGCCAGATTTTCACCTGCTGATAATTGCCAGTCTGGTAG  >  NC_000913/2291430‑2291916

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