Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 284,869 C→T P557S (CCG→TCG)  yagF → CP4‑6 prophage, dehydratase family protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913284,8690CT98.3% 187.2 / NA 60P557S (CCG→TCG) yagFCP4‑6 prophage, dehydratase family protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (30/29);  minor base A (1/0);  total (31/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 = 1.00e+00

ACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGCCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGACGTTAACTGGCAGCGTGAACTTCATCGGCACCGCGGACAACCCGCTGACGCCGGAAGAGGGCGCGCGCGA  >  NC_000913/284734‑284993
                                                                                                                                       |                                                                                                                            
aCCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGa                                                                                                                           >  2:355296/1‑139 (MQ=255)
     ccAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGc                                                                                                                      >  1:777879/1‑139 (MQ=255)
     ccAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGc                                                                                                                      >  1:1490410/1‑139 (MQ=255)
        gCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAg                                                                                                                   <  2:1558259/139‑1 (MQ=255)
          tCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCt                                                                                                                 >  1:1609929/1‑139 (MQ=255)
              cTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAg                                                                                                                   >  1:1783280/1‑133 (MQ=255)
              cTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAg                                                                                                                   <  2:1783280/133‑1 (MQ=255)
              cTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTgcgc                                                                                                             <  2:777879/139‑1 (MQ=255)
              cTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTgcgc                                                                                                             >  1:1647789/1‑139 (MQ=255)
                   cgcTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATaa                                                                                                        <  2:1609929/139‑1 (MQ=255)
                       aaaGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGAc                                                                                                    <  1:355296/139‑1 (MQ=255)
                            atatCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGAcatcat                                                                                               <  2:717851/139‑1 (MQ=255)
                                 tcGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCgaga                                                                                          <  2:1490410/139‑1 (MQ=255)
                                    tGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTgcg                                                                                                              >  1:1169312/1‑116 (MQ=255)
                                    tGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTgcg                                                                                                              <  2:1169312/116‑1 (MQ=255)
                                    tGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACg                                                                                                      >  1:140993/1‑124 (MQ=255)
                                    tGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACg                                                                                                      <  2:140993/124‑1 (MQ=255)
                                                 tgtCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTgcgc                                                                                                             >  1:1430453/1‑104 (MQ=255)
                                                 tgtCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTgcgc                                                                                                             <  2:1430453/104‑1 (MQ=255)
                                                 tgtCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCgaga                                                                                          <  2:241836/123‑1 (MQ=255)
                                                 tgtCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGACGATTGGCAAGCTGCGCGATAACGACATCATCgaga                                                                                          >  1:241836/1‑123 (MQ=255)
                                                   tCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGAt                                                                                                          <  2:652494/105‑1 (MQ=255)
                                                   tCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGAt                                                                                                          >  1:652494/1‑105 (MQ=255)
                                                   tCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGa                                                                        <  2:298617/139‑1 (MQ=255)
                                                      cTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGACGt                                                                     >  2:1707564/1‑139 (MQ=255)
                                                           caCCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGACGTTAACt                                                                >  1:129249/1‑139 (MQ=255)
                                                                   cgcgcTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGACGTTAACTGGCAGCGt                                                        >  2:1269971/1‑139 (MQ=255)
                                                                     cgcTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGACGTTAACTGGCAGCGTGa                                                      <  2:792526/139‑1 (MQ=255)
                                                                     cgcTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGACGTTAACTGGCAGCGTGa                                                      >  1:792526/1‑139 (MQ=255)
                                                                               gCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTgcg                                                                                                              >  1:1751937/1‑73 (MQ=255)
                                                                               gCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTgcg                                                                                                              <  2:1751937/73‑1 (MQ=255)
                                                                                  tgtCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCgaga                                                                                          >  2:935071/1‑90 (MQ=255)
                                                                                  tgtCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCgaga                                                                                          <  1:1307505/90‑1 (MQ=255)
                                                                                  tgtCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCgaga                                                                                          >  2:1307505/1‑90 (MQ=255)
                                                                                  tgtCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCgaga                                                                                          <  1:935071/90‑1 (MQ=255)
                                                                                  tgtCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGa                                                                        >  2:417854/1‑108 (MQ=255)
                                                                                  tgtCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGa                                                                        <  1:417854/108‑1 (MQ=255)
                                                                                  tgtCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGACGTTAACTGGCAGCGTGAACTTCATCGGCAc                                         >  1:571232/1‑139 (MQ=255)
                                                                                    tCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCgaga                                                                                          <  2:220177/88‑1 (MQ=255)
                                                                                    tCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCgaga                                                                                          >  1:220177/1‑88 (MQ=255)
                                                                                       aCGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGcgcg                                                                                                            <  1:754338/67‑1 (MQ=255)
                                                                                       aCGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGcgcg                                                                                                            >  2:754338/1‑67 (MQ=255)
                                                                                         gggCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGAcatcat                                                                                               <  1:914143/78‑1 (MQ=255)
                                                                                         gggCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGAcatcat                                                                                               >  2:914143/1‑78 (MQ=255)
                                                                                              ccTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGACGTTAACTGGCAGCGTGAACTTCATCGGCACCGCGGACAAccc                             >  1:876076/1‑139 (MQ=255)
                                                                                                 gCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGACGTTAACTGGCAGCGTGAACTTCATCGGCACCGCGGACAACCCGCt                          <  2:12792/139‑1 (MQ=255)
                                                                                                            gtgtCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGAcatcat                                                                                               <  1:733585/59‑1 (MQ=255)
                                                                                                            gtgtCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGAcatcat                                                                                               >  2:733585/1‑59 (MQ=255)
                                                                                                            gtgtCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGACGTTAACTGGCAGCGTGAACTTCATCGGCACCGCGGACAACCCGCTGACGCCGGAag               <  2:730226/139‑1 (MQ=255)
                                                                                                             tgtCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGACGTTAACTg                                                               >  2:1573157/1‑90 (MQ=255)
                                                                                                             tgtCGCCGGAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGACGTTAACTg                                                               <  1:1573157/90‑1 (MQ=255)
                                                                                                                     gAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCgaga                                                                                          >  2:1148714/1‑55 (MQ=255)
                                                                                                                     gAGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCgaga                                                                                          <  1:1148714/55‑1 (MQ=255)
                                                                                                                      aGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGcgcg                                                                                                            >  2:546540/1‑36 (MQ=255)
                                                                                                                      aGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGcgcg                                                                                                            <  1:546540/36‑1 (MQ=255)
                                                                                                                      aGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCgaga                                                                                          >  1:1616240/1‑54 (MQ=255)
                                                                                                                      aGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCgaga                                                                                          <  2:1616240/54‑1 (MQ=255)
                                                                                                                      aGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGACGTTAACTg                                                               <  2:597980/81‑1 (MQ=255)
                                                                                                                      aGGCGCTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGACGTTAACTg                                                               >  1:597980/1‑81 (MQ=255)
                                                                                                                         cgcTGGCGGGCGGGTCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGACGTTAACTGGCAGCGTGAACTTCATCGGCACCGCGGACAACCCGCTGACGCCGGAAGAGGGCGCGCGCGa  >  2:22156/1‑139 (MQ=255)
                                                                                                                                       |                                                                                                                            
ACCTACCAGCTCACCTCCGCGCTAAAGCATATCTCGTGGGGCAAGACGGTGTCGCTCATCACCGATGCGCGCTTCTCGGGCGTGTCGACGGGCGCCTGCTTCGGCCACGTGTCGCCGGAGGCGCTGGCGGGCGGGCCGATTGGCAAGCTGCGCGATAACGACATCATCGAGATTGCCGTGGATCGTCTGACGTTAACTGGCAGCGTGAACTTCATCGGCACCGCGGACAACCCGCTGACGCCGGAAGAGGGCGCGCGCGA  >  NC_000913/284734‑284993

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