Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,257,472 T→A 10.3% T161S (ACC→TCC) ‡ PP_2854 ← hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,257,4720TA10.3% 157.6 / 8.6 58T161S (ACC→TCC) ‡PP_2854hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (26/26);  new base A (3/3);  total (29/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 = 3.17e-01

CAAGCAGCGTCGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCGGTGGTG  >  NC_002947/3257325‑3257611
                                                                                                                                                    |                                                                                                                                           
cAAGCAGCGTCGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTc                                                                                                                                            >  7:39636/1‑149 (MQ=255)
  agcagcGTCGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTg                                                                                                                                          >  1:194137/1‑149 (MQ=255)
    cagcGTCGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTgg                                                                                                                                         <  1:80138/149‑1 (MQ=255)
    cagcGTCGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGc                                                                                                                                        <  1:289944/149‑1 (MQ=255)
    cagcGTCGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGc                                                                                                                                        <  3:32197/149‑1 (MQ=255)
         tCGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCcggcg                                                                                                                                   >  1:160601/1‑149 (MQ=255)
         tCGAGGTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCcggcg                                                                                                                                   <  7:170801/149‑1 (MQ=255)
          cGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCcggcgg                                                                                                                                  >  3:192710/1‑149 (MQ=255)
                 tGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGcctcgt                                                                                                                            >  6:203724/1‑144 (MQ=255)
                 tGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTAc                                                                                                                           >  4:125319/1‑149 (MQ=255)
                 tGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTAc                                                                                                                           >  2:1908/1‑149 (MQ=255)
                    agcagGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGTACTGGCCGCCGGCCAGACCCTg                                                                                                                        >  7:94045/1‑149 (MQ=255)
                      cagGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGc                                                                                                                      >  6:37239/1‑149 (MQ=255)
                        ggCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGTACTGGCCGCCGGCCAGACCCTggtcg                                                                                                                    <  8:94045/149‑4 (MQ=255)
                        ggCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGTACTGGCCGCCGGCCAGACCCTggtcg                                                                                                                    <  4:39527/149‑4 (MQ=255)
                         gCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACTAGGGTCTGGCCggcggc                                                                                                                                 <  3:21967/135‑1 (MQ=255)
                         gCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCggcggc                                                                                                                                 >  4:21967/1‑135 (MQ=255)
                          ccACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTg                                                                                                                  <  8:137754/149‑1 (MQ=255)
                           cACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGc                                                                                                                 >  1:156261/1‑149 (MQ=255)
                             ccTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCtg                                                                                                               <  2:156261/149‑1 (MQ=255)
                               ttGGGCCAGGAGTTGCACTCTGCTGCCCTGGTTGTTGGTGGCATTGATGTTGAGCTGCACACCCAGACCGCTGGATGACTCACTGAGCGCGCCGGCGCCATTTTCCGCGTCCTGTGTCTGGGCGGCTGCCAGTACCTGGGCTTGCtgtt                                                                                                             <  8:212803/149‑1 (MQ=255)
                                      aGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTAcggc                                                                                                      <  1:1908/149‑1 (MQ=255)
                                         cGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCggg                                                                                                   <  4:192710/149‑1 (MQ=255)
                                           ttGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGTACTGGCCGCCGGCCAGACCCTggtcg                                                                                                                    <  3:203460/130‑4 (MQ=39)
                                           ttGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGTACTGGCCGCCGGCCAGACCCTggtcg                                                                                                                    >  4:203460/1‑127 (MQ=39)
                                               aCGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTg                                                                                             >  4:202818/1‑149 (MQ=255)
                                                 gctgctGTCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTg                                                                                                                        >  8:227672/1‑120 (MQ=255)
                                                 gctgctGTCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTg                                                                                                                        <  7:227672/120‑1 (MQ=255)
                                                        cccTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCttt                                                                                    <  2:99270/149‑1 (MQ=255)
                                                         ccTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTg                                                                                   <  2:194137/149‑1 (MQ=255)
                                                         ccTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTg                                                                                   <  8:189583/149‑1 (MQ=255)
                                                         ccTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTg                                                                                   <  5:37239/149‑1 (MQ=255)
                                                          cTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCtgtt                                                                                                             <  4:145355/122‑1 (MQ=255)
                                                          cTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCtgtt                                                                                                             >  3:145355/1‑122 (MQ=255)
                                                                 ttgCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGt                                                                            <  8:39636/148‑1 (MQ=255)
                                                                   gCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGt                                                                                                  <  2:72555/124‑1 (MQ=255)
                                                                   gCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGt                                                                                                  >  1:72555/1‑124 (MQ=255)
                                                                     tGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGata                                                                       >  7:112206/1‑149 (MQ=255)
                                                                         aTTGATGTTGAGCTTCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCGGGTTGGCTTTGGTCACGTTGATATCg                                                                    >  6:254812/1‑148 (MQ=255)
                                                                         aTTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGtt                                                                                          <  6:185218/126‑1 (MQ=255)
                                                                         aTTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGtt                                                                                          >  5:185218/1‑126 (MQ=255)
                                                                                ttGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACgttgt                                                             <  1:240664/148‑1 (MQ=255)
                                                                                   aGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTtgttgt                                                          <  8:112206/148‑1 (MQ=255)
                                                                                       gCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGTGTGCCTGGTTGGCTTTGGTCACGTTGATATCGCCGTGGTTGTAAgcg                                                     >  8:121317/1‑149 (MQ=255)
                                                                                        cacacccacaccGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCcggcgg                                                                                                                                  <  1:97631/71‑1 (MQ=255)
                                                                                        cacacccacaccGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCcggcgg                                                                                                                                  >  2:97631/1‑71 (MQ=255)
                                                                                               acaccGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGTCCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCgttgtt                                             >  1:229601/1‑149 (MQ=255)
                                                                                                 accGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGACAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGttgttg                                            >  5:295545/1‑148 (MQ=255)
                                                                                                               acacTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACgttgtt                                                            >  6:289931/1‑119 (MQ=255)
                                                                                                               acacTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACgttgtt                                                            <  5:289931/119‑1 (MQ=255)
                                                                                                                         cgcCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTgg                    <  5:254812/149‑1 (MQ=255)
                                                                                                                          gcCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTggg                   <  3:125319/149‑1 (MQ=255)
                                                                                                                                         cccGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCggtgg    >  1:301770/1‑149 (MQ=255)
                                                                                                                                          ccGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCggtggt   >  3:51805/1‑149 (MQ=255)
                                                                                                                                           cGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAgcgg                                                    <  6:333212/99‑1 (MQ=255)
                                                                                                                                           cGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAgcgg                                                    >  5:333212/1‑99 (MQ=255)
                                                                                                                                           cGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGACTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCggtggtg  >  5:89084/1‑149 (MQ=255)
                                                                                                                                           aggccAGGTACTGGCCGCCGGCCAGACCCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCGCCGGCTGCACGGACCACCTGGGTTACACCTTCggtggtg  >  7:253284/4‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                           
CAAGCAGCGTCGAGTTCTGCAGCAGGCCACCTTGGGCCAGGCGTTGCACGCTGCTGCCCTGGTTGTTGCTGGCATTGATGTTGAGCTGCACACCCACACCGCTGGAAGACACACTGAGCGCGCCGGCGCCATTTTCCCCGACCAGGGTCTGGCCGGCGGCCAGTACCTGGCCTTGCTGTTGTACGGCGGGTGCCTGGTTGGCTTTGGTCACGTTGATATCGACGTTGTTGTAAGCGGCGTTGTTGTCACCGGCTGCACGGACCACCTGGGTTACACCTTCGGTGGTG  >  NC_002947/3257325‑3257611

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: