Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,981,016 C→T 11.1% G450D (GGC→GAC)  PP_2607 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,981,0160CT11.1% 115.7 / 6.0 45G450D (GGC→GAC) PP_2607hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (19/21);  new base T (2/3);  total (21/24)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.51e-01

AACCGAGGTCGCCGGTGTGAATGCGGATGGCGTAGAACGGCCCCTTGGTAAGGGGTGCCAAGCATGGGTTGGGCGTATGTTGTGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGCCGATGGCCAGGGCCAGTGCCTGCGATGTGGTGGCCCGGTACAGGTAGCCTGCGTCGATCAACGCCTGGTTGTTGAC  >  NC_002947/2980870‑2981146
                                                                                                                                                  |                                                                                                                                  
aaCCGAGGTCGCCGGTGTGAATGCGGATGGCGTAGAACGGCCCCTTGGTAAGGGGTGCCAAGCATGGGTTGGGCGTATGTTGTGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGcc                                                                                                                                   <  2:267791/148‑1 (MQ=255)
  ccGAGGTCGCCGGTGTGAATGCGGATGGCGTAGAACGGCCCCTTGGTAAGGGGTGCCAAGCATGGGTTGGGCGTATGTTGTGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCgg                                                                                                                                 >  8:162770/1‑148 (MQ=255)
   cGAGGTCGCCGGTGTGAATGCGGATGGCGTAGAACAGCCCCTTGGTAAGGGGTGCCAAGCATGGGTTGGGCGTATGTTGTGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCa                                                                                                                               >  5:277897/1‑149 (MQ=255)
                   aaTGCGGATGGCGTAGAACGGCCCCTTGGTAAGGGGTGCCAAGCATGGGTTGGGCGTATGTTGTGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGTCGGCATCGCTGTTGAACtgt                                                                                                                >  4:242689/1‑148 (MQ=255)
                   aaTGCGGATGGCGTAGAACGGCCCCTTGGTAAGGGGTGCCAAGCATGGGTTGGGCGTATGTTGTGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGTCGGCATCGCTGTTGAACtgt                                                                                                                >  4:238506/1‑148 (MQ=255)
                      gCGGATGGCGTAGAACGGCCCCTTGGTAAGGGGTGCCAAGCATGGGTTGGGCGTATGTTGTGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGa                                                                                                              >  3:69429/1‑147 (MQ=255)
                           tGGCGTAGAACGGCCCCTTGGTAAGGGGGGCCAAGCATGGGTTGGGCGTATGTTGTGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGtt                                                                                                       <  6:107265/149‑1 (MQ=255)
                             gCGTAGAACGGCCCCTTGGTAAGGGGTGCCAAGCATGGGTTGGGCGTATGTTGTGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGttt                                                                                                      <  7:162770/148‑1 (MQ=255)
                                tAGAACGGCCCCTTGGTAAGGGGTGCCAAGCATGGGTTGGGCGTATGTTGTGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTgcg                                                                                                   >  5:72760/1‑148 (MQ=255)
                                         cccTTGGTAAGGGGTGCCAAGCATGGGTTGGGCGTATGTTGTGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGttt                                                                                                      >  1:173302/1‑136 (MQ=255)
                                         cccTTGGTAAGGGGTGCCAAGCATGGGTTGGGCGTATGTTGTGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGttt                                                                                                      <  2:173302/136‑1 (MQ=255)
                                          ccTTGGTAAGGGGTGCCAAGCATGGGTTGGGCGTATGTTGTGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCtgt                                                                                        <  5:155081/149‑1 (MQ=255)
                                          ccTTGGTAAGGGGTGCCAAGCATGGGTTGGGCGTATGTTGTGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCtgt                                                                                        <  4:69429/149‑1 (MQ=255)
                                              ggTAAGGGGTGCCAAGCATGGGTTGGGCGTATGTTGTGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGa                                                                                     <  1:1282/148‑1 (MQ=255)
                                                    gggTGCCAAGCATGGGTTGGGCGTATGTTGTGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGc                                                                              >  2:216081/1‑149 (MQ=255)
                                                         ccAAGCATGGGTTGGGCGTATGTTGTGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGCCGATg                                                                         >  4:113480/1‑149 (MQ=255)
                                                             gCATGGGTTGGGCGTATGTTGTGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGTCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGCCGATggcca                                                                     <  1:95106/149‑1 (MQ=255)
                                                              cATGGGTTGGGCGTATGTTGTGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGCCGATggccag                                                                    >  6:181118/1‑149 (MQ=255)
                                                              cATGGGTTGGGCGTATGTTGTGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTATTTGCGCAAACACCTGTGGATCGACGCCGATggccag                                                                    >  4:220679/1‑149 (MQ=255)
                                                                  ggTTGGGCGTATGTTGTGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGTCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGCCGATggccagggc                                                                 <  3:242689/148‑1 (MQ=255)
                                                                  ggTTGGGCGTATGTTGTGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGTCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGCCGATggccagggc                                                                 <  3:238506/148‑1 (MQ=255)
                                                                      gggCGTATGTTGTGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACg                                                                               <  8:228227/130‑1 (MQ=255)
                                                                            aTGTTGTGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGCCGATGGCCAGGGCCAGTGCCTGc                                                       <  6:190165/148‑1 (MQ=255)
                                                                               ttgtGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCtgt                                                                                        >  3:79412/1‑112 (MQ=255)
                                                                               ttgtGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCtgt                                                                                        <  4:79412/112‑1 (MQ=255)
                                                                                  tgGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGCCGATGGCCAGGGCCAGTGCCTGCGAtgtg                                                 <  5:181972/148‑1 (MQ=255)
                                                                                  tgGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGCCGATGGCCAGGGCCAGTGCCTGCGAtgtg                                                 >  6:181972/1‑148 (MQ=255)
                                                                                   gggTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCtgt                                                                                        >  6:321918/1‑108 (MQ=255)
                                                                                   gggTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCtgt                                                                                        <  5:321918/108‑1 (MQ=255)
                                                                                        gCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACATGTGGATCGACGCCGATGGCCAGGGCCAGTGCCTGCGATGTGGTGGccc                                          >  6:91938/1‑149 (MQ=255)
                                                                                          ccATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGCCGATGGCCAGGGCCAGTGCCTGCGATGTGGTGGCCCgg                                        >  5:307648/1‑149 (MQ=255)
                                                                                                 gCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTTTTTGCGCAAACACCTGTGGATCGACGCCGATGGCCAGGGCCAGTGCCTGCGATGTGGTGGCCCGGTACAGGt                                 <  3:220679/149‑1 (MQ=255)
                                                                                                 gCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGGTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGCCGATGGCCAGGGCCAGTGCCTGCGATGTGGTGGCCCGGTACAgg                                  <  6:277897/148‑1 (MQ=255)
                                                                                                  cGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGCCGATGGCCAGGGCCAGTGCCTGCGATGTGGTGGCCCGGTACAGGTa                                <  6:307648/149‑1 (MQ=255)
                                                                                                                  cgccTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGCCGATGGCCAGGGCCAGTGCCTGCGATGTGGTGGCCCGGTACAGGTAGCCTGCGTCGATCAAc                <  3:113480/149‑1 (MQ=255)
                                                                                                                           cAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGCCGATggccagg                                                                   <  6:316479/89‑1 (MQ=255)
                                                                                                                           cAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGCCGATggccagg                                                                   >  5:316479/1‑89 (MQ=255)
                                                                                                                           cAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGCCGATGGCCAGGGCCAGTGCCTGCGATGTGGTGGCCCGGTACAGGTAGCCTGCGt                        >  4:39239/1‑132 (MQ=255)
                                                                                                                           cAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGCCGATGGCCAGGGCCAGTGCCTGCGATGTGGTGGCCCGGTACAGGTAGCCTGCGTCGATCAACGCCTGgttg       >  2:104394/1‑149 (MQ=255)
                                                                                                                            aaaGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGCCGATGGCCAGGGCCAGTGCCTGCGATGTGGTGGCCCGGTACAGGTAGCCTGCGTCGATCAACGCCTGgttgt      <  6:72760/149‑1 (MQ=255)
                                                                                                                             aaGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGCCGATGGCCAGGGCCAGTGCCTGCGATGTGGTGGCCCGGTACAGGTAGCCTGCGTCGATCAACGCCTGgttgtt     <  1:104394/149‑1 (MQ=255)
                                                                                                                              aGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGCCGATGGCCAGGGCCAGTGCCTGCGATGTGGTGGCCCGGTACAGGTAGCCTGCGTCGATCAACGCCTGGttgttg    <  5:181118/149‑1 (MQ=255)
                                                                                                                                        gATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGCCGATGGCCAGGGCCAGTGCCTGCGATGtggtgg                                             >  2:247297/1‑98 (MQ=255)
                                                                                                                                        gATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGCCGATGGCCAGGGCCAGTGCCTGCGATGtggtgg                                             <  1:247297/98‑1 (MQ=255)
                                                                                                                                         aTGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGCCGATGGCCAGGGCCAGTGCCTGCGATGTGGTGGCCCGGTACAGGTAGCCTGCGTCGATCAACGCCTGGTTGTTGAc  >  7:216501/1‑140 (MQ=255)
                                                                                                                                         aTGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGCCGATGGCCAGGGCCAGTGCCTGCGATGTGGTGGCCCGGTACAGGTAGCCTGCGTCGATCAACGCCTGGTTGTTGAc  <  8:216501/140‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                  
AACCGAGGTCGCCGGTGTGAATGCGGATGGCGTAGAACGGCCCCTTGGTAAGGGGTGCCAAGCATGGGTTGGGCGTATGTTGTGGGTCGCCCATGTAGCGGTTGTAACTGTTGCCGCCTTTGCCAAAGGCGTGGTCGATGCCATTGCCGGCATCGCTGTTGAACTGTGAAAGTGTTTGCGCAAACACCTGTGGATCGACGCCGATGGCCAGGGCCAGTGCCTGCGATGTGGTGGCCCGGTACAGGTAGCCTGCGTCGATCAACGCCTGGTTGTTGAC  >  NC_002947/2980870‑2981146

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