Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,120,082 G→A 19.0% A152T (GCG→ACG)  adhB → alcohol dehydrogenase cytochrome c subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,120,0820GA19.0% 38.5 / 4.6 21A152T (GCG→ACG) adhBalcohol dehydrogenase cytochrome c subunit
Reads supporting (aligned to +/- strand):  ref base G (6/11);  new base A (2/2);  total (8/13)
Fisher's exact test for biased strand distribution p-value = 6.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.71e-01

ACTTCGACCGCGCCGTGCGCCAGGGTGTCGCGCCGGGTGGCCGACGGCTGTACCCGGCGATGCCTTACCCGTCCTACGCCAAGCTCAGCGATGATGATGTCAAGGCGCTGTATGCGTTCTTCATGCACGGCGTGCAACCGGCGCGGCAGGCCAACCTGGGCAGCGACATCCCCTGGCCGCTGAACCTGCGCTGGCCCATCGCCTTGTGGAACGGCCTGTTCGCCGCCACTACGCCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATAT  >  NC_002947/4119942‑4120227
                                                                                                                                            |                                                                                                                                                 
aCTTCGACCGCGCCGTGCGCCAGGGTGTCGCGCCGGGTGGCCGACGGCTGTACCCGGCGATGCCTTACCCGTCCTACGCCAAGCTCAGCGATGATGATGTCAAGGCGCTGTATGCGTTCTTCATGCACGGCGTGCAACCGGCGCGGCAgg                                                                                                                                          <  2:425590/150‑1 (MQ=255)
       ccGCGCCGTGCTCCAGGGTGTCGCGCCGGGTGGCCGACGGCTGTACCCGGCGATGCCTTACCCGTCCTACGCCAAGCTCAGCGATGATGATGTCAAGGCGCTGTATGCGTTCTTCATGCACGGCGTGCAACCGGCGCGGCAGGCCAACCt                                                                                                                                   <  1:148454/150‑1 (MQ=255)
        cgcgCCGTGCGCCAGGGTGTCGCGCCGGGTGGCCGACGGCTGTACCCGGCGATGCCTTACCCGTCCTACGCCAAGCTCAGCGATGATGATGTCAAGGCGCTGTATGCGTTCTTCATGCACGGCGTGCAACCGGCGCGGCAGGCCAACCTg                                                                                                                                  <  2:229499/150‑1 (MQ=255)
           gcCGTGCGCCAGGGTGTCGCGCCGGGTGGCCGACGGCTGTACCCGGCGATGCCTTACCCGTCCTACGCCAAGCTCAGCGATGATGATGTCAAGGCGCTGTATGCGTTCTTCATGCACGGCGTGCAACCGGCGCGGCAGGCCAACCTGGGc                                                                                                                               >  2:657182/1‑150 (MQ=255)
                        gtgtCGCGCCGGGTGGCCGACGGCTGTACCCGGCGATGCCTTACCCGTCCTACGCCAAGCTCAGCGATGATGATGTCAAGGCGCTGTATGCGTTCTTCATGCACGGCGTGCAACCGGCGCGGCAGGCCAACCTGGGCAGCGACATCCCCt                                                                                                                  <  1:398792/150‑1 (MQ=255)
                           tCGCGCCGGGTGGCCGACGGCTGTACCCGGCGATGCCTTACCCGTCCTACGCCAAGCTCAGCGATGATGATGTCAAGGCGCTGTATGCGTTCTTCATGCACGGCGTGCAACCGGCGCGGCAGGCCAACCTGGGCAGCGACATCCCCTGGc                                                                                                               <  1:171624/150‑1 (MQ=255)
                               gcCGGGTGGCCGACGGCTGTACCCGGCGATGCCTTACCCGTCCTACGCCAAGCTCAGCGATGATGATGTCAAGGCGCTGTATGCGTTCTTCATGCACGGCGTGCAACCGGCGCGGCAGGCCAACCTGGGCAGCGACATCCCCTGGCCGCt                                                                                                           <  2:563353/150‑1 (MQ=255)
                                           aCGGCTGTACCCGGCGCTGGACTCCCTGTCGTACGCCAATCTCAGCGATGATGATGTCAAGTCTTTTTATCCTTTCTTCATGCACGGCGTGCAACCGGCGCGGCAGGCCAACCTGGGCAGCGACATCCCCTGGCCGCTGAACCTGCGCTg                                                                                               <  1:203193/150‑1 (MQ=255)
                                                tGTACCCGGCGATGCCTTACCCGTCCTACGCCAAGCTCAGCGATGATGATGTCAAGGCGCTGTATGCGTTCTTCATGCACGGCGTGCAACCGGCGCGGCAGGCCAACCTGGGCAGCGACATCCCCTGGCCGCTGAACCTGCGCTGGCCCa                                                                                          >  2:414103/1‑150 (MQ=255)
                                                     ccGGCGATGCCTTACCCGTCCTACGCCAAGCTCAGCGATGATGATGTCAAGGCGCTGTATGCGTTCTTCATGCACGGCGTGCAACCGGCACGGCAGGCCAACCTGGGCAGCGACATCCCCTGGCCGCTGAACCTGCGCTGGCCCATCGcc                                                                                     <  2:527233/150‑1 (MQ=255)
                                                                    ccGTCCTACGCCAAGCTCAGCGATGATGATGTCAAGGCGCTGTATGCGTTCTTCATGCACGGCGTGCAACCGGCGCGGCAGGCCAACCTGGGCAGCGACATCCCCTGGCCGCTGAACCTGCGCTGGCCCATCGCCTTGTGGAACGGCCTg                                                                      <  1:211733/150‑1 (MQ=255)
                                                                    ccGTCCTACGCCAAGCTCAGCGATGATGATGTCAAGGCGCTGTATGCGTTCTTCATGCACGGCGTGCAACCGGCGCGGCAGGCCAACCTGGGCAGCGACATCCCCTGGCCGCTGAACCTGCGCTGGCCCATCGCCTTGTGGAACGGCCTg                                                                      >  1:18946/1‑150 (MQ=255)
                                                                              ccAAGCTCAGCGATGATGATGTCAAGGCGCTGTATGCGTTCTTCATGCACGGCGTGCAACCGGCGCGGCAGGCCAACCTGGGCAGCGACATCCCCTGGCCGCTGAACCTGCGCTGGCCCATCGCCTTGTGGAACGGCCTGTTCGCCGCCa                                                            >  1:120441/1‑150 (MQ=255)
                                                                                                  gTCAAGGCGCTGTATGCGTTCTTCATGCACGGCGTGCAACCGACGCGGCAGGCCAACCTGGGCAGCGACATCCCCTGGCCGCTGAACCTGCGCTGGCCCATCGCCTTGTGGAACGGCCTGTTCGCCGCCACTACGCCTTACACCGACAAg                                        >  1:301756/1‑150 (MQ=255)
                                                                                                  gTCAAGGCGCTGTATGCGTTCTTCATGCACGGCGTGCAACCGACGCGGCAGGCCAACCTGGGCAGCGACATCCCCTGGCCGCTGAACCTGCGCTGGCCCATCGCCTTGTGGAACGGCCTGTTCGCCGCCACTACGCCTTACACCGACAAg                                        >  1:301712/1‑150 (MQ=255)
                                                                                                            tGTATGCGTTCTTCATGCACGGCGTGCAACCGACGCGGCAGGCCAACCTGGGCAGCGACATCCCCTGGCCGCTGAACCTGCGCTGGCCCATCGCCTTGTGGAACGGCCTGTTCGCCGCCACTACGCCTTACACCGACAAGGCCGGGCAgg                              <  2:301712/150‑1 (MQ=255)
                                                                                                            tGTATGCGTTCTTCATGCACGGCGTGCAACCGACGCGGCAGGCCAACCTGGGCAGCGACATCCCCTGGCCGCTGAACCTGCGCTGGCCCATCGCCTTGTGGAACGGCCTGTTCGCCGCCACTACGCCTTACACCGACAAGGCCGGGCAgg                              <  2:301756/150‑1 (MQ=255)
                                                                                                                         cATGCACGGCGTGCAACCGGCGCGGCAGGCCAACCTGGGCAGCGACATCCCCTGGCCGCTGAACCTGCGCTGGCCCATCGCCTTGTGGAACGGCCTGTTCGCCGCCACTACGCCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGaa                 <  1:414103/150‑1 (MQ=255)
                                                                                                                           tGCACGGCGTGCAACCGGCGCGGCAGGCCAACCTGGGCAGCGACATCCCCTGGCCGCTGAACCTGCGCTGGCCCATCGCCTTGTGGAACGGCCTGTTCGCCGCCACTACGCCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAAcc               >  1:267596/1‑150 (MQ=255)
                                                                                                                             cACGGCGTGCAACCGGCGCGGCAGGCCAACCTGGGCAGCGACATCCCCTCGCCGCTGAACCTGCGCTGGCCCATCGCCTTGTGGAACGGCCTGTTCGCCGCCACTACGCCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGt             <  1:657182/150‑1 (MQ=255)
                                                                                                                                        aCCGGCGCGGCAGGCCAACCTGGGCAGCGACATCCCCTGGCCGCTGAACCTGCGCTGGCCCATCGCCTTGTGGAACGGCCTGTTCGCCGCCACTACGCCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTatat  >  1:630554/1‑150 (MQ=255)
                                                                                                                                            |                                                                                                                                                 
ACTTCGACCGCGCCGTGCGCCAGGGTGTCGCGCCGGGTGGCCGACGGCTGTACCCGGCGATGCCTTACCCGTCCTACGCCAAGCTCAGCGATGATGATGTCAAGGCGCTGTATGCGTTCTTCATGCACGGCGTGCAACCGGCGCGGCAGGCCAACCTGGGCAGCGACATCCCCTGGCCGCTGAACCTGCGCTGGCCCATCGCCTTGTGGAACGGCCTGTTCGCCGCCACTACGCCTTACACCGACAAGGCCGGGCAGGATGCACAGTGGAACCGTGGCGCCTATAT  >  NC_002947/4119942‑4120227

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