Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,769,486 G→A *274* (TAG→TAA)  yfjQ → CP4‑57 prophage, uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,769,4860GA100.0% 188.0 / NA 59*274* (TAG→TAA) yfjQCP4‑57 prophage, uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (31/28);  total (31/28)

AGGGCGGACTCAGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAGACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTG  >  NC_000913/2769353‑2769600
                                                                                                                                     |                                                                                                                  
aGGGCGGACTCAGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGtt                                                                                                               >  1:1181365/1‑139 (MQ=255)
      gACTCAGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGt                                                                                                         >  2:55105/1‑139 (MQ=255)
          cAGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGcc                                                                                                     <  2:934154/139‑1 (MQ=255)
          cAGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGcc                                                                                                     <  1:55105/139‑1 (MQ=255)
             tGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACg                                                                                                  <  2:451747/139‑1 (MQ=255)
                   tAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTCCTCCCGAAACTGCAGTAAACGTTTCAGGTGGCCACGTTTTTa                                                                                            >  1:1158998/1‑139 (MQ=255)
                     aTGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAat                                                                                          >  1:1039399/1‑139 (MQ=255)
                      tGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAata                                                                                         >  2:674204/1‑139 (MQ=255)
                      tGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAata                                                                                         <  2:1181365/139‑1 (MQ=255)
                      tGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAata                                                                                         <  2:1311156/139‑1 (MQ=255)
                         tAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCg                                                                                      <  1:559975/139‑1 (MQ=255)
                                ggACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCa                                                                                                            >  2:730122/1‑110 (MQ=255)
                                ggACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCa                                                                                                            <  1:730122/110‑1 (MQ=255)
                                 gACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATCTGAAGCTTAACCGTGCGCTGTGGGTGATGGCGGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGAcaccac                                                                              >  1:827255/1‑139 (MQ=255)
                                           tACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCa                                                                                                            >  1:211895/1‑99 (MQ=255)
                                           tACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCa                                                                                                            <  2:211895/99‑1 (MQ=255)
                                           tACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCAt                                                                    <  2:787837/139‑1 (MQ=255)
                                                gtgccgtgcGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCa                                                                                                    <  1:1339334/102‑1 (MQ=255)
                                                gtgccgtgcGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCa                                                                                                    >  2:1339334/1‑102 (MQ=255)
                                                gtgccgtgcGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTa                                                               >  2:1073154/1‑139 (MQ=255)
                                                  gccgtgcGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTAtg                                                             >  1:1381/1‑139 (MQ=255)
                                                     gtgcgtgGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGc                                                          <  2:1381/139‑1 (MQ=255)
                                                        cgtgGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                                   <  2:197843/95‑1 (MQ=255)
                                                        cgtgGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                                   >  1:197843/1‑95 (MQ=255)
                                                        cgtgGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACg                                                                                                  <  2:198156/96‑1 (MQ=255)
                                                        cgtgGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACg                                                                                                  >  1:198156/1‑96 (MQ=255)
                                                          tgGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCgtg                                                     <  1:1073154/139‑1 (MQ=255)
                                                           ggCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCgtgt                                                    >  2:1075867/1‑139 (MQ=255)
                                                                  aCGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGcc                                                                                                     >  1:349636/1‑83 (MQ=255)
                                                                  aCGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGcc                                                                                                     <  2:349636/83‑1 (MQ=255)
                                                                      ggATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATccc                                         <  2:1039399/139‑1 (MQ=255)
                                                                             aaGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                                   >  2:392460/1‑74 (MQ=255)
                                                                             aaGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                                   <  1:392460/74‑1 (MQ=255)
                                                                                  tAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                                   <  1:1043624/69‑1 (MQ=255)
                                                                                  tAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                                   >  2:1043624/1‑69 (MQ=255)
                                                                                  tAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCt                             >  2:765889/1‑139 (MQ=255)
                                                                                    aCCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCa                                                                                                    >  1:1297963/1‑66 (MQ=255)
                                                                                    aCCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCa                                                                                                    <  2:1297963/66‑1 (MQ=255)
                                                                                       gTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                                   >  2:861539/1‑64 (MQ=255)
                                                                                       gTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                                   <  1:861539/64‑1 (MQ=255)
                                                                                       gTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTaa                        >  2:408716/1‑139 (MQ=255)
                                                                                        tGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATccc                                         >  2:38807/1‑121 (MQ=255)
                                                                                        tGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATccc                                         <  1:38807/121‑1 (MQ=255)
                                                                                          cgcTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTa                                                                                            >  2:1282712/1‑68 (MQ=255)
                                                                                          cgcTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTa                                                                                            <  1:1282712/68‑1 (MQ=255)
                                                                                              gtgGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGt                 <  1:1075867/139‑1 (MQ=255)
                                                                                                       ggCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACg                                                                                                  >  1:728521/1‑49 (MQ=255)
                                                                                                       ggCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACg                                                                                                  <  2:728521/49‑1 (MQ=255)
                                                                                                       ggCAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTgg        >  1:959629/1‑139 (MQ=255)
                                                                                                         cAGAAACACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCt      >  1:554841/1‑139 (MQ=255)
                                                                                                            aaaCACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                                   <  1:1278769/43‑1 (MQ=255)
                                                                                                            aaaCACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCAc                                                                                                   >  2:1278769/1‑43 (MQ=255)
                                                                                                             aaCACTGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTg  >  1:287741/1‑139 (MQ=255)
                                                                                                                  tGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACgttgt                                                                                              <  1:500414/42‑1 (MQ=255)
                                                                                                                  tGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACgttgt                                                                                              >  2:500414/1‑42 (MQ=255)
                                                                                                                  tGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCt                                                                <  1:511222/72‑1 (MQ=255)
                                                                                                                  tGCTCACGCAACTGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCt                                                                >  2:511222/1‑72 (MQ=255)
                                                                                                                              tGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTtat                                 <  1:1156424/91‑1 (MQ=255)
                                                                                                                              tGCAGTAAACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTtat                                 >  2:1156424/1‑91 (MQ=255)
                                                                                                                                     |                                                                                                                  
AGGGCGGACTCAGTGGGCGTAATGCTAAAGGCGGACGAACTCATACTCGTGCCGTGCGTGGCATTGACGGGGATGTGAAGCTTAACCGTGCGCTGTGGGTGATGGCAGAAACACTGCTCACGCAACTGCAGTAGACGTTTCATGTTGCCACGTTGTTAATATCGGACACCACCTGTCCGCATCGCTATGTGCTCGTGTGCCTCAATCCCCCGGTTATAGCTTTTAACCCCCGTTACATCTGGCTTTTG  >  NC_000913/2769353‑2769600

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: ‑