Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,765,829 (C)5→4 intergenic (‑53/‑89) yfjM ← / → rnlA CP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, RNase LS

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,765,8250C.100.0% 206.3 / NA 49intergenic (‑49/‑93)yfjM/rnlACP4‑57 prophage, uncharacterized protein/CP4‑57 prophage, RNase LS
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (26/23);  total (26/23)

TGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACT  >  NC_000913/2765689‑2765942
                                                                                                                                         |                                                                                                                     
tGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGccc                                                                                                                     >  1:232677/1‑139 (MQ=255)
            aataatTACTAAGTATCGAGTTGCTGGTAGGTGAGATCAGGAGGAGGTTACCGAAACGTTGGACATCATCCCCCTATACGTGGACTGCCCCCCCGGTAAACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTgg                                                                                                        >  2:136493/1‑139 (MQ=255)
            aataatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTgg                                                                                                        <  2:232677/139‑1 (MQ=255)
            aataatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTgg                                                                                                        >  1:421010/1‑139 (MQ=255)
                                 tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATac                                                                                   <  2:421010/139‑1 (MQ=255)
                                   ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATacac                                                                                 >  1:82137/1‑139 (MQ=255)
                                          gTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                          <  2:421495/139‑1 (MQ=255)
                                           tGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCaa                                                                                      >  1:662502/1‑126 (MQ=255)
                                           tGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCaa                                                                                      <  2:662502/126‑1 (MQ=255)
                                                    gAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTCTTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTattatt                                                                >  2:817118/1‑139 (MQ=255)
                                                    gAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTattatt                                                                >  1:107260/1‑139 (MQ=255)
                                                       gagGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTa                                                             <  1:817118/139‑1 (MQ=255)
                                                         ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTaga                                                           >  2:1085621/1‑139 (MQ=255)
                                                         ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTaga                                                           >  2:68156/1‑139 (MQ=255)
                                                          gTTACCGAAACGATGGACATCATCCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCa                                                                                                             >  1:1002450/1‑89 (MQ=255)
                                                          gTTACCGAAACGATGGACATCATCCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCa                                                                                                             <  2:1002450/89‑1 (MQ=255)
                                                            tACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta                                                                     >  1:569087/1‑126 (MQ=255)
                                                            tACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta                                                                     <  2:569087/126‑1 (MQ=255)
                                                               cGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTcc                                                     >  1:24608/1‑139 (MQ=255)
                                                                 aaaCGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTcccc                                                   <  2:386168/139‑1 (MQ=255)
                                                                   aCGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCat                                                 <  2:107260/139‑1 (MQ=255)
                                                                                       aTTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGc                                                                                                                       >  1:780501/1‑51 (MQ=255)
                                                                                       aTTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGc                                                                                                                       <  2:780501/51‑1 (MQ=255)
                                                                                           tgtgGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGaaa                          <  2:24608/139‑1 (MQ=255)
                                                                                           tgtgGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGaaa                          <  1:1085621/139‑1 (MQ=255)
                                                                                                  tGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATaca                                                                                  <  2:120711/76‑1 (MQ=255)
                                                                                                  tGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATaca                                                                                  >  1:120711/1‑76 (MQ=255)
                                                                                                  tGCCACCCCGGTAGACGATGCTGCCCTATAGTTGGAGTGCCCCTATTCAGATGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAAt                   >  2:843476/1‑139 (MQ=255)
                                                                                                   gCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                                                    >  2:476548/1‑73 (MQ=255)
                                                                                                   gCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                                                    >  2:293893/1‑73 (MQ=255)
                                                                                                   gCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                                                    <  1:293893/73‑1 (MQ=255)
                                                                                                   gCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                                                    <  1:476548/73‑1 (MQ=255)
                                                                                                         ccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGt            >  1:656708/1‑139 (MQ=255)
                                                                                                            gTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                                                            >  1:710801/1‑72 (MQ=255)
                                                                                                            gTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                                                            <  2:710801/72‑1 (MQ=255)
                                                                                                                   aTCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTGCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACt  <  2:656708/139‑1 (MQ=255)
                                                                                                                   aTCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACt  >  1:572512/1‑139 (MQ=255)
                                                                                                                   aTCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACt  >  1:1082068/1‑139 (MQ=255)
                                                                                                                       tGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                                                       <  2:542680/66‑1 (MQ=255)
                                                                                                                       tGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                                                       >  1:542680/1‑66 (MQ=255)
                                                                                                                       tGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                          >  1:501931/1‑63 (MQ=255)
                                                                                                                       tGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                          <  2:501931/63‑1 (MQ=255)
                                                                                                                       tGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATc                                 >  2:956501/1‑104 (MQ=255)
                                                                                                                       tGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATc                                 <  1:956501/104‑1 (MQ=255)
                                                                                                                              taGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                                                       <  2:740005/59‑1 (MQ=255)
                                                                                                                              taGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                                                       >  1:864662/1‑59 (MQ=255)
                                                                                                                              taGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                                                       <  2:864662/59‑1 (MQ=255)
                                                                                                                              taGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                                                       >  1:740005/1‑59 (MQ=255)
                                                                                                                                   ggAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                                                    <  1:14957/41‑1 (MQ=255)
                                                                                                                                   ggAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                                                    >  2:14957/1‑41 (MQ=255)
                                                                                                                                   ggAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCa                                                  >  2:651691/1‑75 (MQ=255)
                                                                                                                                   ggAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCa                                                  <  1:651691/75‑1 (MQ=255)
                                                                                                                                         |                                                                                                                     
TGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACT  >  NC_000913/2765689‑2765942

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: