Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	573,440	C→T	R24C (CGC→TGC)	rusA →	DLP12 prophage, endonuclease RUS

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	573,440	0	C	T	100.0%	229.4 / NA	66	R24C (CGC→TGC)	rusA	DLP12 prophage, endonuclease RUS
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (37/29);  total (37/29)

AGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCG  >  NC_000913/573305‑573567                                                                                                                                        |                                                                                                                                aGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGtgca                                                                                                                              >  2:38221/1‑139 (MQ=255) aGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGtgca                                                                                                                              >  1:852164/1‑139 (MQ=255)     tGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGtgcacgca                                                                                                                          >  2:879526/1‑139 (MQ=255)       cgcgcTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGcacgcacg                                                                                                                        >  1:1551225/1‑139 (MQ=255)           cTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAg                                                                                                                    <  1:38221/139‑1 (MQ=255)                 ggTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCaga                                                                                                              <  2:1482497/139‑1 (MQ=255)                     tagCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCATATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAgggg                                                                                                          <  1:1186476/137‑1 (MQ=255)                      ggCGAGAACACAGG‑TATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGggca                                                                                                        >  1:105847/1‑139 (MQ=255)                                 aGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACcgc                                                                                              >  1:1252351/1‑139 (MQ=255)                                             tGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGcc                                                                                  >  1:712783/1‑139 (MQ=255)                                              gAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGccc                                                                                 <  2:725615/139‑1 (MQ=255)                                                aaGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGa                                                                               <  1:879026/139‑1 (MQ=255)                                                 agagGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGaa                                                                              <  2:1551225/139‑1 (MQ=255)                                                     ggggTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCAt                                                                          >  2:1604821/1‑139 (MQ=255)                                                     ggggTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCAt                                                                          >  2:1540447/1‑139 (MQ=255)                                                             taAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGggca                                                                                                        >  2:247392/1‑101 (MQ=255)                                                             taAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGggca                                                                                                        <  1:247392/101‑1 (MQ=255)                                                             taAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACg                                                                  <  2:852164/139‑1 (MQ=255)                                                              aaGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACg                                                                  >  1:378798/1‑138 (MQ=255)                                                              aaGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACg                                                                  <  2:378798/138‑1 (MQ=255)                                                                 gCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGggcaggca                                                                                                    <  1:1495517/101‑1 (MQ=255)                                                                 gCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGggcaggca                                                                                                    >  2:1495517/1‑101 (MQ=255)                                                                 gCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATc                                                                            >  2:1139415/1‑125 (MQ=255)                                                                 gCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATc                                                                            <  1:1139415/125‑1 (MQ=255)                                                                  cGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATg                                                             <  1:577299/139‑1 (MQ=255)                                                                    tGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACg                                                                  >  1:1671295/1‑132 (MQ=255)                                                                    tGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACg                                                                  <  2:1671295/132‑1 (MQ=255)                                                                    tGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCa                                                                >  1:1713709/1‑134 (MQ=255)                                                                    tGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCa                                                                <  2:1713709/134‑1 (MQ=255)                                                                    tGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCt                                                           >  1:1258494/1‑139 (MQ=255)                                                                      aaTGCCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCag                                                                                                               <  1:1054659/85‑1 (MQ=255)                                                                      aaTACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAGTCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTgg                                                         >  2:1089117/1‑139 (MQ=255)                                                                      aaTACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCag                                                                                                               >  2:1054659/1‑85 (MQ=255)                                                                         aCCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGata                                                      <  2:105847/139‑1 (MQ=255)                                                                              cAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGc                                                 <  2:712783/139‑1 (MQ=255)                                                                                    cacaTTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCa                                                                <  1:1456139/118‑1 (MQ=255)                                                                                    cacaTTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCa                                                                >  2:1456139/1‑118 (MQ=255)                                                                                      caTTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTAt                                         >  2:1415570/1‑139 (MQ=255)                                                                                         tACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGcc                                      >  1:511332/1‑139 (MQ=255)                                                                                                 ccTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGggcagg                                                                                                      >  1:1657421/1‑67 (MQ=255)                                                                                                 ccTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGggcagg                                                                                                      <  2:1657421/67‑1 (MQ=255)                                                                                                 ccTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGc                                                 >  2:100388/1‑120 (MQ=255)                                                                                                 ccTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGc                                                 <  1:100388/120‑1 (MQ=255)                                                                                                     cGAGCAATAATCGCTATTACCGCCATAATCGCGGGtgcacgca                                                                                                                          <  2:115078/43‑1 (MQ=255)                                                                                                     cGAGCAATAATCGCTATTACCGCCATAATCGCGGGtgcacgca                                                                                                                          >  1:115078/1‑43 (MQ=255)                                                                                                     cGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCCGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCg                          >  2:380629/1‑139 (MQ=255)                                                                                                     cGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGAt                                                                                           <  1:563333/74‑1 (MQ=255)                                                                                                     cGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGAt                                                                                           >  2:563333/1‑74 (MQ=255)                                                                                                     cGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATaa                                                                                         >  2:1768035/1‑76 (MQ=255)                                                                                                     cGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATaa                                                                                         <  1:1768035/76‑1 (MQ=255)                                                                                                     cGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCg                          >  2:333404/1‑139 (MQ=255)                                                                                                     cGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCg                          >  2:512885/1‑139 (MQ=255)                                                                                                     cGAGCAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCACGATAACGTCGCCAGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCg                          >  2:950524/1‑139 (MQ=255)                                                                                                      gAGCAATAATCGCTATTACCGCCCTAATCGCGGGTGCACGCACGTCAGCGCAGCGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGc                         >  1:935243/1‑139 (MQ=255)                                                                                                         cAATAATCGCTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCAtt                      <  2:511332/139‑1 (MQ=255)                                                                                                            taatCGCTATTACCGCCATAATCGCGGGTGCACGCACGt                                                                                                                       <  1:583581/39‑1 (MQ=255)                                                                                                            taatCGCTATTACCGCCATAATCGCGGGTGCACGCACGt                                                                                                                       >  2:583581/1‑39 (MQ=255)                                                                                                                  cTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACg                                                                  >  2:270212/1‑86 (MQ=255)                                                                                                                  cTATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACg                                                                  <  1:270212/86‑1 (MQ=255)                                                                                                                   tATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCcaca            >  2:1520423/1‑139 (MQ=255)                                                                                                                   tATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCcaca            <  2:1258494/139‑1 (MQ=255)                                                                                                                   tATTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCcaca            <  1:380629/139‑1 (MQ=255)                                                                                                                    aTTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGc                                                 >  1:834893/1‑101 (MQ=255)                                                                                                                    aTTACCGCCATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGc                                                 <  2:834893/101‑1 (MQ=255)                                                                                                                           ccATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGTAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGaaaa                              <  2:1076058/113‑1 (MQ=255)                                                                                                                           ccATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGaaaa                              >  1:1076058/1‑113 (MQ=255)                                                                                                                           ccATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGAt    >  2:687986/1‑139 (MQ=255)                                                                                                                           ccATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGAt    >  2:1719068/1‑139 (MQ=255)                                                                                                                            cATAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGAtc   >  2:1052274/1‑139 (MQ=255)                                                                                                                             aTAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGAtcg  >  2:1264227/1‑139 (MQ=255)                                                                                                                             aTAATCGCGGGTGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGAtcg  >  1:1575621/1‑139 (MQ=255)                                                                                                                                        |                                                                                                                                AGAATGCGCGCTGGAAGGTATGGCGAGAACACAGGTTATCTGGCTGAAAGAGGGGGTTATTAAGGCGTGAATACCTACAGCATCACATTACCCTGGCCTCCGAGCAATAATCGCTATTACCGCCATAATCGCGGGCGCACGCACGTCAGCGCAGAGGGGCAGGCATACCGCGATAACGTCGCCCGAATCATTAAAAACGCAATGCTGGATATCGGCCTGGCTATGCCTGTGAAAATCCGCATTGAGTGCCACATGCCGGATCG  >  NC_000913/573305‑573567

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