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% 230.7 / NA 53intergenic (‑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/27);  total (26/27)

AATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGC  >  NC_000913/2765701‑2765958
                                                                                                                            |                                                                                                                                     
aataatTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTgg                                                                                                                        <  1:709654/139‑1 (MQ=255)
           aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGc                                                                                                                                       <  2:765499/114‑1 (MQ=255)
           aGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGc                                                                                                                                       >  1:765499/1‑114 (MQ=255)
             tATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCaa                                                                                                           >  1:689276/1‑139 (MQ=255)
               tCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                                                                         >  1:604164/1‑139 (MQ=255)
                cGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATc                                                                                                        >  2:251756/1‑139 (MQ=255)
                       ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATacac                                                                                                 >  2:271083/1‑139 (MQ=255)
                            tGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAg                                                                                                              <  2:536514/121‑1 (MQ=255)
                            tGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAg                                                                                                              >  1:536514/1‑121 (MQ=255)
                              gTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                                          <  1:251756/139‑1 (MQ=255)
                              gTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                                          <  2:395441/139‑1 (MQ=255)
                              gTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                                          <  2:362513/139‑1 (MQ=255)
                              gTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                                          <  2:689276/139‑1 (MQ=255)
                                 agaTCAGGAAGAGGTTACCGAAACCATGGACATCATCCCCCATGTTGTGGAGTACCCTCACGGGAGAAGATCCTGCCCTATAGTTGGAGTG‑CTCCTATTCAGGTGGCTAAATTTAGCCAATCAATACACTTCAgggtgg                                                                                       <  2:909190/139‑1 (MQ=255)
                                      aGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCAAAATTTAGCAAATCAATACACTTCAgggggg                                                                                       >  1:395751/1‑134 (MQ=255)
                                      aGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCAAAATTTAGCAAATCAATACACTTCAgggggg                                                                                       <  2:395751/134‑1 (MQ=255)
                                       ggAAGAGGTTACCGAAACGATGGACATCATCCCCTATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACATCAGGGGGGtattat                                                                                 <  2:604164/139‑1 (MQ=255)
                                       ggAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                                                                         >  2:170241/1‑115 (MQ=255)
                                       ggAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                                                                         <  1:170241/115‑1 (MQ=255)
                                        gAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                                                                    <  2:848723/119‑1 (MQ=255)
                                        gAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                                                                    >  1:848723/1‑119 (MQ=255)
                                             ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTaga                                                                           <  1:186220/139‑1 (MQ=255)
                                             ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTaga                                                                           >  1:473504/1‑139 (MQ=255)
                                              gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                                                                       <  1:846629/126‑1 (MQ=255)
                                              gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                                                                       >  2:846629/1‑126 (MQ=255)
                                              gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                                          <  2:53806/123‑1 (MQ=255)
                                              gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                                          >  1:53806/1‑123 (MQ=255)
                                               ttACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGt                                                                         >  2:1405/1‑139 (MQ=255)
                                                tACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCACTATACTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCCAATCAATACACTTCAGGGGGGTATTATTGTAGAGtt                                                                        >  1:584664/1‑139 (MQ=255)
                                                 aCCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGttt                                                                       <  1:85726/139‑1 (MQ=255)
                                                  ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                                                                         >  1:577136/1‑104 (MQ=255)
                                                  ccGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                                                                         <  2:577136/104‑1 (MQ=255)
                                                      aaCGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCa                                                                  <  2:653890/139‑1 (MQ=255)
                                                        cGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCata                                                                <  1:1405/139‑1 (MQ=255)
                                                                     cccccTTTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggg                                                                                         <  2:376586/101‑1 (MQ=255)
                                                                     cccccATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGCTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggg                                                                                         >  1:376586/1‑101 (MQ=255)
                                                                       ccctgtttGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCTATATGTTTCTATGGGATc                                                 <  1:612573/134‑1 (MQ=255)
                                                                       cccATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATc                                                 <  1:271083/139‑1 (MQ=255)
                                                                                tgGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATg                                        >  2:427362/1‑139 (MQ=255)
                                                                                    cTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACac                                    <  2:584664/139‑2 (MQ=255)
                                                                                       cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGt                                                             >  2:674008/1‑111 (MQ=255)
                                                                                       cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGt                                                             <  1:674008/111‑1 (MQ=255)
                                                                                       cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCa                                 >  2:702201/1‑139 (MQ=255)
                                                                                       cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCa                                 >  1:882877/1‑139 (MQ=255)
                                                                                          ccccGGTAGACGATCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGa                              >  1:315327/1‑139 (MQ=255)
                                                                                                      aTCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTcc                                                                     >  2:694696/1‑88 (MQ=255)
                                                                                                      aTCCTGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTcc                                                                     <  1:694696/88‑1 (MQ=255)
                                                                                                          tGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                                                                            <  2:299313/61‑1 (MQ=255)
                                                                                                          tGCCCTATAGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                                                                            >  1:299313/1‑61 (MQ=255)
                                                                                                                 taGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCa                                                                                                            <  1:800555/38‑1 (MQ=255)
                                                                                                                 taGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCa                                                                                                            >  2:800555/1‑38 (MQ=255)
                                                                                                                 taGTTGGAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTc       >  1:26472/1‑139 (MQ=255)
                                                                                                                      ggAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta                                                                                     <  1:927062/56‑1 (MQ=255)
                                                                                                                      ggAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta                                                                                     >  2:927062/1‑56 (MQ=255)
                                                                                                                      ggAGTG‑CCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGc  >  1:770607/1‑139 (MQ=255)
                                                                                                                            |                                                                                                                                     
AATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTCAGGGC  >  NC_000913/2765701‑2765958

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