Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,489,613 C→G 9.5% P204A (CCC→GCC)  PP_2184 → formate dehydrogenase subunit beta

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,489,6130CG 9.5% 112.0 / 4.1 42P204A (CCC→GCC) PP_2184formate dehydrogenase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (19/19);  new base G (2/2);  total (21/21)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.29e-01

CCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCTTCGCCATCGCACGTGACGCTGGCTACCTGGGTG  >  NC_002947/2489465‑2489760
                                                                                                                                                    |                                                                                                                                                   
ccTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGGACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGAcc                                                                                                                                                     <  7:303256/149‑1 (MQ=255)
      tGGCCGGGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTcc                                                                                                                                               <  4:219430/149‑1 (MQ=255)
       ggCCGGGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGGGTCGCCTTcca                                                                                                                                              <  1:58858/149‑2 (MQ=255)
        gCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGGGTCGCCTTccat                                                                                                                                             <  3:161673/149‑3 (MQ=255)
                       cgTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATg                                                                                                                              <  6:331406/149‑1 (MQ=255)
                        gTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGa                                                                                                                             >  7:154690/1‑149 (MQ=255)
                         ttCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGAACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGa                                                                                                                             >  3:123984/1‑148 (MQ=255)
                                ccggcATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGcc                                                                                                                     <  4:343757/149‑1 (MQ=255)
                                      tCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCTGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTc                                                                                                               <  1:313016/149‑1 (MQ=255)
                                      tCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTc                                                                                                               >  3:135050/1‑149 (MQ=255)
                                        aaGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGc                                                                                                             <  4:328198/149‑1 (MQ=255)
                                         aGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCgg                                                                                                                   >  2:31626/1‑142 (MQ=255)
                                         aGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTACCTGCTGATCGAGGGCATGATCATTGCCgg                                                                                                                   <  1:31626/142‑1 (MQ=255)
                                               gTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGt                                                                                                      >  2:61173/1‑149 (MQ=255)
                                                         tgATCCTGGCGCCGGGCAGACGGCCGTGGTGGCCAACGCCTACGAGGGGTACACCGGCACGCTCGCTGATCGCAGTCGTATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAgg                                                                                            <  4:242928/149‑1 (MQ=255)
                                                              cTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTTATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTTGGTGCCGACAAGGGCtat                                                                                       <  8:154690/149‑1 (MQ=255)
                                                                 gcgcCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATc                                                                                    >  4:105365/1‑149 (MQ=255)
                                                                    ccGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTAt                                                                                 <  5:211651/149‑1 (MQ=255)
                                                                             aGTACGTGGTGTGCAACTCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTc                                                                         <  7:306575/148‑1 (MQ=255)
                                                                               tACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGa                                                                      <  4:288476/149‑1 (MQ=255)
                                                                                    ggtgTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTAt                                                                  <  4:123984/148‑1 (MQ=255)
                                                                                    ggtgTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATc                                                                 >  7:131828/1‑149 (MQ=255)
                                                                                          cAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGAt                                                            <  1:61173/148‑1 (MQ=255)
                                                                                          cAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATg                                                           >  7:176806/1‑149 (MQ=255)
                                                                                                    gaGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCggg                                                  >  8:191178/1‑148 (MQ=255)
                                                                                                      gggTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTAc                                               >  3:302347/1‑149 (MQ=255)
                                                                                                             tCGGGCACGGTCGCTTATCCCATGCTTATGGAAGCCGACCCCTTCCTGCTGATCCCGGGCATGATCATGGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTAtgt                                                                               <  2:80813/110‑1 (MQ=255)
                                                                                                               gggCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAgg                                      <  4:302347/149‑1 (MQ=255)
                                                                                                               gggCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCAGCTGATCGAGGGAATGATCATTGCCGGCATCGCCGTGGGTCCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCacg                                      >  5:315563/1‑147 (MQ=255)
                                                                                                                    cGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCTTc                                 <  8:176806/149‑1 (MQ=255)
                                                                                                                     gTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCTTCg                                >  7:264198/1‑149 (MQ=255)
                                                                                                                           tGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCTTCGCCATCg                          >  3:309530/1‑149 (MQ=255)
                                                                                                                            gATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCTTCGCCATCg                          >  4:31004/1‑148 (MQ=255)
                                                                                                                                 cATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTCCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCTTCGCCATCGCACGt                     >  6:250359/1‑148 (MQ=255)
                                                                                                                                      tGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCTTCGCCATCGCACGTGACGCt               >  2:230835/1‑149 (MQ=255)
                                                                                                                                      agcaggAAGGGGTCGCCTTCCTGCTGATTGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCTTCGCCATCGCACGTGACGCt               >  3:70943/5‑149 (MQ=255)
                                                                                                                                         tGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCTTCGCCATCGCACGTGACGCTgg             <  3:31004/148‑1 (MQ=255)
                                                                                                                                         aggAAGGGGTCGCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCTTCGCCATCGCACGTGACGCTGGc            >  1:208451/2‑149 (MQ=255)
                                                                                                                                           gAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCTTCGCCATCGCACGTGACGCTGGCTa          >  8:288461/1‑149 (MQ=255)
                                                                                                                                              ggCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCtt                                  >  6:70457/1‑122 (MQ=255)
                                                                                                                                              ggCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCtt                                  <  5:70457/122‑1 (MQ=255)
                                                                                                                                                cGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCTTCGCCATCGCACGTGACGCTGGCTACCTgg     <  7:191178/149‑1 (MQ=255)
                                                                                                                                                cGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCTTCGCCATCGCACGTGACGCTGGCTACCTgg     >  8:102535/1‑149 (MQ=255)
                                                                                                                                                   ccccTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCTTCGCCATCGCACGTGACGCTGGCTACCTGGGTg  <  1:230835/149‑1 (MQ=255)
                                                                                                                                                    cccTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAg                                       <  2:282276/111‑1 (MQ=255)
                                                                                                                                                    cccTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAg                                       >  1:282276/1‑111 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCTTCGCCATCGCACGTGACGCTGGCTACCTGGGTG  >  NC_002947/2489465‑2489760

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