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% 242.4 / NA 59intergenic (‑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 . (27/32);  total (27/32)

CTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTC  >  NC_000913/2765687‑2765953
                                                                                                                                          |                                                                                                                                
cTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGc                                                                                                                                  <  2:452339/139‑1 (MQ=255)
                  atTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTaa                                                                                                               >  1:383040/1‑139 (MQ=255)
                   ttACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAg                                                                                                                       >  1:1511043/1‑130 (MQ=255)
                   ttACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAg                                                                                                                       <  2:1511043/130‑1 (MQ=255)
                        aaGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                                                         <  2:878493/139‑1 (MQ=255)
                        aaGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                                                         <  2:28950/139‑1 (MQ=255)
                          gTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCa                                                                                                       <  2:1237266/139‑1 (MQ=255)
                           tATCGAGTTGCTGCTTGGTGAGATCAGGAAAAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGc                                                                                                                                  >  1:1408413/1‑112 (MQ=255)
                           tATCGAGTTGCTGCTTGGTGAGATCAGGAAAAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGc                                                                                                                                  <  2:1408413/112‑1 (MQ=255)
                             tCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                                                         <  2:234324/134‑1 (MQ=255)
                             tCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                                                         >  1:234324/1‑134 (MQ=255)
                             tCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                                                                    <  1:612547/139‑1 (MQ=255)
                                aGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCaa                                                                                                 >  2:1207516/1‑139 (MQ=255)
                                  ttGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                                                               <  1:937927/139‑1 (MQ=255)
                                            gTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                                     <  1:1207516/139‑1 (MQ=255)
                                                  tCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGCCTGCCCCCCCGGTAGACGATCCTGCCCTCTAGTTGGAGTGCCCCTATTCAGGTGGCTAATTTTAGCa                                                                                                       >  1:170717/1‑115 (MQ=255)
                                                  tCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGc                                                                                                                                  <  1:122829/89‑1 (MQ=255)
                                                  tCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGc                                                                                                                                  >  2:122829/1‑89 (MQ=255)
                                                  tCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCa                                                                                                       <  2:170717/115‑1 (MQ=255)
                                                      gAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                                                                  <  1:723985/132‑1 (MQ=255)
                                                      gAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                                                                  >  2:723985/1‑132 (MQ=255)
                                                           ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                                                         >  2:1089132/1‑104 (MQ=255)
                                                           ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                                                         <  1:1089132/104‑1 (MQ=255)
                                                           ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTaga                                                                      <  2:736615/139‑1 (MQ=255)
                                                           ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGACTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTaga                                                                      >  2:392398/1‑139 (MQ=255)
                                                            gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                                                                  >  2:688958/1‑126 (MQ=255)
                                                            gTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                                                                  <  1:688958/126‑1 (MQ=255)
                                                                                catcCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta                                                                                <  1:832781/108‑1 (MQ=255)
                                                                                catcCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta                                                                                >  2:832781/1‑108 (MQ=255)
                                                                                catcCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATg                                                 <  1:717109/139‑1 (MQ=255)
                                                                                catcCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGACTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATg                                                 <  1:392398/139‑1 (MQ=255)
                                                                                 atcCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCaa                                                                                                 >  2:309887/1‑90 (MQ=255)
                                                                                 atcCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCaa                                                                                                 <  1:309887/90‑1 (MQ=255)
                                                                                        aTTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAg                                         >  1:451128/1‑139 (MQ=255)
                                                                                          tttGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGAt                                             <  1:220371/133‑1 (MQ=255)
                                                                                          tttGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGAt                                             >  2:220371/1‑133 (MQ=255)
                                                                                            tgtgGACTGCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAAt                                    >  2:331799/1‑139 (MQ=255)
                                                                                                 aCTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAgg                                        >  2:1207622/1‑131 (MQ=255)
                                                                                                 aCTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAgg                                        <  1:1207622/131‑1 (MQ=255)
                                                                                                   tGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                                     <  1:899433/84‑1 (MQ=255)
                                                                                                   tGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                                                                     >  2:899433/1‑84 (MQ=255)
                                                                                                   tGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATg                                   >  2:1266856/1‑134 (MQ=255)
                                                                                                   tGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATg                                   <  1:1266856/134‑1 (MQ=255)
                                                                                                     cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTGc                                                                                                                                  <  1:87928/38‑1 (MQ=255)
                                                                                                     cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTGc                                                                                                                                  >  2:87928/1‑38 (MQ=255)
                                                                                                     cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                                                                    <  2:1282463/67‑1 (MQ=255)
                                                                                                     cccccccGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                                                                    >  1:1282463/1‑67 (MQ=255)
                                                                                                              tAGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACa                   <  2:451128/139‑1 (MQ=255)
                                                                                                               aGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTa                                                                        <  2:621456/85‑1 (MQ=255)
                                                                                                               aGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTa                                                                        >  1:621456/1‑85 (MQ=255)
                                                                                                               aGACGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACaa                  <  1:331799/139‑1 (MQ=255)
                                                                                                                  cGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                                                         >  2:944022/1‑49 (MQ=255)
                                                                                                                  cGATCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                                                         <  1:944022/49‑1 (MQ=255)
                                                                                                                    aTCCTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACt             <  2:383040/139‑1 (MQ=255)
                                                                                                                       cTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                                                                    >  2:579323/1‑49 (MQ=255)
                                                                                                                       cTGCCCTATAGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                                                                    <  1:579323/49‑1 (MQ=255)
                                                                                                                               taGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                                                         >  2:1367776/1‑36 (MQ=255)
                                                                                                                               taGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAg                                                                                                         <  1:1367776/36‑1 (MQ=255)
                                                                                                                               taGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCaa                                                                                                 >  2:477119/1‑44 (MQ=255)
                                                                                                                               taGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCaa                                                                                                 <  1:477119/44‑1 (MQ=255)
                                                                                                                               taGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                                                               >  2:1169751/1‑46 (MQ=255)
                                                                                                                               taGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                                                               <  1:1169751/46‑1 (MQ=255)
                                                                                                                               taGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTc  >  2:208824/1‑139 (MQ=255)
                                                                                                                               taGTTGGAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTc  >  2:563685/1‑139 (MQ=255)
                                                                                                                                    ggAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggg                                                                                    >  2:1318888/1‑52 (MQ=255)
                                                                                                                                    ggAGTGCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggg                                                                                    <  1:1318888/52‑1 (MQ=255)
                                                                                                                                          |                                                                                                                                
CTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATCCAGGAAATGACAATCAGGAGTTACAAAAACTTAAATCTGGTC  >  NC_000913/2765687‑2765953

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