Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,770,404 G→A *234* (TAG→TAA)  yfjR → CP4‑57 prophage, putative DNA‑binding transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,770,4040GA100.0% 179.2 / NA 56*234* (TAG→TAA) yfjRCP4‑57 prophage, putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (29/27);  total (29/27)

GTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCGGAGTGTCCGACAGCC  >  NC_000913/2770295‑2770517
                                                                                                             |                                                                                                                 
gTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACAt                                                                                      >  2:765704/1‑139 (MQ=255)
     tgtCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATggagg                                                                                 >  1:1481309/1‑139 (MQ=255)
          ccTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAg                                                                                                              <  1:74453/105‑1 (MQ=255)
          ccTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAg                                                                                                              >  2:74453/1‑105 (MQ=255)
            tGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCg                                                                                                           >  1:811600/1‑106 (MQ=255)
            tGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCg                                                                                                           <  2:811600/106‑1 (MQ=255)
            tGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAg                                                                          >  2:655430/1‑139 (MQ=255)
             ggTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATggag                                                                                  >  1:1001764/1‑130 (MQ=255)
             ggTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATggag                                                                                  <  2:1001764/130‑1 (MQ=255)
              gTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTg                                                                        <  2:1481309/139‑1 (MQ=255)
                    gACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCgt                                                                                                          <  2:1109620/99‑1 (MQ=255)
                    gACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCgt                                                                                                          >  1:1109620/1‑99 (MQ=255)
                        cGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGc                                                                                           >  2:884192/1‑110 (MQ=255)
                        cGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGc                                                                                           <  1:884192/110‑1 (MQ=255)
                                 ttcaGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCg                                                     <  1:765704/139‑1 (MQ=255)
                                  tcaGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACATCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGc                                                    <  2:1083969/139‑1 (MQ=255)
                                  tcaGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGc                                                    <  2:1267554/139‑1 (MQ=255)
                                     ggCCTTATCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTGTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCgg                                                                                            <  1:954058/96‑1 (MQ=255)
                                     ggCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCgg                                                                                            >  2:954058/1‑96 (MQ=255)
                                           aCCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCgg                                                                                            <  2:147261/90‑1 (MQ=255)
                                           aCCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCgg                                                                                            >  1:147261/1‑90 (MQ=255)
                                                       cTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCgg                                                                                            >  1:1104991/1‑78 (MQ=255)
                                                       cTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCgg                                                                                            <  2:1104991/78‑1 (MQ=255)
                                                       cTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGc                                                                       >  1:1191216/1‑99 (MQ=255)
                                                       cTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGc                                                                       <  2:1191216/99‑1 (MQ=255)
                                                           aTCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAg                                                                                                              >  1:769993/1‑56 (MQ=255)
                                                           aTCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAg                                                                                                              <  2:769993/56‑1 (MQ=255)
                                                           aTCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCg                                                                                                           >  2:933134/1‑59 (MQ=255)
                                                           aTCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCg                                                                                                           <  1:933134/59‑1 (MQ=255)
                                                                    tCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCgg                                                                                            <  2:849040/65‑1 (MQ=255)
                                                                    tCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCgg                                                                                            >  1:849040/1‑65 (MQ=255)
                                                                      aCTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCa                                                                                          >  1:1217369/1‑65 (MQ=255)
                                                                      aCTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCa                                                                                          <  2:1217369/65‑1 (MQ=255)
                                                                      aCTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCgg                >  1:1190332/1‑139 (MQ=255)
                                                                      aCTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCgg                >  1:662733/1‑139 (MQ=255)
                                                                      aCTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAA‑‑‑‑GAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCa                                                  >  1:153676/1‑101 (MQ=255)
                                                                      aCTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAA‑‑‑‑GAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCa                                                  <  2:153676/101‑1 (MQ=255)
                                                                           ttACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCGGAgtgt           >  2:1326065/1‑139 (MQ=255)
                                                                            tACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATggag                                                                                  >  2:1499273/1‑67 (MQ=255)
                                                                            tACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATggag                                                                                  <  1:1499273/67‑1 (MQ=255)
                                                                             aCGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCGGAGTGTcc         <  1:655430/139‑1 (MQ=255)
                                                                                    ttCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTa                                                                                                     <  2:843883/40‑1 (MQ=255)
                                                                                    ttCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTa                                                                                                     >  1:843883/1‑40 (MQ=255)
                                                                                    ttCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACa                                                                                       >  2:1231651/1‑54 (MQ=255)
                                                                                    ttCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACa                                                                                       <  1:1231651/54‑1 (MQ=255)
                                                                                    ttCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGcc                                                                             >  2:811215/1‑64 (MQ=255)
                                                                                    ttCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGcc                                                                             <  1:811215/64‑1 (MQ=255)
                                                                                    ttCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCa                                                  >  1:821396/1‑91 (MQ=255)
                                                                                    ttCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCa                                                  >  1:1160418/1‑91 (MQ=255)
                                                                                    ttCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCa                                                  <  2:821396/91‑1 (MQ=255)
                                                                                    ttCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCa                                                  <  2:1160418/91‑1 (MQ=255)
                                                                                    ttCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCGGAGTGTCCGACAGcc  >  2:36900/1‑139 (MQ=255)
                                                                                             aCAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGc                                                                       >  2:348089/1‑61 (MQ=255)
                                                                                             aCAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGc                                                                       <  1:348089/61‑1 (MQ=255)
                                                                                                 cccACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGcc                                                                             <  2:1124423/51‑1 (MQ=255)
                                                                                                 cccACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGcc                                                                             >  1:1124423/1‑51 (MQ=255)
                                                                                                             |                                                                                                                 
GTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCGGAGTGTCCGACAGCC  >  NC_000913/2770295‑2770517

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