Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,446,664 C→T 8.1% intergenic (+14/‑23) gpD → / → PP_3066 tail formation/lytic enzyme

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,446,6640CT 8.1% 139.4 / 2.2 49intergenic (+14/‑23)gpD/PP_3066tail formation/lytic enzyme
Reads supporting (aligned to +/- strand):  ref base C (24/21);  new base T (2/2);  total (26/23)
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.50e-01

GGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATGGGGAAGTTCGCCATCATCACGCGCCTGCGTATGGC  >  NC_002947/3446521‑3446811
                                                                                                                                               |                                                                                                                                                    
ggTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGc                                                                                                                                                 >  7:286969/1‑149 (MQ=255)
ggTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGc                                                                                                                                                 >  6:234569/1‑149 (MQ=255)
     aGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTGTACGGccgtgg                                                                                                                                            <  8:15536/149‑5 (MQ=255)
          gTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGga                                                                                                                                       <  6:313221/149‑1 (MQ=255)
             tttACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTaca                                                                                                                                                     >  4:119011/1‑132 (MQ=255)
             tttACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTaca                                                                                                                                                     <  3:119011/132‑1 (MQ=255)
                  ccAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCaa                                                                                                                               >  4:24041/1‑149 (MQ=255)
                   cAGTCGGGCTGGACTTCGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAAt                                                                                                                              <  8:119206/149‑1 (MQ=255)
                           cTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTAtct                                                                                                                      <  5:234569/149‑1 (MQ=255)
                              gACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAg                                                                                                                   >  7:312534/1‑149 (MQ=255)
                              gACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAg                                                                                                                   <  6:72990/149‑1 (MQ=255)
                                cTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGtt                                                                                                                 <  8:77077/149‑1 (MQ=255)
                                   cGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACCGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCaa                                                                                                              >  4:30368/1‑149 (MQ=255)
                                       cGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCATGGCCGTGTACGGccgtggg                                                                                                                                           >  8:94174/1‑111 (MQ=255)
                                       cGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCATGGCCGTGTACGGccgtggg                                                                                                                                           <  7:94174/116‑6 (MQ=255)
                                       cGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGc                                                                                                          >  7:138748/1‑149 (MQ=255)
                                                cAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGc                                                                                                                           <  2:220753/123‑1 (MQ=255)
                                                cAACGGCGGCAAGAAGGGCAAGGCCAAGGCCACGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGc                                                                                                                           >  1:220753/1‑123 (MQ=255)
                                                  aCGGCGGCAAGAAGGGCAAGGCCAAGGCCTCTGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACGCCACGGCCGTACACGGCCATCACTGGAGAAAc                                                                                                                                  <  3:327287/114‑1 (MQ=255)
                                                  aCGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCAAGGCCGTACACGGCCATCACTGGAGAAAc                                                                                                                                  >  4:327287/1‑114 (MQ=255)
                                                         cAAGAAGGGCAAGGCCAAAGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTcccc                                                                                        >  7:260969/1‑148 (MQ=255)
                                                                  caaggccaaggcCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTGaa                                                                                                              <  3:321043/118‑1 (MQ=255)
                                                                  caaggccaaggcCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCaa                                                                                                              >  4:321043/1‑118 (MQ=255)
                                                                       ccaaggccTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCa                                                                           >  5:108501/1‑147 (MQ=255)
                                                                          aggccTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGcc                                                                      >  2:313927/1‑149 (MQ=255)
                                                                           ggccTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCg                                                                     >  3:178783/1‑149 (MQ=255)
                                                                             ccTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGc                                                                   <  8:312534/149‑1 (MQ=255)
                                                                              cTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCg                                                                  >  5:84039/1‑149 (MQ=255)
                                                                                 gggCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCTTCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGc                                                                       <  5:267511/141‑1 (MQ=255)
                                                                                 gggCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGc                                                                       >  6:267511/1‑141 (MQ=255)
                                                                                 gggCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGtt                                                                <  3:30368/148‑1 (MQ=255)
                                                                                     aagaaaaagaaaGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGt                                                           >  2:328189/1‑149 (MQ=255)
                                                                                     aagaaaaagaaaGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGt                                                           <  8:286969/149‑1 (MQ=255)
                                                                                          aaagaaagACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCg                                                      >  2:309497/1‑149 (MQ=255)
                                                                                             gaaagaCACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGAc                                                   <  1:309497/149‑1 (MQ=255)
                                                                                              aaagaCACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACt                                                  >  6:33906/1‑149 (MQ=255)
                                                                                                      cAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGt                                                                 >  1:122165/1‑126 (MQ=255)
                                                                                                      cAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGt                                                                 <  2:122165/126‑1 (MQ=255)
                                                                                                        aGCTGCTCAACCTTGAGCGGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGctac                                         >  6:315764/1‑148 (MQ=255)
                                                                                                           cGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATg                                     <  1:313927/149‑1 (MQ=255)
                                                                                                            gCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATgg                                    <  8:138748/149‑1 (MQ=255)
                                                                                                                aaGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATgggg                                  <  4:293008/147‑1 (MQ=255)
                                                                                                                      gAGCAGCTGTAACCCCACGGCCGTACACGGCAATCACTGGAGAAACCAATGGCTATCTCAGTTAAACAGCTGCACC‑AGATCCTCCCCAACGCCGGCCGCAACGCCGGCGTTTATGTTCCCGGCCTCCAGGCTAGAAGGGGGAAGTTCGc                          >  4:129068/1‑149 (MQ=255)
                                                                                                                       agcagcTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATGGGGAAGTTCGcc                         <  8:260969/149‑1 (MQ=255)
                                                                                                                       agcagcTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATGGGGAAGTTCGcc                         <  4:178783/149‑1 (MQ=255)
                                                                                                                       agcagcTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATGGGGAAGTTCGc                          >  3:347564/1‑148 (MQ=255)
                                                                                                                           gcTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATGGGGAAGTTCGCcatca                     <  4:347564/149‑1 (MQ=255)
                                                                                                                             tGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATgg                                    >  4:225732/1‑132 (MQ=255)
                                                                                                                             tGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATgg                                    <  3:225732/132‑1 (MQ=255)
                                                                                                                                 aCCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATGGGGAAGTTCGCCATCATCAcgc               >  1:182923/1‑149 (MQ=255)
                                                                                                                                   gtgatggCCGTGTACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAACGAAAACATCCCCAACCCCGGCCGCAAAGCCGGCGTTAATTTTCCCGGAAAAAACGATAGTATGGGGAAGTTCGCCATCATCAcgcc              >  1:287475/6‑148 (MQ=255)
                                                                                                                                              acacGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCAGGCGTTTATGTTCCCGGACTCAACGCTACTATGGGGAAGTTCGCCATCATCACGCGCCTGCGTATGGc  >  7:230246/1‑149 (MQ=255)
                                                                                                                                               cacGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATGGGGAAGTTCGCCATCATCACGCGCCTGCGTATGGc  <  6:108501/148‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                    
GGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATGGGGAAGTTCGCCATCATCACGCGCCTGCGTATGGC  >  NC_002947/3446521‑3446811

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