Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,646,976 C→T 17.0% Q141Q (CAG→CAA)  gcd ← quinoprotein glucose dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,646,9760CT17.0% 45.3 / 4.5 24Q141Q (CAG→CAA) gcdquinoprotein glucose dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (10/9);  new base T (2/2);  total (13/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.33e-01

GGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCA  >  NC_002947/1646828‑1647124
                                                                                                                                                    |                                                                                                                                                    
gggTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAAct                                                                                                                                                     >  2:495223/1‑150 (MQ=255)
    gATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAActggct                                                                                                                                                 <  1:475945/150‑1 (MQ=255)
        tGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAATTGGCTGGcc                                                                                                                                             >  1:475771/1‑150 (MQ=255)
        tGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAATTGGCTGGcc                                                                                                                                             >  1:475858/1‑150 (MQ=255)
                          cGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGaa                                                                                                                           >  2:493524/1‑150 (MQ=255)
                               gCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAATTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGAc                                                                                                                      <  2:475858/150‑1 (MQ=255)
                               gCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAATTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGAc                                                                                                                      <  2:475771/150‑1 (MQ=255)
                                  aTGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTCCTGTCGCGGCCCAGTTCGCCGAACACCTCACCGGGATGGGTGAAATGGCTGGCCAGCGCGCAAGCGCGCGAAGCGACCAc                                                                                                                   >  1:590060/1‑150 (MQ=255)
                                      tCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGTCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGcc                                                                                                               >  2:500442/1‑150 (MQ=255)
                                        ggTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCAc                                                                                                             >  2:334774/1‑150 (MQ=255)
                                         gTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACg                                                                                                            >  2:152896/1‑150 (MQ=255)
                                          tACGGCCGTAGGCCTGCCATTCGCCCTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGc                                                                                                           >  1:545397/1‑150 (MQ=255)
                                                     gCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAgcagtgca                                                                                                <  2:616307/150‑1 (MQ=255)
                                                        tGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGcagtgcagtg                                                                                             >  2:95498/1‑150 (MQ=255)
                                                                    tCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGcc                                                                                                               >  2:361356/1‑120 (MQ=255)
                                                                    tCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGcc                                                                                                               <  1:361356/120‑1 (MQ=255)
                                                                        gCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAActgg                                                                                                                                                   <  1:108880/80‑1 (MQ=255)
                                                                        gCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAActgg                                                                                                                                                   >  2:108880/1‑80 (MQ=255)
                                                                                            gCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAActgg                                                                                                                                                   <  1:529618/60‑1 (MQ=255)
                                                                                            gCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAActgg                                                                                                                                                   >  2:529618/1‑60 (MQ=255)
                                                                                                                  ccAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGTTTTGCCCGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAAcacca                                   >  2:478516/1‑150 (MQ=255)
                                                                                                                    aGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACg                                 <  1:495223/150‑1 (MQ=255)
                                                                                                                        cGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGa                             <  2:545397/150‑1 (MQ=255)
                                                                                                                          ccGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATg                           <  1:451274/150‑1 (MQ=255)
                                                                                                                           cGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATgg                          <  1:95498/150‑1 (MQ=255)
                                                                                                                                acTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAAc                     >  1:467414/1‑150 (MQ=255)
                                                                                                                                    cACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGa                 <  1:388732/150‑1 (MQ=255)
                                                                                                                                             gggTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATagccagcc        <  1:334774/150‑1 (MQ=255)
                                                                                                                                                tGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTg     >  1:460118/1‑150 (MQ=255)
                                                                                                                                                   aCTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCa  <  1:493524/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
GGGTGATCTGGCGCAGCGGCGAGTAGCGGTCGCCATGCTCGGTACGGCCGTAGGCCTGCCATTCGCCGTCGGGCATGGCCGGGGCGGCGCTGGCCATTTCGCTGCTGTCGCGGCCCAGTTCGCCGAACACTTCACCGGGATGGGTGAACTGGCTGGCCAGCGCGCAAGCGCCCGAAGCGACCACTGCCACGCCGAGCAGTGCAGTGTTGGCTTTGCTGGCTGGGCCGATCAGCGGGCGACGTGCCCACGGCAGCAGCAACACCACGCCGATGGCGAACCAGATAGCCAGCCGTGGCA  >  NC_002947/1646828‑1647124

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: ‑