Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,389,508 T→A 10.1% T411T (ACA→ACT)  PP_3865 ← tail protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,389,5080TA10.1% 144.5 / 3.3 54T411T (ACA→ACT) PP_3865tail protein
Reads supporting (aligned to +/- strand):  ref base T (22/26);  new base A (4/2);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 4.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGC‑CGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCG  >  NC_002947/4389364‑4389652
                                                                                                                                                 |                                                                                                                                                
cTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGC‑CGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTgg                                                                                                                                              >  5:159332/1‑149 (MQ=255)
    aTCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGC‑CGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCcgc                                                                                                                                          <  8:84406/149‑1 (MQ=255)
       aGCTCCCACATGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGC‑CGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGCGGCCACAGTGGCcacgc                                                                                                                                        <  7:50458/148‑2 (MQ=255)
       aGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGC‑CGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTgg                                                                                                                                              >  8:215065/1‑142 (MQ=255)
       aGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGC‑CGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTgg                                                                                                                                              <  7:215065/142‑1 (MQ=255)
       aGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGC‑CGCCCCCGACCCCCCGGCGCACACGCCAAAACGCGATGCCCCGCAGGCCCCACCCAGCCCCCCCCAGGCCCCACTGGCCGCGa                                                                                                                                        >  4:308415/1‑148 (MQ=255)
                           gaaCGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGC‑CGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCgg                                                                                                                   <  2:17535/149‑1 (MQ=255)
                           gaaCGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGC‑CGCCGCGATCCCCGCGGCGATAAGGCCAATACGCGTGGCCGCGAAGGCCACCCCGAGCACCCGCCTGGCCACTGTGGCCGCGAAAACCCCCGCACGCACCGCgg                                                                                                                   >  1:945/1‑149 (MQ=255)
                             aCGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGC‑CGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTg                                                                                                                 <  4:4479/149‑1 (MQ=255)
                                  ccGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGC‑CGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCTCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGgc                                                                                                            <  1:59583/149‑1 (MQ=255)
                                     aCCTTCTCCCAGTTGGCGACGATCAGACCGGCCGC‑CGCCGCGATACCCCCCGGGAAGAGGCCAAAAGGCGTTGCCCCGAACGCCACACCCAGCCCCCCCCTGGCCACAGTGGCCGCGCAAACCGCCCCACCCCGCCCGGTGAAAgccgc                                                                                                         >  2:8112/1‑149 (MQ=255)
                                       cTTCTCCCAGTTGGCGACGATCAGACCGGCCGC‑CGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGccc                                                                                                       <  6:159332/149‑1 (MQ=255)
                                              cAGTTGGCGACGATCAGACCGGCCGC‑CGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAAt                                                                                                <  6:186966/149‑1 (MQ=255)
                                                      gacgaTCAGACCGGCCGCGCGCCCCCAACCCCCCGGCCAAAAGGCCAAAAAGCGTGGCCGCGAAGACCACACCCACCCCCCCCCAGGCAACAGTGGCCCCGAACACCCCCCCCCCCACCCCGGGGAAGGCCCCCCCGCCaaaaaca                                                                                            >  2:287473/1‑141 (MQ=255)
                                                        cgaTCAGACCGGCCGC‑CGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCcg                                                                                      <  1:300076/149‑1 (MQ=255)
                                                         gaTCAGACCGGCCGC‑CGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCcgc                                                                                     >  1:53260/1‑149 (MQ=255)
                                                             aGACCGGCCGC‑CGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGa                                                                                                                <  3:197293/118‑1 (MQ=255)
                                                             aGACCGGCCGC‑CGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCa                                                                                  <  2:191094/148‑1 (MQ=255)
                                                                 cGGCCGC‑CGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCCCACCGCGCCCCCGCGTGGCCCCTGTGGCCGCGAACACCGCCGCACGCCGCGCGGTGAAGGCCGCCCCGGCAATGACGCTGCCCCGCACCCAc                                                                             >  7:213741/1‑149 (MQ=255)
                                                                      gc‑cgccgcGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCTcc                                                                                >  2:111371/1‑141 (MQ=255)
                                                                          gccgcGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGtt                                                                     >  2:78644/1‑149 (MQ=255)
                                                                             gcgATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTc                                                                  <  3:308415/149‑1 (MQ=255)
                                                                             gcgATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTc                                                                  <  2:53260/149‑1 (MQ=255)
                                                                                 tCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTg                                                              <  2:945/149‑1 (MQ=255)
                                                                                  ccccGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGa                                                             >  1:291859/1‑149 (MQ=255)
                                                                                     cgcgGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGatca                                                          >  3:125916/1‑149 (MQ=255)
                                                                                             tCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTgcgc                                                   >  6:207566/1‑148 (MQ=255)
                                                                                                ggCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCata                                               >  8:1153/1‑149 (MQ=255)
                                                                                                ggCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCata                                               >  8:29140/1‑149 (MQ=255)
                                                                                                 gCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCatat                                              >  4:101375/1‑149 (MQ=255)
                                                                                                 gCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCatat                                              >  6:232159/1‑149 (MQ=255)
                                                                                                   cAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATgg                                            >  5:155476/1‑149 (MQ=255)
                                                                                                    aaTCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCCCGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCTTATGGc                                           >  5:189738/1‑149 (MQ=255)
                                                                                                    aaTCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGc                                           <  1:287473/149‑1 (MQ=255)
                                                                                                     aTCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGtt                                                                     >  3:96816/1‑122 (MQ=255)
                                                                                                     aTCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGtt                                                                     <  4:96816/122‑1 (MQ=255)
                                                                                                     aTCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCt                                          <  7:29140/149‑1 (MQ=255)
                                                                                                             ggCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACg                                  >  8:268461/1‑149 (MQ=255)
                                                                                                              gCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACgg                                 <  6:155476/149‑1 (MQ=255)
                                                                                                               ccGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGc                                <  2:185323/149‑1 (MQ=255)
                                                                                                                          aaacCGAGGACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTc                     <  6:189738/147‑1 (MQ=255)
                                                                                                                            acCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGa                   <  1:78644/149‑1 (MQ=255)
                                                                                                                                 gCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGAt                                       <  2:319300/124‑1 (MQ=255)
                                                                                                                                 gCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGAt                                       >  1:319300/1‑124 (MQ=255)
                                                                                                                                  cACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACAAACTGCGGGATGTCTCTGATCATCTGCGCCATATGGCTGATCAACGGCCGCAGTTGCTCTAGCACCg             >  2:149130/1‑149 (MQ=255)
                                                                                                                                      ggtgTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCa           <  8:177474/144‑1 (MQ=255)
                                                                                                                                      cgcgTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCAtt         <  7:1153/149‑1 (MQ=255)
                                                                                                                                         gcggCCACAGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTg        >  4:236166/3‑147 (MQ=255)
                                                                                                                                         gcggCCACAGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTg        <  3:236166/145‑1 (MQ=255)
                                                                                                                                             ccACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCg  <  2:232210/149‑1 (MQ=255)
                                                                                                                                             ccACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCg  >  1:232210/1‑149 (MQ=255)
                                                                                                                                             ccACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCg  <  3:101375/149‑1 (MQ=255)
                                                                                                                                              cACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTc                     >  1:4584/1‑129 (MQ=255)
                                                                                                                                              cACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTc                     <  2:4584/129‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                
CTTGATCAGCTCCCACAGGGCGCTGAAGAACGGCCCGACCTTCTCCCAGTTGGCGACGATCAGACCGGCCGC‑CGCCGCGATCCCCGCGGCGATCAGGCCAATCGGCGTGGCCGCGAAGGCCACACCGAGCACCCGCGTGGCCACTGTGGCCGCGAACACCGCCGCACGCAGCGCGGTGAAGGCCGCCCCGGCAATGACGATGCCGCGCACCAACTGCGGGTTGTCGCTGATCATCTGCGCCATATGGCTGATCCACGGCCGCAGTTGCTCGAGCACCGCATTGATGCCG  >  NC_002947/4389364‑4389652

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