Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,160,643 C→T 10.2% G138S (GGC→AGC)  PP_4543 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,160,6430CT10.2% 133.2 / 5.3 49G138S (GGC→AGC) PP_4543hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (19/25);  new base T (2/3);  total (21/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

ACCCAGGCAGGCCCGCCAACTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCG  >  NC_002947/5160496‑5160780
                                                                                                                                                   |                                                                                                                                         
aCCCAGGCAGGCCCGCCAACTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCg                                                                                                                                          >  5:183068/1‑149 (MQ=255)
  ccAGGCAGGACCGCCAACTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTGGCCCAGATACTCCGGCAGGGCCTGCCGCATCTACAAACCGGCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCt                                                                                                                                        >  5:267655/1‑149 (MQ=255)
        aggcCCGCCAACTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCgg                                                                                                                                  <  1:352214/149‑1 (MQ=255)
             cGCCAACTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTg                                                                                                                              >  4:222195/1‑148 (MQ=255)
               ccAACTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAAc                                                                                                                           <  7:339782/149‑1 (MQ=255)
                 aaCTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCTGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCGGCCCACGGCTGCTGAcgg                                                                                                                                   <  8:250848/139‑3 (MQ=255)
                 aaCTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCGGCCCACGGCTGCTGAcgg                                                                                                                                   >  7:250848/1‑137 (MQ=255)
                 aaCTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACgg                                                                                                                         <  4:170616/149‑1 (MQ=255)
                   cTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGca                                                                                                                       <  5:294497/149‑1 (MQ=255)
                      cgcgTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGcaca                                                                                                                     <  4:226002/148‑1 (MQ=255)
                           ccAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGcc                                                                                                                >  3:310791/1‑148 (MQ=255)
                             aGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAg                                                                                                             >  1:318450/1‑149 (MQ=255)
                                  cGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACgg                                                                                                                         >  1:136488/1‑132 (MQ=255)
                                  cGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACgg                                                                                                                         <  2:136488/132‑1 (MQ=255)
                                    ggTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGtttt                                                                                                      <  4:310791/149‑1 (MQ=255)
                                      tCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTgg                                                                                                    <  6:6361/149‑1 (MQ=255)
                                      tCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTgg                                                                                                    <  2:347412/149‑1 (MQ=255)
                                       cGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTggg                                                                                                   <  5:89247/149‑1 (MQ=255)
                                         gTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCGGCCCACGGCTGCTGAcgg                                                                                                                                   <  6:103390/115‑3 (MQ=255)
                                         gTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCGGCCCACGGCTGCTGAcgg                                                                                                                                   >  5:103390/1‑113 (MQ=255)
                                            aaGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGgtagt                                                                                               <  6:132908/148‑1 (MQ=255)
                                                  aTGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGcac                                                                                                                      >  6:313949/1‑119 (MQ=255)
                                                   tgtgCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATa                                                                                       <  1:226824/149‑1 (MQ=255)
                                                      gcgcTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCg                                                                                    <  5:90584/149‑1 (MQ=255)
                                                        gcgccGGAAAGCCCACCCTTGCCCAGATCCTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCCGCCGTGGGCCGCTGACCGGGCTGTACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGc                                                                                  <  6:267655/146‑1 (MQ=255)
                                                              gAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAAGCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTc                                                                            <  2:8255/149‑1 (MQ=255)
                                                              gAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCt                                                                             >  2:169238/1‑148 (MQ=255)
                                                              gAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCt                                                                             >  1:49953/1‑148 (MQ=255)
                                                                  gCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACaa                                                                        <  1:169238/149‑1 (MQ=255)
                                                                  gCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACaa                                                                        <  3:222195/149‑1 (MQ=255)
                                                                    ccACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACaa                                                                        >  4:51947/1‑147 (MQ=255)
                                                                    ccACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACaa                                                                        <  3:51947/147‑1 (MQ=255)
                                                                              ccAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCaa                                                            >  7:337886/1‑149 (MQ=255)
                                                                                              ggCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACAt                                            <  6:183068/149‑1 (MQ=255)
                                                                                                ccTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCg                                                                  <  3:226842/125‑1 (MQ=255)
                                                                                                ccTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCg                                                                  >  4:226842/1‑125 (MQ=255)
                                                                                                ccTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGaa                                                 <  6:153090/142‑1 (MQ=255)
                                                                                                ccTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGaa                                                 >  5:153090/1‑142 (MQ=255)
                                                                                                  tGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTgg                                                                                                    >  7:251809/1‑89 (MQ=255)
                                                                                                  tGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTgg                                                                                                    <  8:251809/89‑1 (MQ=255)
                                                                                                    cagcaTCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGc                                      <  2:127404/149‑1 (MQ=255)
                                                                                                         tCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAAc                                 <  2:318450/149‑1 (MQ=255)
                                                                                                                    tCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGGCACtt                      >  1:80316/1‑149 (MQ=255)
                                                                                                                                gCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCggg          >  7:287395/1‑149 (MQ=255)
                                                                                                                                 cggTCAGCGGCCCACGGCTGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGc                                                                <  4:366432/92‑1 (MQ=38)
                                                                                                                                 ccGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCggg          >  8:34504/1‑148 (MQ=255)
                                                                                                                                    tCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCCGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCaa       >  3:239717/1‑148 (MQ=255)
                                                                                                                                      agcagcCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACa     <  2:49953/148‑1 (MQ=255)
                                                                                                                                        cagcCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAgcg  >  4:173312/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                         
ACCCAGGCAGGCCCGCCAACTGCGCGTCCAGCGCCGGGTCGGTGAAGTTGATGTGCGCTCCGGAAAGCCCACCCTTGCCCAGATACTCCGGCAGGGCCTGCAGCATCAACAAACCGTCAGCCGGATCGGCCGTCAGCAGCCGTGGGCCGCTGACCGGGCTGAACGGCACAGCGCCAAGCAGTTTTGGGTAGTAGGCGATACCGGCCCGCTCACAAGCATCGGCCCAACCATGATCGAACACATACTCGCCAAACGAGTGCCACTTGCGGTACGCGGGCAACAGCG  >  NC_002947/5160496‑5160780

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: