Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,765,794	(C)7→6	intergenic (‑18/‑124)	yfjM ← / → rnlA	CP4‑57 prophage; protein YfjM/CP4‑57 prophage; RNase LS, toxin of the RnlAB toxin‑antitoxin system

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,765,788	0	C	.	100.0%	210.0 / NA	50	intergenic (‑12/‑130)	yfjM/rnlA	CP4‑57 prophage; protein YfjM/CP4‑57 prophage; RNase LS, toxin of the RnlAB toxin‑antitoxin system
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base . (12/38);  total (12/38)

GGGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATC  >  NC_000913/2765649‑2765911                                                                                                                                            |                                                                                                                            gggCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGc                                                                                                                             <  1:552754/140‑1 (MQ=255)  ggCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGcc                                                                                                                            <  1:311182/140‑1 (MQ=255)   gCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccc                                                                                                                           <  1:794413/140‑1 (MQ=255)    ctctctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGcccc                                                                                                                          <  1:705076/140‑1 (MQ=255)     tctctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccc                                                                                                                         <  2:290221/140‑1 (MQ=255)     tctctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccc                                                                                                                         <  1:150171/140‑1 (MQ=255)     tctctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccc                                                                                                                         <  2:148008/140‑1 (MQ=255)       tctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccccg                                                                                                                       >  2:35168/1‑139 (MQ=255)       tctGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGccccccg                                                                                                                       >  1:557008/1‑139 (MQ=255)            tCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGa                                                                                                                 <  1:246643/140‑1 (MQ=255)            tCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGa                                                                                                                 <  1:220894/140‑1 (MQ=255)                 cTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATc                                                                                                            <  1:362015/140‑1 (MQ=255)                 cTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATc                                                                                                            <  2:557268/140‑1 (MQ=255)                 cTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATc                                                                                                            <  2:727639/140‑1 (MQ=255)                 cTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATc                                                                                                            <  1:170952/140‑1 (MQ=255)                   gTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCt                                                                                                          <  1:415540/140‑1 (MQ=255)                       cGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGccc                                                                                                      <  2:555242/140‑1 (MQ=255)                       cGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGccc                                                                                                      <  1:416665/140‑1 (MQ=255)                       cGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGccc                                                                                                      <  1:72217/140‑1 (MQ=255)                              gAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGcc                                                                                                       <  1:744892/132‑1 (MQ=255)                              gAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGcc                                                                                                       >  2:744892/1‑132 (MQ=255)                              gAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGtt                                                                                               <  2:241809/140‑1 (MQ=255)                               aaTTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTg                                                                                              <  1:82423/140‑1 (MQ=255)                                   ttttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGt                                                                                          <  1:58653/140‑1 (MQ=255)                                   ttttCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGt                                                                                          <  2:763640/140‑1 (MQ=255)                                      tCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcc                                                                                       <  1:549859/140‑1 (MQ=255)                                      tCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGcc                                                                                       <  2:336354/140‑1 (MQ=255)                                               ctgGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCa                                                                              <  1:625729/140‑1 (MQ=255)                                                 ggCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATATTTGGAGTGCCCCCTATTCAgg                                                                            <  2:588137/140‑1 (MQ=255)                                                              tAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAg                                                               <  1:202228/140‑1 (MQ=255)                                                                   aTCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAAt                                                          <  2:368507/140‑1 (MQ=255)                                                                    tCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATc                                                         >  1:381524/1‑140 (MQ=255)                                                                     cGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCa                                                        >  1:625313/1‑140 (MQ=255)                                                                        gTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                     <  2:159224/140‑1 (MQ=255)                                                                        gTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                     <  1:35168/140‑1 (MQ=255)                                                                          tgctgcTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATaca                                                   >  2:498216/1‑140 (MQ=255)                                                                           gctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATacac                                                  >  2:804554/1‑140 (MQ=255)                                                                            ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACt                                                 >  1:681283/1‑140 (MQ=255)                                                                            ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACt                                                 >  1:227375/1‑140 (MQ=255)                                                                            ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACt                                                 >  2:748406/1‑140 (MQ=255)                                                                            ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCAT‑CCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                     <  1:321990/135‑1 (MQ=255)                                                                            ctgctTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCAT‑CCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATa                                                     >  2:321990/1‑135 (MQ=255)                                                                                ttGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                             <  2:381524/140‑1 (MQ=255)                                                                                ttGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                             <  1:550561/140‑1 (MQ=255)                                                                                ttGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAg                                             <  2:227375/140‑1 (MQ=255)                                                                                  ggTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAggg                                           <  2:249828/140‑1 (MQ=255)                                                                                   gTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggg                                          <  2:625313/140‑1 (MQ=255)                                                                                     gagaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAgggggg                                        <  2:557008/140‑1 (MQ=255)                                                                                       gaTCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta                                      <  1:529639/140‑1 (MQ=255)                                                                                             gAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTg                                <  2:661071/140‑1 (MQ=255)                                                                                                  ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta                                      <  2:413160/129‑1 (MQ=255)                                                                                                  ggTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGta                                      >  1:413160/1‑129 (MQ=255)                                                                                                       ccGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTcc                      >  2:605383/1‑140 (MQ=255)                                                                                                       ccGAAACGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTcc                      >  2:651579/1‑140 (MQ=255)                                                                                                          aaaCGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCa                   <  2:95603/140‑1 (MQ=255)                                                                                                           aaCGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCat                  <  1:230934/140‑1 (MQ=255)                                                                                                             cGATGGACATCATCCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCatat                <  2:409600/140‑1 (MQ=255)                                                                                                                     atcatcCCCCATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTAt        <  1:498216/140‑1 (MQ=255)                                                                                                                           ccccATTTTGTGGACTG‑CCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATc  <  2:681283/140‑1 (MQ=255)                                                                                                                                            |                                                                                                                            GGGCTCTCTGCTCGTTCTGTTTCGACGTAGAATTTTTTCTTGTCTGCTGGCGAATAATTACTAAGTATCGAGTTGCTGCTTGGTGAGATCAGGAAGAGGTTACCGAAACGATGGACATCATCCCCCATTTTGTGGACTGCCCCCCCGGTAGACGATCCTGCCCTATAGTTGGAGTGCCCCCTATTCAGGTGGCTAAATTTAGCAAATCAATACACTTCAGGGGGGTATTATTGTAGAGTTTCCCCATATGTTTCTATGGGATC  >  NC_000913/2765649‑2765911

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑