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% 142.7 / NA 45*234* (TAG→TAA) yfjRCP4‑57 prophage, putative DNA‑binding transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (23/22);  total (23/22)

GCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCGGAGTGTCCGACAGCC  >  NC_000913/2770269‑2770517
                                                                                                                                       |                                                                                                                 
gCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAacgc                                                                                                                <  2:526819/139‑1 (MQ=255)
                    aTGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCgg                                                                                            <  1:814439/139‑1 (MQ=255)
                       aCAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCaa                                                                                         <  2:1822394/139‑1 (MQ=255)
                               tgtCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATggagg                                                                                 >  1:643189/1‑139 (MQ=255)
                                    ccTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCg                                                                                                           >  1:1131587/1‑108 (MQ=255)
                                    ccTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCg                                                                                                           <  2:1131587/108‑1 (MQ=255)
                                            aaGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCg                                                                    >  1:1327641/1‑139 (MQ=255)
                                                      cTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCAt                                                          <  2:1327641/139‑1 (MQ=255)
                                                           ttcaGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCg                                                     <  2:643189/139‑1 (MQ=255)
                                                                     aCCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAg                                                                                                              <  1:426615/72‑1 (MQ=255)
                                                                     aCCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAg                                                                                                              >  2:426615/1‑72 (MQ=255)
                                                                     aCCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAAcc                                                                                               <  2:1529472/87‑1 (MQ=255)
                                                                     aCCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAAcc                                                                                               >  1:1529472/1‑87 (MQ=255)
                                                                          aTTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTg                                                              >  1:198109/1‑115 (MQ=255)
                                                                          aTTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTg                                                              <  2:198109/115‑1 (MQ=255)
                                                                          aTTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTAt                                                       <  2:918798/122‑1 (MQ=255)
                                                                          aTTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTAt                                                       >  1:918798/1‑122 (MQ=255)
                                                                           ttttCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCgg                                                                                            <  2:1547262/84‑1 (MQ=255)
                                                                           ttttCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCgg                                                                                            >  1:1547262/1‑84 (MQ=255)
                                                                               cGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGcc                                                                      <  1:1793694/102‑1 (MQ=255)
                                                                               cGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGcc                                                                      >  2:1793694/1‑102 (MQ=255)
                                                                                     aTCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGcc                                                   <  2:1659528/115‑1 (MQ=255)
                                                                                     aTCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGcc                                                   >  1:1659528/1‑115 (MQ=255)
                                                                                     aTCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTg                           >  2:1357430/1‑139 (MQ=255)
                                                                                      tCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGc                                                                       >  2:274327/1‑94 (MQ=255)
                                                                                      tCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGc                                                                       <  1:274327/94‑1 (MQ=255)
                                                                                       ccAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGaa                         <  1:29899/139‑1 (MQ=255)
                                                                                                aCTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGctcctcca                                                                <  2:399091/91‑8 (MQ=255)
                                                                                                aCTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGctcctcca                                                                >  1:399091/1‑84 (MQ=255)
                                                                                                aCTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCg                                                     >  1:1373183/1‑102 (MQ=255)
                                                                                                aCTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCg                                                     <  2:1373183/102‑1 (MQ=255)
                                                                                                aCTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGaa                         >  1:1580971/1‑130 (MQ=255)
                                                                                                aCTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGaa                         <  2:1580971/130‑1 (MQ=255)
                                                                                                  tGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAacg                                                                                                                 <  2:1105361/40‑1 (MQ=255)
                                                                                                   gCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACa                                                                                       <  2:304010/65‑1 (MQ=255)
                                                                                                   gCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACa                                                                                       >  1:304010/1‑65 (MQ=255)
                                                                                                    cTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCGGAgtg            >  2:2522/1‑139 (MQ=255)
                                                                                                     ttACGTTTGTTCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCGGAgtgt           >  1:1331045/1‑139 (MQ=255)
                                                                                                              ttCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATggag                                                                                  >  2:1560004/1‑59 (MQ=255)
                                                                                                              ttCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATggag                                                                                  <  1:1560004/59‑1 (MQ=255)
                                                                                                              ttCCGGCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCGGAGTGTCCGACAGcc  >  2:382014/1‑139 (MQ=255)
                                                                                                                   gCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGt                                                                         >  1:951993/1‑63 (MQ=255)
                                                                                                                   gCCGACAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGt                                                                         <  2:951993/63‑1 (MQ=255)
                                                                                                                        cAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCgg                                                                                            <  2:1611288/39‑1 (MQ=38)
                                                                                                                        cAGCCCACCGGAATAACGCAGGCGTTGTTATCAACCCgg                                                                                            >  1:1611288/1‑39 (MQ=255)
                                                                                                                                cGGAATAACGCAGGCGTTGTTATCAACCCGGCAACa                                                                                       <  1:1035972/36‑1 (MQ=255)
                                                                                                                                cGGAATAACGCAGGCGTTGTTATCAACCCGGCAACa                                                                                       >  2:1035972/1‑36 (MQ=37)
                                                                                                                                       |                                                                                                                 
GCAGGAGACTTTTCGCCGCAATGACAGTCTGTGTCGCCTGGTTGAAGACCCGGTCTTCATTCAGGCCTTACCCCATTTTCGCTTTATCCAGCATTCACTGCTTACGTTTGTTCCGGCCGACAGCCCACCGGAATAGCGCAGGCGTTGTTATCAACCCGGCAACATGGAGGAGCCCAGTGCCCGTTATTGCCATTATCGCCATTGTTATCATCGTCATCATTCTGAACAAAACCGGAGTGTCCGACAGCC  >  NC_000913/2770269‑2770517

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