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

CTGCCGTCTTCCTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGT  >  NC_000913/280911‑281177
                                                                                                                                        |                                                                                                                                     
cTGCCGTCTTCCTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCgggg                                                                                                                                   >  1:77718/1‑139 (MQ=255)
  gCCGTCTTCCTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCgggggg                                                                                                                                 >  1:1413589/1‑138 (MQ=255)
          ccTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTgtgt                                                                                                                         <  1:859167/139‑1 (MQ=255)
             gCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGc                                                                                                                      <  1:335494/139‑1 (MQ=255)
               aTCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCg                                                                                                                    <  1:942767/139‑1 (MQ=255)
                            cGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCcgcg                                                                                                       <  1:550205/139‑1 (MQ=255)
                              aCCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCt                                                                                                     >  2:827574/1‑139 (MQ=255)
                              aCCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCt                                                                                                     >  1:210608/1‑139 (MQ=255)
                              aCCCGTTCCCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCt                                                                                                     >  1:980901/1‑139 (MQ=255)
                                             ttcGCCTTTCCTGCGCTCACACTTACCCCTTTCACGCTCATCTTCCCGCTGATCTCATCTTTCTTAACCATCACCCCTTCACCGTATTCTGGGGGGCGTTGTGCTGCCCCTGTACTGCCGCGCTGACGGCTGATACcgc                                                                                      <  2:76279/137‑1 (MQ=255)
                                                gCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAg                                                                                   <  2:77718/139‑1 (MQ=255)
                                                gCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAg                                                                                   <  2:607428/139‑1 (MQ=255)
                                                  cTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGTTGATATCCACTTTCCTGACCATCACCCCTTCTTCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGcc                                                                                 <  2:980901/139‑1 (MQ=255)
                                                       cTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTa                                                                            <  2:696712/139‑1 (MQ=255)
                                                       cTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTa                                                                            <  1:867630/139‑1 (MQ=255)
                                                          cgcTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCg                                                                         >  2:964355/1‑139 (MQ=255)
                                                               aGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTgtgt                                                                                                                         <  1:141/86‑1 (MQ=255)
                                                               aGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTgtgt                                                                                                                         >  2:141/1‑86 (MQ=255)
                                                               aGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCata                                                                    >  1:1101545/1‑139 (MQ=255)
                                                                  cTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATcc                                                                 <  2:1101545/139‑1 (MQ=255)
                                                                      ccccTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAgcg                                                             <  1:1140371/139‑1 (MQ=255)
                                                                           ttCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCg                                                        >  1:778787/1‑139 (MQ=255)
                                                                           ttCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCg                                                        >  1:273792/1‑139 (MQ=255)
                                                                               cGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCg                                                    >  1:1199452/1‑139 (MQ=255)
                                                                                gCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCgg                                                   >  2:550441/1‑139 (MQ=255)
                                                                                 cTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCggg                                                  >  2:694401/1‑139 (MQ=255)
                                                                                    aGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGt                                               >  2:1341100/1‑139 (MQ=255)
                                                                                       tttCCGCTGATATCCACTTTCCTGACCGTCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa                                            <  2:1177171/139‑1 (MQ=255)
                                                                                       tttCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa                                            <  1:1228624/139‑1 (MQ=255)
                                                                                       tttCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa                                            <  2:1413589/139‑1 (MQ=255)
                                                                                       tttCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa                                            <  2:1199452/139‑1 (MQ=255)
                                                                                       tttCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa                                            <  2:1242807/139‑1 (MQ=255)
                                                                                         tCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGtt                                              <  1:783100/135‑1 (MQ=255)
                                                                                         tCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGtt                                              >  2:783100/1‑135 (MQ=255)
                                                                                                atatCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAg                                   >  2:1369069/1‑139 (MQ=255)
                                                                                                atatCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAg                                   >  1:745366/1‑139 (MQ=255)
                                                                                                 tatCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGa                                  <  2:778787/139‑1 (MQ=255)
                                                                                                           ccTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTg                                                      <  2:1380310/109‑1 (MQ=255)
                                                                                                           ccTGACCATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTg                                                      >  1:1380310/1‑109 (MQ=255)
                                                                                                                 cATCACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa                                            >  1:970217/1‑113 (MQ=255)
                                                                                                                 cATCACCCCTTCATCGTATCCCGGGGGGCGTTGTGTTACCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGATCa                                            <  2:970217/113‑1 (MQ=255)
                                                                                                                    cgccccTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGaa               >  2:740775/3‑139 (MQ=255)
                                                                                                                    cACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGaa               >  1:1491778/1‑139 (MQ=255)
                                                                                                                    cACCCCTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGaa               >  2:54894/1‑139 (MQ=255)
                                                                                                                        ccTTCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGgcg           <  1:827574/139‑1 (MQ=255)
                                                                                                                          ttCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGgcgcg         <  1:1341100/139‑1 (MQ=255)
                                                                                                                           tCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa                                            >  1:1061557/1‑103 (MQ=255)
                                                                                                                           tCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa                                            <  2:1061557/103‑1 (MQ=255)
                                                                                                                           tCATCGTATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGgcgcgc        <  1:54894/139‑1 (MQ=255)
                                                                                                                             ttctatttCCGGGGGGGTTTGTTTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGAGAGCCCGGTTCACCCATATCCAGCGCTTCGTGCGGGCGTTCAAGGTTATAGGCCGTCCGCCAGTTGTCGAAGGCGCGCTg      <  2:761869/133‑1 (MQ=255)
                                                                                                                                 tATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa                                            >  2:784930/1‑97 (MQ=255)
                                                                                                                                 tATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCa                                            <  1:784930/97‑1 (MQ=255)
                                                                                                                                 tATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGt  <  1:1201459/139‑1 (MQ=255)
                                                                                                                                 tATTCCGGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGt  >  2:1178332/1‑139 (MQ=255)
                                                                                                                                        |                                                                                                                                     
CTGCCGTCTTCCTGCATCTCCTTCAGCCCGACCCGTTCTCCCCTGAACGCCTTGCCTGCGCTCAGACTTACCCCTTTCACGCTCAGCTTTCCGCTGATATCCACTTTCCTGACCATCACCCCTTCATCGTATTCCGGGGGCGTTGTGTTGCCGCTGTACTGCCGCGCTGACGGCTGATACCGCGAGCCCGGTACCGCCATATCCAGCGCCTCGTGCGGGCGTTCAAGGTTATAGACCGTCCGCCAGTGGTCGAAGGCGCGCTGCAGT  >  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: