Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 281,050 (G)5→6 coding (934/1155 nt) yagA ← CP4‑6 prophage, putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913281,0451.G98.0% 186.1 / ‑3.6 51coding (939/1155 nt)yagACP4‑6 prophage, putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base G (23/27);  total (23/28)
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

CTGCCGTCTTCCTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑C‑GGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGT  >  NC_000913/280911‑281177
                                                                                                                                        |                                                                                                                                     
cTGCCGTCTTCCTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑C‑gggg                                                                                                                                   >  1:77718/1‑139 (MQ=255)
  gCCGTCTTCCTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑C‑gggggg                                                                                                                                 >  1:1413589/1‑138 (MQ=255)
          ccTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTgtgt                                                                                                                         <  1:859167/139‑1 (MQ=255)
             gCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGc                                                                                                                      <  1:335494/139‑1 (MQ=255)
               aTCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCg                                                                                                                    <  1:942767/139‑1 (MQ=255)
                            cGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCcgcg                                                                                                       <  1:550205/139‑1 (MQ=255)
                              aCCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCt                                                                                                     >  2:827574/1‑139 (MQ=255)
                              aCCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCt                                                                                                     >  1:210608/1‑139 (MQ=255)
                              aCCCGTTCCCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCt                                                                                                     >  1:980901/1‑139 (MQ=255)
                                             ttcGCCTTTCCTGCGCTCACACTTACCCCTTTCACGCTCATCTTCCCGCTGATCTCATCTTTCTTAACCATCACCCCTTCACCGTATTCTG‑GGGGGCG‑TTGTGCTGCCCCTGTACTGCCGCGCTGACGGCTGATACcgc                                                                                      <  2:76279/137‑1 (MQ=255)
                                                gCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAg                                                                                   <  2:77718/139‑1 (MQ=255)
                                                gCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAg                                                                                   <  2:607428/139‑1 (MQ=255)
                                                  cTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGTTGATATCCACTTTCCTGACCATCACCCCTTCTTCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGcc                                                                                 <  2:980901/139‑1 (MQ=255)
                                                       cTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTa                                                                            <  2:696712/139‑1 (MQ=255)
                                                       cTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTa                                                                            <  1:867630/139‑1 (MQ=255)
                                                          cgcTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCg                                                                         >  2:964355/1‑139 (MQ=255)
                                                               aGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTgtgt                                                                                                                         <  1:141/86‑1 (MQ=255)
                                                               aGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTgtgt                                                                                                                         >  2:141/1‑86 (MQ=255)
                                                               aGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCata                                                                    >  1:1101545/1‑139 (MQ=255)
                                                                  cTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATcc                                                                 <  2:1101545/139‑1 (MQ=255)
                                                                      ccccTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAgcg                                                             <  1:1140371/139‑1 (MQ=255)
                                                                           ttCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCg                                                        >  1:778787/1‑139 (MQ=255)
                                                                           ttCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCg                                                        >  1:273792/1‑139 (MQ=255)
                                                                               cGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCg                                                    >  1:1199452/1‑139 (MQ=255)
                                                                                gCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCgg                                                   >  2:550441/1‑139 (MQ=255)
                                                                                 cTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCggg                                                  >  2:694401/1‑139 (MQ=255)
                                                                                    aGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGt                                               >  2:1341100/1‑139 (MQ=255)
                                                                                       tttCCGCTGATATCCACTTTCCTGACCGTCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa                                            <  2:1177171/139‑1 (MQ=255)
                                                                                       tttCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa                                            <  1:1228624/139‑1 (MQ=255)
                                                                                       tttCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa                                            <  2:1413589/139‑1 (MQ=255)
                                                                                       tttCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa                                            <  2:1199452/139‑1 (MQ=255)
                                                                                       tttCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa                                            <  2:1242807/139‑1 (MQ=255)
                                                                                         tCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGtt                                              <  1:783100/135‑1 (MQ=255)
                                                                                         tCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGtt                                              >  2:783100/1‑135 (MQ=255)
                                                                                                atatCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAg                                   >  2:1369069/1‑139 (MQ=255)
                                                                                                atatCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAg                                   >  1:745366/1‑139 (MQ=255)
                                                                                                 tatCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGa                                  <  2:778787/139‑1 (MQ=255)
                                                                                                           ccTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTg                                                      <  2:1380310/109‑1 (MQ=255)
                                                                                                           ccTGACCATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTg                                                      >  1:1380310/1‑109 (MQ=255)
                                                                                                                 cATCACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa                                            >  1:970217/1‑113 (MQ=255)
                                                                                                                 cATCACCCCTTCATCGTATCC‑CGGGGGGCG‑TTGTGTTACCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGATCa                                            <  2:970217/113‑1 (MQ=255)
                                                                                                                    cgccccTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGaa               >  2:740775/3‑139 (MQ=255)
                                                                                                                    cACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGaa               >  1:1491778/1‑139 (MQ=255)
                                                                                                                    cACCCCTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGaa               >  2:54894/1‑139 (MQ=255)
                                                                                                                        ccTTCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGgcg           <  1:827574/139‑1 (MQ=255)
                                                                                                                          ttCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGgcgcg         <  1:1341100/139‑1 (MQ=255)
                                                                                                                           tCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa                                            >  1:1061557/1‑103 (MQ=255)
                                                                                                                           tCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa                                            <  2:1061557/103‑1 (MQ=255)
                                                                                                                           tCATCGTATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGgcgcgc        <  1:54894/139‑1 (MQ=255)
                                                                                                                             ttctatttC‑C‑GGGGGGGTTTGTTTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGAGAGCCCGGTTCACCCATATCCAGCGCTTCGTGCGGGCGTTCAAGGTTATAGGCCGTCCGCCAGTTGTCGAAGGCGCGCTg      <  2:761869/133‑1 (MQ=255)
                                                                                                                                 tATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa                                            >  2:784930/1‑97 (MQ=255)
                                                                                                                                 tATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa                                            <  1:784930/97‑1 (MQ=255)
                                                                                                                                 tATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGt  <  1:1201459/139‑1 (MQ=255)
                                                                                                                                 tATTC‑CGGGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGt  >  2:1178332/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                     
CTGCCGTCTTCCTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTC‑C‑GGGGGCG‑TTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGT  >  NC_000913/280911‑281177

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