Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,404,284 C→A 7.7% P329T (CCG→ACG)  hflC → HflB protease modulator specific for phage lambda cII repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,404,2840CA 7.7% 168.7 / 3.3 52P329T (CCG→ACG) hflCHflB protease modulator specific for phage lambda cII repressor
Reads supporting (aligned to +/- strand):  ref base C (25/23);  new base A (2/2);  total (27/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.23e-01

TTTGCTGATGCATTCAGTAAAGATCCGGACTTCTACGCATTCATCCGTAGCCTGCGTGCTTATGAGAACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTGGTACTG  >  NC_000913/4404146‑4404430
                                                                                                                                          |                                                                                                                                                  
tttGCTGATGCATTCAGTAAAGATCCGGACTTCTACGCATTCATCCGTAGCCTGCGTGCTTATGAGAACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTcc                                                                                                                                            <  5:215268/147‑1 (MQ=255)
tttGCTGATGCATTCAGTAAAGATCCGGACTTCTACGCATTCATCCGTAGCCTGCGTGCTTATGAGAACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCg                                                                                                                                           >  6:79525/1‑148 (MQ=255)
     tGATGCATTCAGTAAAGATCCGGACTTCTACGCATTCATCCGTAGCCTGCGTGCTTATGAGAACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAAcgc                                                                                                                                     >  6:47196/1‑149 (MQ=255)
        tGCATTCAGTAAAGATCCGGACTTCTACGCATTCATCCGTAGCCTGCGTGCTTATGAGAACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGtt                                                                                                                                  <  7:36613/149‑1 (MQ=255)
            ttCAGTAAAGATCCGGACTTCTACGCATTCATCCGTAGCCTGCGTGCTTATGAGAACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAata                                                                                                                              >  5:76911/1‑149 (MQ=255)
            ttCAGTAAAGATCCGGACTTCTACGCATTCATCCGTAGCCTGCGTGCTTATGAGAACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAata                                                                                                                              >  1:50578/1‑149 (MQ=255)
                    agctCCGGACTTCTACGCATTCCTCCGTAGCCTGCGTGCTTATGAGAACAGCTTCGCTGGCACTCAGGCCGTGATGGGCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGAc                                                                                                                       <  6:35949/145‑1 (MQ=255)
                        ccGGACTTCTACGCATTCATCCGTAGCCTGCGTGCTTATGAGAACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTc                                                                                                                                             >  5:183484/1‑122 (MQ=255)
                        ccGGACTTCTACGCATTCATCCGTAGCCTGCGTGCTTATGAGAACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTc                                                                                                                                             <  6:183484/122‑1 (MQ=255)
                            aCTTCTACGCATTCATCCGTAGCCTGCGTGCTTATGAGAACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACa                                                                                                              <  2:224298/149‑1 (MQ=255)
                             cTTCTACGCATTCATCCGTAGCCTGCGTGCTTATGAGAACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAg                                                                                                             <  3:48275/149‑1 (MQ=255)
                                cTACGCATTCATCCGTAGCCTGCGTGCTTATGAGAACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGACGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGt                                                                                                           >  4:228918/1‑148 (MQ=255)
                                cTACGCATTCATCCGTAGCCTGCGTGCTTATGAGAACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGACGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGt                                                                                                           >  4:225469/1‑148 (MQ=255)
                                       ttCATCCGTAGCCTGCGTGCTTATGAGAACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAg                                                                                                   >  1:45494/1‑149 (MQ=255)
                                                 gCCTGCGTGCTTATGAGAACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATc                                                                                         >  1:193904/1‑149 (MQ=255)
                                                  ccTGCGTGCTTATGAGAACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATcc                                                                                        >  1:94443/1‑149 (MQ=255)
                                                      cGTGCTTATGAGAACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCa                                                                                     <  5:79525/148‑1 (MQ=255)
                                                            tATGAGAACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGt                                                                              >  1:780/1‑149 (MQ=255)
                                                                agaACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCggt                                                                          <  7:221915/149‑1 (MQ=255)
                                                                 gaACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCggtg                                                                         >  4:52763/1‑149 (MQ=255)
                                                                    cAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGtt                                                                      >  1:122454/1‑149 (MQ=255)
                                                                         tctcTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATaaa                                                                                                    >  7:26578/1‑114 (MQ=255)
                                                                         tctcTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATaaa                                                                                                    <  8:26578/114‑1 (MQ=255)
                                                                             tGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGttttctttttc                                                             >  3:227367/1‑149 (MQ=255)
                                                                             tGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGttttcttt                                                                <  2:50578/146‑1 (MQ=255)
                                                                               gCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGttttcttttt                                                              <  7:47054/146‑1 (MQ=255)
                                                                               gCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGACGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGttttcttttt                                                              <  3:228918/146‑1 (MQ=255)
                                                                               gCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGACGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGttttcttttt                                                              <  3:225469/146‑1 (MQ=255)
                                                                               gCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTGCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGGTAATCTAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGttttcttttt                                                              <  1:125691/146‑1 (MQ=255)
                                                                                cAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCta                                                           >  6:112548/1‑148 (MQ=255)
                                                                                 aaTCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGCTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGATTTCTTTTTCtat                                                          >  5:113367/1‑148 (MQ=255)
                                                                                  aTCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATaa                                                        <  5:224092/149‑1 (MQ=255)
                                                                                         cGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATaat                                                 >  7:72604/1‑149 (MQ=255)
                                                                                           tGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATaatga                                               <  5:47196/149‑1 (MQ=255)
                                                                                            gATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATaatgaa                                              >  1:175369/1‑149 (MQ=255)
                                                                                                         ccGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAAt                                 >  8:80529/1‑149 (MQ=255)
                                                                                                         ccGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAAt                                 <  5:112548/149‑1 (MQ=255)
                                                                                                         ccGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAAt                                 >  1:154539/1‑149 (MQ=255)
                                                                                                                gCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATctggctg                          <  2:780/149‑1 (MQ=255)
                                                                                                                 cGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGa                                     >  4:38483/1‑137 (MQ=255)
                                                                                                                 cGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGa                                     <  3:38483/137‑1 (MQ=255)
                                                                                                                 cGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATctggctgg                         >  2:20962/1‑149 (MQ=255)
                                                                                                                   aTTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCggt                                                                          <  6:156374/98‑1 (MQ=255)
                                                                                                                   aTTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCggt                                                                          >  5:156374/1‑98 (MQ=255)
                                                                                                                   aTTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGgcg                       <  8:72604/149‑1 (MQ=255)
                                                                                                                    tttCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGgcgc                      <  2:122454/149‑1 (MQ=255)
                                                                                                                         tcCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGcc                 >  5:16080/1‑149 (MQ=255)
                                                                                                                            gCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTg              <  4:227367/149‑1 (MQ=255)
                                                                                                                                    aaGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTgg       <  2:193904/148‑1 (MQ=255)
                                                                                                                                      gACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTGGTAc    >  7:241342/1‑149 (MQ=255)
                                                                                                                                       aCGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTGGTACt   <  6:76911/149‑1 (MQ=255)
                                                                                                                                        cGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTGGTACTg  >  3:161981/1‑149 (MQ=255)
                                                                                                                                          |                                                                                                                                                  
TTTGCTGATGCATTCAGTAAAGATCCGGACTTCTACGCATTCATCCGTAGCCTGCGTGCTTATGAGAACAGCTTCTCTGGCAATCAGGACGTGATGGTCATGAGCCCGGATAGCGATTTCTTCCGCTACATGAAGACGCCGACTTCCGCAACGCGTTAATATAACGACTGCGGTACAGGTCAATAAAGCCACCGCATCCTCAGGGATGTCGGTGGTTTTCTTTTTCTATAAGGATAATGAATGAATTCGACAATCTGGCTGGCGCTTGCCCTGGTTTTGGTACTG  >  NC_000913/4404146‑4404430

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