Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,446,663 A→G 8.1% intergenic (+13/‑24) 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,6630AG 8.1% 131.0 / 2.9 49intergenic (+13/‑24)gpD/PP_3066tail formation/lytic enzyme
Reads supporting (aligned to +/- strand):  ref base A (24/21);  new base G (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 = 5.42e-01

GACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATGGGGAAGTTCGCCATCATCACGCGCCTGCGTATGGC  >  NC_002947/3446516‑3446811
                                                                                                                                                   |                                                                                                                                                     
gACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTa                                                                                                                                                       <  1:150243/148‑1 (MQ=255)
     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                                                                                                                                                     <  3:119011/132‑1 (MQ=255)
                  tttACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTaca                                                                                                                                                     >  4:119011/1‑132 (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                                                           <  8:286969/149‑1 (MQ=255)
                                                                                          aagaaaaagaaaGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGt                                                           >  2:328189/1‑149 (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                                                                 <  2:122165/126‑1 (MQ=255)
                                                                                                           cAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGt                                                                 >  1:122165/1‑126 (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                         <  4:178783/149‑1 (MQ=255)
                                                                                                                            agcagcTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATGGGGAAGTTCGcc                         <  8:260969/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)
                                                                                                                                                   |                                                                                                                                                     
GACGGGGTCGAGCAGGTATTTACCCAGTCCGGCTGGACTACGACCGTCGAGTGCAACGGCGGCAAGAAGGGCAAGGCCAAGGCCTCGGGCAAGAAAAAGAAAGACACCAAGCCGCTCAAGGTTGAGCAGCTGTAACCCCACGGCCGTACACGGCCATCACTGGAGAAACCAATGGCTATCTCAGTTCAACAGCTGCAAC‑AGATCCTCCCCAACGCCGGCCGCAAAGCCGGCGTTTATGTTCCCGGACTCAACGCTACTATGGGGAAGTTCGCCATCATCACGCGCCTGCGTATGGC  >  NC_002947/3446516‑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: ‑