Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,120,891 G→A 10.3% A421A (GCG→GCA 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,8910GA10.3% 103.0 / 3.9 39A421A (GCG→GCAadhBalcohol dehydrogenase cytochrome c subunit
Reads supporting (aligned to +/- strand):  ref base G (18/17);  new base A (2/2);  total (20/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.50e-01

TATCACGCTCAATGGTTCGCAGCGGGTGGTGGCTGCAGGGGTGCCGGACGCCTATCGCATGCCGGCTTTGCGCGAGCAGTTGAGTGATCAGGAGATTGCCGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCCAGCAGCAGCCCTATCATCTTGCATATGCGTTGAATTGGAGTGGCACATGTACCGACGCGAGGCTGCCCCCAACAAAGACCGCGCAGCACATGAGCCAC  >  NC_002947/4120744‑4121037
                                                                                                                                                   |                                                                                                                                                  
tatCACGCTCAATGGTTCGCAGCGGGTGGTGGCTGCAGGGGTGCCGGACGCCTATCGCATGCCGGCTTTGCGCGAGCAGTTGAGTGATCAGGAGATTGCCGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGt                                                                                                                                                  <  2:238211/150‑1 (MQ=255)
    aCGCTCAATGGTTCGCAGCGGGTGGTGGCTGCAGGGGTGCCGGACGCCTATCGCATGCCGGCTTTGCGCGAGCAGTTGAGTGATCAGGAGATTGCCGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGAt                                                                                                                                              >  2:107107/1‑150 (MQ=255)
         cAATGGTTCGCAGCGGGTGGTGGCTGCAGGGGTGCCGGACGCCTATCGCATGCCGGCTTTGCGCGAGCAGTTGAGTGATCAGGAGATTGCCGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCa                                                                                                                                         >  2:324953/1‑150 (MQ=255)
                   cAGCGGGTGGTGGCTGCAGGGGTGCCGGACGCCTATCGCATGCCGGCTTTGCGCGAGCAGTTGAGTGATCAGGAGATTGCCGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGc                                                                                                                               >  2:695796/1‑150 (MQ=255)
                     gCGGGTGGTGGCTGCAGGGGTGCCGGACGCCTATCGCATGCCGGCTTTGCGCGAGCAGTTGAGTGATCAGGAGATTGCCGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCaa                                                                                                                             <  1:273391/150‑1 (MQ=255)
                     gCGGGTGGTGGCTGCAGGGGTGCCGGACGCCTATCGCATGCCGGCTTTGCGCGAGCAGTTGAGTGATCAGGAGATAGCCGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCaa                                                                                                                             <  1:695796/150‑1 (MQ=255)
                        ggtggtGGCTGCAGGGGTGCCGGACGCCTATCGCATGCCGGCTTTGCGCGAGCAGTTGAGTGATCAGGAGATTGCCGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCt                                                                                                                          <  1:324953/150‑1 (MQ=255)
                             tggCTGCAGGGGTGCCGGACGCCTATCGCATGCCGGCTTTGCGCGAGCAGTTGAGTGATCAGGAGATTGCCGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTg                                                                                                                     >  2:761497/1‑150 (MQ=255)
                             tggCTGCAGGGGTGCCGGACGCCTATCGCATGCCGGCTTTGCGCGAGCAGTTGAGTGATCAGGAGATTGCCGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTg                                                                                                                     <  1:540127/150‑1 (MQ=255)
                              ggCTGGAGGGGTTCCGGACGCCTATCGCATGCCGGCTTTGCGCGAGCAGTTGAGTGATCAGGAGATTGCCGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTgg                                                                                                                    <  2:643269/150‑1 (MQ=255)
                                                cGCCTATCGCATGCCGGCTTTGCGCGAGCAGTTGAGTGATCAGGAGATTGCCGGCGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCcag                                                                                                  <  1:761497/150‑1 (MQ=255)
                                                     aTCGCATGCCGGCTTTGCGCGAGCAGTTGAGTGATCAGGAGATTGCCGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCcagcagca                                                                                             >  1:622491/1‑150 (MQ=255)
                                                       cGCATGCCGGCTTTGCGCGAGCAGTTGAGTGATCAGGAGATTGCCGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCCagcagcagc                                                                                           <  2:196132/150‑1 (MQ=255)
                                                              cGGCTTTGCGCGAGCAGTTGAGTGATCAGGAGATTGCCGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCCAGCAGCAGCCCTatca                                                                                    >  2:254319/1‑150 (MQ=255)
                                                               ggCTTTGTGCGAGCAGTTGAGTGATCAGGAGATTGCCGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCCAGCAGCAGCCCTatcat                                                                                   >  1:468812/1‑150 (MQ=255)
                                                                     gcgcGAGCAGTTGAGTGATCAGGAGATTGCCGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCCAGCAGCAGCCCTATCATCTTGCa                                                                             >  2:586257/1‑150 (MQ=255)
                                                                                   gTGATCAGGAGATTGCCGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCCAGCAGCAGCCCTATCATCTTGCATATGCGTTGAAttc                                                               <  2:546776/150‑2 (MQ=255)
                                                                                       tCAGGAGATTGCCGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCCAGCAGAAGACCTAGCTTGGTGCATATGCGTTGAATTGGAGt                                                           >  2:498388/1‑150 (MQ=255)
                                                                                       tCAGGAGATTGCCGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTACGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCCAGCAGCAGCCCTATCATCTTGCATATGCGTTGAATTGGAGt                                                           >  2:263425/1‑150 (MQ=255)
                                                                                         aGGAGATTGCCGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCAGTGGATGCGCAGGCGGTAGGCAAGCTGCGTgg                                                                                                                    <  1:695028/91‑1 (MQ=255)
                                                                                         aGGAGATTGCCGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCAGTGGATGCGCAGGCGGTAGGCAAGCTGCGTgg                                                                                                                    >  2:695028/1‑91 (MQ=255)
                                                                                           gagaTTGCCGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCCAGCAGCAGCCCTATCATCTTGCATATGCGTTGAATTGGAGTGGca                                                       >  2:4242/1‑150 (MQ=255)
                                                                                           gagaTTGCCGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCCAGCAGCAGCCCTATCATCTTGCATATGCGTTGAATTGGAGTGGca                                                       >  2:55307/1‑150 (MQ=255)
                                                                                           gagaTTGCCGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCCAGCAGCAGCCCTATCATCTTACATATGCGTTGAATTGGAGTGGca                                                       <  1:254319/150‑1 (MQ=255)
                                                                                                   cGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCCAGCAGCAGCCCTATCATCTTGCATATGc                                                                        >  2:711936/1‑125 (MQ=255)
                                                                                                   cGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCCAGCAGCAGCCCTATCATCTTGCATATGc                                                                        <  1:711936/125‑1 (MQ=255)
                                                                                                        tGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCAGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCCAGCAGCAGCCCTATCATCTTGCATATGCGTTGAATTGGa                                                             <  1:677444/131‑1 (MQ=255)
                                                                                                        tGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCAGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCCAGCAGCAGCCCTATCATCTTGCATATGCGTTGAATTGGa                                                             >  2:677444/1‑131 (MQ=255)
                                                                                                          cTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCCAGCAGCAGCCCTATCATCTTGCATATGCGTTGAATTGGAGTGGCACATGTACCGACGCGa                                        >  2:595668/1‑150 (MQ=255)
                                                                                                                tttGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCcagcagca                                                                                             >  1:67678/1‑91 (MQ=255)
                                                                                                                tttGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCcagcagca                                                                                             <  2:67678/91‑1 (MQ=255)
                                                                                                                      cGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCCAGCAGCAGCCCTATCATCTTGCATATGCGTTGAATTGGAGTGGCACATGTACCGACGCGAGGCTGCCCCcaa                            <  1:464469/150‑1 (MQ=255)
                                                                                                                       gTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCCAGCAGCAGCCCTATCATCTTGCATATGCGTTGAATTGGAGTGGCACATGTACCGACGCGAGGCTGCCCCcaac                           >  1:354524/1‑150 (MQ=255)
                                                                                                                       gTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCCAGCAGCAGCCCTATAATCTTGCATATGCGTTGAATTGGAGTGGCACATGTACCGACGCGa                                        >  1:378181/1‑137 (MQ=255)
                                                                                                                       gTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCCAGCAGCAGCCCTATAATCTTGCATATGCGTTGAATTGGAGTGGCACATGTACCGACGCGa                                        <  2:378181/137‑1 (MQ=255)
                                                                                                                            gCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCCAGCAGCAGCCCTATCATCTTGCATATGCGTTGAATTGGAGTGGCACATGTACCGACGCGAGGCTGCCCCCAACAAAGa                      <  1:595668/150‑1 (MQ=255)
                                                                                                                                 gggCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCCAGCAGCAGCCCTATCATCTTGCATATGCGTTGAATTGGAGTGGCACATGTACCGACGCGAGGCTGCCCCCAACAAAGACcgcg                 >  1:236784/1‑150 (MQ=255)
                                                                                                                                               gTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCCAGCAGCAGCCCTATCATCTTGCATATGCGTTGAATTGGAGTGGCACATGTACCGACGCGAGGCTGCCCCCAACAAAGACCGCGCAGCACATGAGCCa   <  1:586257/150‑1 (MQ=255)
                                                                                                                                                tGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCCAGCAGCAGCCGTATCATCTTGCATATGCGTTGAATTGGAGTGGCACATGTACCGACGCGAGGCTGCCCCCAACAAAGACCGCGCAGCACATGAGCCAc  <  2:622491/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
TATCACGCTCAATGGTTCGCAGCGGGTGGTGGCTGCAGGGGTGCCGGACGCCTATCGCATGCCGGCTTTGCGCGAGCAGTTGAGTGATCAGGAGATTGCCGACGTGCTCAGCTTTGTGCGTACGGCCTGGGGCAACCAGGGGGGTGCGGTGGATGCGCAGGCGGTAGGCAAGCTGCGTGGGCATACCGACCCGGCCAGCAGCAGCCCTATCATCTTGCATATGCGTTGAATTGGAGTGGCACATGTACCGACGCGAGGCTGCCCCCAACAAAGACCGCGCAGCACATGAGCCAC  >  NC_002947/4120744‑4121037

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