Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,996,246 C→A 8.2% S94Y (TCC→TAC)  PP_3523 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,996,2460CA 8.2% 151.6 / 3.7 49S94Y (TCC→TAC) PP_3523hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (18/27);  new base A (2/2);  total (20/29)
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.78e-01

GACACCGTCGAGCACATTGAAACCACGCTGCTTGGCATGATCACACTCAATGGCACCACCCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAGGCCGCAAACCGTCATTTTTCCTGCAACACC  >  NC_002947/3996107‑3996385
                                                                                                                                           |                                                                                                                                           
gACACCGTCGAGCACATTGAAACCACGCTGCTTGGCATGATCACACTCAATGGCACCACCCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTg                                                                                                                                      >  4:175302/1‑147 (MQ=255)
 acacCGTCGAGCACATTGAAACCACGCTGCTTGGCATGATCACACTCAATGGCACCACCCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAg                                                                                                                                   >  7:168474/1‑149 (MQ=255)
 acacCGTCGAGCACATAGAAACCACGCTGCTTGGCATGATCACACTCAATGGCACCACCCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGaa                                                                                                                                    <  2:164041/148‑1 (MQ=255)
           gCACATTGAAACCACGCTGCTTGGCATGATCACACTCAATGGCACCACCCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCt                                                                                                                         >  2:311234/1‑149 (MQ=255)
            cacaTTGAACCCACGCTGCTTGGCATGATCACACTCAATGGCACCACCCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTACACCGt                                                                                                                                       <  7:220048/134‑1 (MQ=255)
            cacaTTGAAACCACGCTGCTTGGCATGATCACACTCAATGGCACCACCCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGGACTGTACACCGt                                                                                                                                       >  8:220048/1‑134 (MQ=255)
            cacaTTGAAACCACGCTGCTTGGCATGATCACACTCAATGGCACCAACCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGAGCTAAAGGAACTGCGCGTACTGTACACAGt                                                                                                                                       >  8:215758/1‑134 (MQ=255)
            cacaTTGAAACCACGCTGCTTGGCATCATCACACTCAATGGCACCAACCTCTCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTACACCga                                                                                                                                       <  7:215758/134‑2 (MQ=255)
             acaTTGAAACCACGCTGCTTGGCATGATCACACTCACTGGCACCACCCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGc                                                                                                                       <  5:13871/149‑1 (MQ=255)
               aTTGAAACCACGCTGCTTGGCATGATCACACTCAATGGCACCACCCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGc                                                                                                                          >  8:148233/1‑144 (MQ=255)
               aTTGAAACCACGCTGCTTGGCATGATCACACTCAATGGCACCACCCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGc                                                                                                                          <  7:148233/144‑1 (MQ=255)
                ttGAAACCACGCTGCTTGGCATGATCACACTCAATGGCACCACCCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCt                                                                                                                     >  4:44752/1‑148 (MQ=255)
                     aCCACGCTGCTTGGCATGATCACACTCAATGGCACCACCCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGAtggt                                                                                                               <  5:86354/149‑1 (MQ=255)
                          gctgctTGGCATGATCACACTCAATGGCACCACCCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACg                                                                                                          <  3:175302/149‑1 (MQ=255)
                              ctTGGCATGATCACACTCAATGGCACCACCCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCt                                                                                                      <  3:44752/149‑1 (MQ=255)
                              ctTGGCATGATCACACTCAATGGCACCACCCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCt                                                                                                      <  8:168474/149‑1 (MQ=255)
                               ttGGCATGATCACACTCAATGGCACCACCCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTg                                                                                                     >  1:324000/1‑149 (MQ=255)
                                  gCATGATCACACTCAATGGCACCACCCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGt                                                                                                  >  7:39367/1‑149 (MQ=255)
                                       aTCACACTCAATGGCACCACCCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTg                                                                                              <  1:67420/148‑1 (MQ=255)
                                        tCACACTCAATGGCACCACCCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGcc                                                                                            <  2:324000/149‑1 (MQ=255)
                                          acacTCAATGGCACCACCCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCAt                                                                                          <  6:241656/149‑1 (MQ=255)
                                            acTCAATGGCACCACCCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTgg                                                                                                    >  6:305455/1‑137 (MQ=255)
                                            acTCAATGGCACCACCCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTgg                                                                                                    <  5:305455/137‑1 (MQ=255)
                                              tCAATGGCACCACCCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGt                                                                                      >  2:279899/1‑149 (MQ=255)
                                                 aTGGCACCACCCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTc                                                                                   >  2:140221/1‑149 (MQ=255)
                                                       ccaccCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATc                                                                             <  1:311234/149‑1 (MQ=255)
                                                           ccTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGc                                                                         <  1:279899/149‑1 (MQ=255)
                                                               gCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCa                                                                      >  5:100392/1‑148 (MQ=255)
                                                                     cGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATc                                                                             <  5:303611/135‑1 (MQ=255)
                                                                     cGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATc                                                                             >  6:303611/1‑135 (MQ=255)
                                                                      gCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGAccc                                                              >  8:202952/1‑149 (MQ=255)
                                                                         aGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGc                                                                                                                       <  2:30901/89‑1 (MQ=255)
                                                                                  cTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTcg                                                  >  8:119699/1‑148 (MQ=255)
                                                                                   tGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCt                                                                                               <  2:327607/103‑1 (MQ=255)
                                                                                   tGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGc                                                                                             <  6:2785/105‑1 (MQ=255)
                                                                                     aaaaGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGgc                                               >  2:299396/1‑149 (MQ=255)
                                                                                      aaaGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGgc                                               <  7:202952/148‑1 (MQ=255)
                                                                                            ggTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCcg                                        >  3:315389/1‑149 (MQ=255)
                                                                                             gTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCcgc                                       >  8:81672/1‑149 (MQ=255)
                                                                                                cgcgATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCgcgc                                     <  1:140221/148‑1 (MQ=255)
                                                                                                   gATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGa                                 <  7:81672/149‑1 (MQ=255)
                                                                                                   gATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGGTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGa                                 <  5:241006/149‑1 (MQ=255)
                                                                                                    aTTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAg                                <  6:100392/149‑1 (MQ=255)
                                                                                                    aTTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAg                                <  1:299396/149‑1 (MQ=255)
                                                                                                    aTTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAg                                <  8:39367/149‑1 (MQ=255)
                                                                                                                        aaGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCa                                                  <  8:30601/111‑1 (MQ=255)
                                                                                                                          ggAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAGGCCGCAAACCGTCATTTTTCCt          <  4:315389/149‑1 (MQ=255)
                                                                                                                                     tACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAGGCCGCAAACCGTCATTTTTCCTGCAACAcc  >  2:206391/1‑146 (MQ=255)
                                                                                                                                     tACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAGGCCGCAAACCGTCATTTTTCCTGCAACAcc  <  1:206391/146‑1 (MQ=255)
                                                                                                                                           |                                                                                                                                           
GACACCGTCGAGCACATTGAAACCACGCTGCTTGGCATGATCACACTCAATGGCACCACCCTCGCTGACCGCAAGCGGGAGGCTGAAAAGCTGGTCCGCGATTACATCCATGACGCGCTCAAGGAACTGCGCGTACTGTCCACCGTGAAGTTCTACGGCTGCCTGATGGTGGACGGGCTGGGTCCTGCCATGCGTTTCGACATCCTGCCCAAATGACCCTCGACGGCCTCAGGCGCCCCCGCGCCCGAGGCCGCAAACCGTCATTTTTCCTGCAACACC  >  NC_002947/3996107‑3996385

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