Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,206,302 G→C 9.7% D536H (GAC→CAC)  PP_2811 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,206,3020GC 9.7% 110.9 / 2.1 39D536H (GAC→CAC) PP_2811hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (16/19);  new base C (2/2);  total (18/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 = 9.00e-01

GCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCTGTTGTCACGCGTCGAGGGTGTGTTGTTCGGGCATCGTCGTCTGGTG  >  NC_002947/3206157‑3206441
                                                                                                                                                  |                                                                                                                                           
gCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGAcc                                                                                                                                          >  3:186592/1‑149 (MQ=255)
  cGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACcgc                                                                                                                                        <  6:278124/149‑1 (MQ=255)
          aGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTAt                                                                                                                                  <  5:330027/147‑1 (MQ=255)
            cgcCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCg                                                                                                                              >  2:200133/1‑149 (MQ=255)
            cgactccTTCACCGAAGTCGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCg                                                                                                                              <  7:51385/144‑1 (MQ=255)
                   ttCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATcgc                                                                                                                       >  2:276019/1‑149 (MQ=255)
                     cACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATcgc                                                                                                                       <  4:186592/147‑1 (MQ=255)
                     cACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGa                                                                                                                     >  5:38032/1‑149 (MQ=255)
                            gACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCg                                                                                                              <  2:262400/149‑1 (MQ=255)
                                     agcaTGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACAGGGATGATGGCGGCATCGCGATTTGCCGCACCTtggg                                                                                                     >  4:130926/1‑147 (MQ=255)
                                          gAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCTACAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGc                                                                                                 >  1:241382/1‑149 (MQ=255)
                                                    cTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTCACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAAc                                                                                      >  8:148242/1‑149 (MQ=255)
                                                    cTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTCACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAAc                                                                                      >  8:144041/1‑149 (MQ=255)
                                                      gAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGc                                                                                             <  3:136527/140‑1 (MQ=255)
                                                      gAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGc                                                                                             >  4:136527/1‑140 (MQ=255)
                                                        aTTTCGTACGAATCCCGGCCAGTTGGCAGATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGtt                                                                                  <  2:262375/149‑1 (MQ=255)
                                                         tttCGTACGAAACCCGGCCAGTTGGCATATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGtt                                                                                  <  7:55584/148‑1 (MQ=255)
                                                           tCGTACGAAACCCGGCCAGTTGGCAGATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCt                                                                                <  6:38032/148‑1 (MQ=255)
                                                                cGAAACCCGGCCAGTTGGCAGATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTCACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTaaa                                                                           <  7:148242/148‑1 (MQ=255)
                                                                cGAAACCCGGCCAGTTGGCAGATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTCACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTaaa                                                                           <  7:144041/148‑1 (MQ=255)
                                                                    aCCCGGCCAGTTGGCAGATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACATTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTa                                                                       >  1:154331/1‑148 (MQ=255)
                                                                     cccGGCCAGTTGGCAGATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATc                                                                     <  3:64689/149‑1 (MQ=255)
                                                                     cccGGCCAGTTGGCAGATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATc                                                                     >  6:289511/1‑149 (MQ=255)
                                                                          ccAGTTGGCAGATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAg                                                                <  2:241382/149‑1 (MQ=255)
                                                                             gtcggCAGATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGc                                                               <  1:276019/144‑1 (MQ=255)
                                                                                         cccT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGcc                                                 >  8:244895/1‑149 (MQ=255)
                                                                                                           ttCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATATCCCGCCCCTTGTGAGGCCTGCCGTGAACGCTTATGTAAAGCTATCTGAAGCCCCGGCACGTGCCCTGTTGTCAAGCGTCGAg                                >  6:262850/1‑149 (MQ=255)
                                                                                                            tCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATcg                                                                                                                        <  4:100059/60‑1 (MQ=255)
                                                                                                            tCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATcg                                                                                                                        >  3:100059/1‑60 (MQ=255)
                                                                                                             cATGGGCGTCTAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCTGTTGTCACGCGTCGAgg                               <  1:200133/148‑1 (MQ=255)
                                                                                                             cATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCTGTTGTCACGCGTCGAgg                               <  5:289511/148‑1 (MQ=255)
                                                                                                                gggCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCTGTTGTCACGCGTCGAGGgtgt                           >  5:97753/1‑149 (MQ=255)
                                                                                                                     gcgAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGccc                                                             >  4:33492/1‑110 (MQ=255)
                                                                                                                     gcgAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGccc                                                             <  3:33492/110‑1 (MQ=255)
                                                                                                                          aCCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCTGTTGTCAcgc                                      >  4:209843/1‑128 (MQ=255)
                                                                                                                          aCCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCTGTTGTCAcgc                                      <  3:209843/128‑1 (MQ=255)
                                                                                                                                tACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCTGTTGTCACGCGTCGAGGGTGTGTTGTTCGGGCAtcgt           >  8:239539/1‑149 (MQ=255)
                                                                                                                                        gCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACGTTGTGAGGCCTGGCGTGAACGGTGCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCTGTTGTCACGCGTCGAGGGTGTGTTGTTCGGGCATCGTCGTCTGGt   <  5:262850/149‑1 (MQ=255)
                                                                                                                                         cTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCTGTTGTCACGCGTCGAGGGTGTGTTGTTCGGGCATCGTCGTCTGGTg  <  7:244895/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                           
GCCGAACCTGAGCGCCACCTTCACCGAAGACGCGAGCAGCATGAACCTGTACCTGAATTTCGTACGAAACCCGGCCAGTTGGCAGATTTCCCT‑CAACTACGCGAAGTTCATGGGCGGCGAAACCCCCTACGACCAGCTTTATCGTGACCGCGATTATGTCGGCATCGCGATTTCCCGCACCTTGTGAGGCCTGCCGTGAACGGTTCTGTAAAGCTATCTGAAGCCCCGGCACGTGGCCTGTTGTCACGCGTCGAGGGTGTGTTGTTCGGGCATCGTCGTCTGGTG  >  NC_002947/3206157‑3206441

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