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% 124.6 / NA 29intergenic (‑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 . (14/15);  total (14/15)

TAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGC  >  NC_000913/2765703‑2765958
                                                                                                                          |                                                                                                                                     
taatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCt                                                                                                                      <  2:224826/139‑1 (MQ=255)
             tCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                                                                         <  1:106593/139‑1 (MQ=255)
                aGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCaa                                                                                                      >  2:231027/1‑139 (MQ=255)
                     ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATacac                                                                                                 >  1:678084/1‑139 (MQ=255)
                           ggTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgg                                                                                           >  2:6013/1‑139 (MQ=255)
                               agaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                                                                       <  1:6013/139‑1 (MQ=255)
                                 aTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTCTTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta                                                                                     <  2:595054/139‑1 (MQ=255)
                                               aCCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATATACTTCAg                                                                                            >  2:400584/1‑118 (MQ=255)
                                               aCCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                                                                            <  1:400584/118‑1 (MQ=255)
                                                           gACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                                                                            <  1:64527/106‑1 (MQ=255)
                                                           gACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                                                                            >  2:64527/1‑106 (MQ=255)
                                                               tcatcCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCa                                                                                                                             <  1:524050/69‑1 (MQ=255)
                                                               tcatcCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCa                                                                                                                             >  2:524050/1‑69 (MQ=255)
                                                                                                    aTCCTGCTCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTcc                                                                     >  2:550270/1‑88 (MQ=255)
                                                                                                    aTCCTGCCCTATAGTTGGAGTGCCCCTATTCGGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTcc                                                                     <  1:550270/88‑1 (MQ=255)
                                                                                                    aTCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATg                                                              >  1:570866/1‑95 (MQ=255)
                                                                                                    aTCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATg                                                              <  2:570866/95‑1 (MQ=255)
                                                                                                     tCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAgg                                             >  1:790368/1‑111 (MQ=255)
                                                                                                     tCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAgg                                             <  2:790368/111‑1 (MQ=255)
                                                                                                               taGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                                                                    <  1:361704/46‑1 (MQ=255)
                                                                                                               taGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                                                                    >  2:361704/1‑46 (MQ=255)
                                                                                                               taGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                                          <  2:538456/56‑1 (MQ=255)
                                                                                                               taGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                                          >  1:538456/1‑56 (MQ=255)
                                                                                                               taGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATg                                                      >  2:418241/1‑92 (MQ=255)
                                                                                                               taGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATg                                                      <  1:418241/92‑1 (MQ=255)
                                                                                                               taGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTc       <  2:790490/139‑1 (MQ=255)
                                                                                                                    ggAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                                                                       <  2:280222/54‑1 (MQ=255)
                                                                                                                    ggAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                                                                       >  1:280222/1‑54 (MQ=255)
                                                                                                                    ggAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGc  >  2:27353/1‑139 (MQ=255)
                                                                                                                          |                                                                                                                                     
TAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGC  >  NC_000913/2765703‑2765958

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: