Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,440,620 G→T 8.3% K121N (AAG→AAT)  PP_3060 → tail sheath protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,440,6200GT 8.3% 150.9 / 4.0 48K121N (AAG→AAT) PP_3060tail sheath protein
Reads supporting (aligned to +/- strand):  ref base G (17/27);  new base T (2/2);  total (19/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 = 5.80e-01

GTTCTTCGGCCTGGCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAA  >  NC_002947/3440473‑3440767
                                                                                                                                                    |                                                                                                                                                   
gTTCTTCGGCCTGGCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAgg                                                                                                                                                    <  3:246694/149‑1 (MQ=255)
    ttcGGCCTGGCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAATGCCgt                                                                                                                                                >  2:238406/1‑148 (MQ=255)
      cGGCCTGGCTGACGGCAACGCCTTCAACGGCGCATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAg                                                                                                                                               <  5:303075/148‑1 (MQ=255)
       ggCCTGGCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGc                                                                                                                                              <  4:189905/148‑1 (MQ=255)
             gCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAg                                                                                                                                         <  2:36557/147‑1 (MQ=255)
              ctgtccgCAACGCCTTCAACGGCGTATTCCGTGTCTGGTTCAACATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGGCGAGCTGAAGCAc                                                                                                                                      <  3:29888/143‑1 (MQ=255)
                   ggCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCAt                                                                                                                                 >  8:157926/1‑149 (MQ=255)
                   ggCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCAt                                                                                                                                 >  4:115535/1‑149 (MQ=255)
                    gCAACGCCTTCAACGGAGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAATCCATTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCAtt                                                                                                                                <  8:255521/149‑1 (MQ=255)
                     cAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTg                                                                                                                               >  7:288969/1‑149 (MQ=255)
                      aaCGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGTGCTTCAGCTCGGCCTTgt                                                                                                                              <  6:241421/149‑2 (MQ=255)
                       aCGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGACGCTGGACATGGGCGACTGGAAGGCCGGCGACAATGCCGTGCTTCAGCTCGGCCTTgt                                                                                                                              <  1:238406/149‑2 (MQ=255)
                       aCGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGcc                                                                                                                             <  3:115535/149‑1 (MQ=255)
                         gCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGTGCTTCAGCTCGGCCTTGTCgc                                                                                                                           <  8:84528/149‑2 (MQ=255)
                          ccTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAATA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTc                                                                                                                          <  4:138114/149‑1 (MQ=255)
                           ctgcatcGGGGTATTGCGAGGCTCGTTCTAGTTCCAGAAGGGCGTATCGCTTGCGGTCGTCGTCACCCTGCGCGGCACCGTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCAGGGCATTGCCGTCa                                                                                                                         <  6:206368/143‑1 (MQ=255)
                           cTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCa                                                                                                                         <  6:146691/149‑1 (MQ=255)
                                      tATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAg                                                                                                              <  2:199056/149‑1 (MQ=255)
                                          ccGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTc                                                                                                           <  3:175283/148‑1 (MQ=255)
                                               gCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTc                                                                                                     >  6:212569/1‑149 (MQ=255)
                                                        aGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAg                                                                                            <  6:87461/149‑1 (MQ=255)
                                                           tCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTc                                                                                         <  2:187132/149‑1 (MQ=255)
                                                               gAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATc                                                                                      >  4:300107/1‑148 (MQ=255)
                                                                 aGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTAc                                                                                   <  7:247394/149‑1 (MQ=255)
                                                                        aaCCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGAtt                                                                             >  3:122189/1‑148 (MQ=255)
                                                                           ccTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGATCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGAcc                                                                         <  7:217515/149‑1 (MQ=255)
                                                                                   tcgtcgtcACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCgg                                                                  >  6:37101/1‑148 (MQ=255)
                                                                                     gtcgtcACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCa                                                                <  8:288969/148‑1 (MQ=255)
                                                                                     gtcgtcACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCAt                                                               <  7:157926/149‑1 (MQ=255)
                                                                                     gtcgtcACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCAt                                                               >  6:171445/1‑149 (MQ=255)
                                                                                        gtcACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGaa                                                                                                        <  5:275071/105‑1 (MQ=255)
                                                                                        gtcACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCAGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGaa                                                                                                        >  6:275071/1‑105 (MQ=255)
                                                                                             ccTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATc                                                                                      >  3:75101/1‑118 (MQ=255)
                                                                                             ccTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATc                                                                                      <  4:75101/118‑1 (MQ=255)
                                                                                             ccTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCgtg                                                       <  5:212569/149‑1 (MQ=255)
                                                                                              cTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCgtgt                                                      >  1:124012/1‑149 (MQ=255)
                                                                                              cTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCgtgt                                                      >  6:103307/1‑149 (MQ=255)
                                                                                                      cACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACg                                              <  5:37101/149‑1 (MQ=255)
                                                                                                          cTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCAc                                          <  4:122189/149‑1 (MQ=255)
                                                                                                          cTGAAAGA‑GCTGGACATGGGCGACTGGAACGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCAc                                          >  6:97819/1‑149 (MQ=255)
                                                                                                           tGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCAcc                                         >  6:146421/1‑149 (MQ=255)
                                                                                                                       gACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTg                                <  5:146421/147‑1 (MQ=255)
                                                                                                                              gCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAg                       >  8:26071/1‑149 (MQ=255)
                                                                                                                                           cGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCg          <  2:124012/149‑1 (MQ=255)
                                                                                                                                            gacggcAATGCCGTGCTTCAGCTCGGCCTTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGtc                                                     <  7:31539/100‑1 (MQ=37)
                                                                                                                                            gacggcAATGCCGTGCTTCAGCTCGGCCTTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGtc                                                     >  8:31539/6‑105 (MQ=37)
                                                                                                                                             gCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGc        <  3:300107/149‑1 (MQ=255)
                                                                                                                                                  aaGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGt    >  2:234600/1‑148 (MQ=255)
                                                                                                                                                   cggCCGAGCCGACGCCCGGCATGGCCCTCACCTCCTCCCACCCCCCAGTCGGCGGCGAGGGCATCTCCGAGATCGACCCGGTCGGCATGAAGCGTGGCAGCCCCGCCACCGCCCCACTCGCCAGTCAGCGCGCCGACCTCGGCCTgcac  <  4:84151/148‑4 (MQ=255)
                                                                                                                                                   aGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTaa  >  2:170347/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
GTTCTTCGGCCTGGCTGACGGCAACGCCTTCAACGGCGTATTCCGAGGCTCGTTCAAGATCCAGAAGGGCGAAACCCTTGCGGTCGTCGTCACCCTGCGCGGCACCCTGAAAGA‑GCTGGACATGGGCGACTGGAAGGCCGGCGACAAGGCCGAGCTGAAGCACGGCATTGCCGTCACCTACTACAAGCTCGAAGTCGGCGGCGAGGTCATCTACGAGATTGACCCGGTCGGCATGAAGCGTGTCATCAACGGCACCGACCAACTGGCCAGTCAGCGCGCCGACCTCGGCCTGTAA  >  NC_002947/3440473‑3440767

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