Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,489,609 C→G 9.7% G202G (GGC→GGG)  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,6090CG 9.7% 119.3 / 4.1 41G202G (GGC→GGG) PP_2184formate dehydrogenase subunit beta
Reads supporting (aligned to +/- strand):  ref base C (18/19);  new base G (2/2);  total (20/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 = 1.81e-01

CCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCTTCGCCATCGCACGTGACGCTGGCTACCTGG  >  NC_002947/2489465‑2489757
                                                                                                                                                |                                                                                                                                                    
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                               <  5:70457/122‑1 (MQ=255)
                                                                                                                                              ggCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCtt                               >  6:70457/1‑122 (MQ=255)
                                                                                                                                                cGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCTTCGCCATCGCACGTGACGCTGGCTACCTgg  <  7:191178/149‑1 (MQ=255)
                                                                                                                                                cGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCTTCGCCATCGCACGTGACGCTGGCTACCTgg  >  8:102535/1‑149 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
CCTGCGTGGCCGTGGTGGTGCGGCGTTCCCGGCCGGCATCAAGTGGCGTACCGTGCGTGATGCTGGCGCCGGGCAGAAGTACGTGGTGTGCAACGCCGACGAGGGTGACTCGGGCACGTTCGCTGATCGCATGCTGATGGAAGGCGACCCCTTCCTGCTGATCGAGGGCATGATCATTGCCGGCCTCGCCGTGGGTGCCGACAAGGGCTATATCTATGTGCGCTCGGAGTATCCCGATGCCATCCGGGTACTCGACCAGGCCTTCGCCATCGCACGTGACGCTGGCTACCTGG  >  NC_002947/2489465‑2489757

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