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% 176.9 / NA 42intergenic (‑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 . (21/21);  total (21/21)

ATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACT  >  NC_000913/2765705‑2765942
                                                                                                                        |                                                                                                                     
atTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGCCATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTaa                                                                                                    >  1:872726/1‑139 (MQ=255)
      aaGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                                              <  2:872726/139‑1 (MQ=255)
           tCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                                                         >  1:504356/1‑139 (MQ=255)
                   ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGc                                                                                                                       <  2:766138/102‑1 (MQ=255)
                   ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGc                                                                                                                       >  1:766138/1‑102 (MQ=255)
                                tCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                                              >  1:119677/1‑113 (MQ=255)
                                tCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                                              <  2:119677/113‑1 (MQ=255)
                                  aGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCa                                                                                            <  1:367081/113‑1 (MQ=255)
                                  aGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCa                                                                                            >  2:367081/1‑113 (MQ=255)
                                    gAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTattatt                                                                >  1:916077/1‑139 (MQ=255)
                                              ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                          <  2:722495/119‑1 (MQ=255)
                                              ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                          >  1:722495/1‑119 (MQ=255)
                                                   aCGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCat                                                 <  2:916077/139‑1 (MQ=255)
                                                              cgtcCCCCACTTTGTGGACTGCCCCCCCGTAAGACGCTCTTGCCCTACAGTCGGAGTGCCCCTATTGGTGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATg                                      <  1:597273/137‑1 (MQ=255)
                                                                 cccccATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGt                                                         <  2:585128/117‑1 (MQ=255)
                                                                 cccccATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCat                                                 >  2:618845/1‑125 (MQ=255)
                                                                 cccccATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCat                                                 <  1:618845/125‑1 (MQ=255)
                                                                 ccaccATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGt                                                         >  1:585128/4‑117 (MQ=255)
                                                                   cccATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGt                                             <  2:141290/127‑1 (MQ=255)
                                                                   cccATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGt                                             >  1:141290/1‑127 (MQ=255)
                                                                                 tGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCa                                                                                                             <  1:468941/49‑1 (MQ=255)
                                                                                 tGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCa                                                                                                             >  2:468941/1‑49 (MQ=255)
                                                                                 tGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCa                                                                                            <  2:246588/66‑1 (MQ=255)
                                                                                 tGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCa                                                                                            >  1:246588/1‑66 (MQ=255)
                                                                                 tGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAAt                   >  2:907229/1‑139 (MQ=255)
                                                                                   cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                                              >  1:357749/1‑62 (MQ=255)
                                                                                   cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                                              <  2:357749/62‑1 (MQ=255)
                                                                                   cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta                                                                     <  1:483842/87‑1 (MQ=255)
                                                                                   cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta                                                                     >  2:483842/1‑87 (MQ=255)
                                                                                   cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTa                                                             <  1:33350/95‑1 (MQ=255)
                                                                                   cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTa                                                             >  2:33350/1‑95 (MQ=255)
                                                                                   cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTcc                                                     <  2:990646/103‑1 (MQ=255)
                                                                                   cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTcc                                                     >  1:990646/1‑103 (MQ=255)
                                                                                   cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGt                                             >  1:889941/1‑111 (MQ=255)
                                                                                   cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGt                                             <  2:889941/111‑1 (MQ=255)
                                                                                        ccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGt            <  2:504356/139‑1 (MQ=255)
                                                                                                  aTCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACt  <  1:907229/139‑1 (MQ=255)
                                                                                                      tGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACaa       <  2:402549/130‑1 (MQ=255)
                                                                                                      tGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACaa       >  1:402549/1‑130 (MQ=255)
                                                                                                        cccTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                                              <  2:483469/41‑1 (MQ=255)
                                                                                                        cccTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                                              >  1:483469/1‑41 (MQ=255)
                                                                                                             taGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACt  >  1:959788/1‑128 (MQ=255)
                                                                                                             taGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACt  <  2:959788/128‑1 (MQ=255)
                                                                                                                  ggAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTAt                                       <  2:703878/86‑1 (MQ=255)
                                                                                                                  ggAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTAt                                       >  1:703878/1‑86 (MQ=255)
                                                                                                                        |                                                                                                                     
ATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACT  >  NC_000913/2765705‑2765942

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