Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,643,501	T→C	T268T (ACA→ACG)	ydfU ←	Qin prophage; protein YdfU

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,643,501	0	T	C	100.0%	191.0 / NA	53	T268T (ACA→ACG)	ydfU	Qin prophage; protein YdfU
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (26/27);  total (26/27)

ATGCTTTCTGCACAAAGGCAACACAAAGAGATCATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTG  >  NC_000913/1643366‑1643637                                                                                                                                        |                                                                                                                                         aTGCTTTCTGCACAAAGGCAACACAAAGAGATCATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACgcgtgt                                                                                                                                      >  1:655758/1‑140 (MQ=255) aTGCTTTCTGCACAAAGGCAACACAAAGAGATCATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACgcgtgt                                                                                                                                      >  1:629024/1‑140 (MQ=255)    cTTTCTGCACAAAGGCAACACAAAGAGATCATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTAtt                                                                                                                                   <  2:856780/140‑1 (MQ=255)        cTGCACAAAGGCAACACAAAGAGATCATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTtctc                                                                                                                               <  1:98112/140‑1 (MQ=255)              aaaGGCAACACAAAGAGATCATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAAc                                                                                                                         >  2:142069/1‑140 (MQ=255)              aaaGGCAACACAAAGAGATCATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAAc                                                                                                                         >  2:237830/1‑140 (MQ=255)                aGGCAACACAAAGAGATCATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAAccc                                                                                                                       >  1:36954/1‑140 (MQ=255)                aGGCAACACAAAGAGATCATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAAccc                                                                                                                       >  1:56341/1‑140 (MQ=255)                      cacaAAGAGATCATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAgcggc                                                                                                                 <  2:50418/140‑1 (MQ=255)                      cacaAAGAGATCATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAgcggc                                                                                                                 <  2:743191/140‑1 (MQ=255)                      cacaAAGAGATCATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAgcggc                                                                                                                 <  2:655758/140‑1 (MQ=255)                      cacaAAGAGATCATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAgcggc                                                                                                                 <  2:629024/140‑1 (MQ=255)                           agagATCATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAgcggc                                                                                                                 >  1:45604/1‑135 (MQ=255)                           agagATCATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAgcggc                                                                                                                 <  2:45604/135‑1 (MQ=255)                            gagaTCATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGtt                                                                                                           >  2:215393/1‑140 (MQ=255)                               aTCATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAg                                                                                                        >  1:330909/1‑140 (MQ=255)                                tCATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAgg                                                                                                       <  2:36954/140‑1 (MQ=255)                                 cATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGc                                                                                                      >  1:476203/1‑140 (MQ=255)                                 cATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGc                                                                                                      >  2:71589/1‑140 (MQ=255)                                 cATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGc                                                                                                      >  1:904127/1‑140 (MQ=255)                                   tgtgCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAg                                                                                                    >  2:452578/1‑140 (MQ=255)                                            gTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGaa                                                                                           <  1:452578/140‑1 (MQ=255)                                                ccATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGAt                                                                                       <  1:71589/140‑1 (MQ=255)                                                  aTTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTc                                                                                     >  1:330772/1‑140 (MQ=255)                                                  aTTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGACTc                                                                                     >  2:292840/1‑140 (MQ=255)                                                     ccGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCgg                                                                                  <  2:330909/140‑1 (MQ=255)                                                         ccTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGAc                                                                              <  2:476203/140‑1 (MQ=255)                                                         ccTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGAc                                                                              <  2:56341/140‑1 (MQ=255)                                                           tGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGActc                                                                            <  2:317705/140‑1 (MQ=255)                                                             cccGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGga                                                                                 >  1:697133/1‑133 (MQ=255)                                                             cccGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGga                                                                                 <  2:697133/133‑1 (MQ=255)                                                                 tGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATc                                                                      <  2:330772/140‑1 (MQ=255)                                                                      aaTCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACgg                                                                 <  1:142069/140‑1 (MQ=255)                                                                       aTCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGCGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACgg                                                                 >  1:779597/1‑139 (MQ=255)                                                                       aTCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGCGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACgg                                                                 <  2:779597/139‑1 (MQ=255)                                                                        tCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTg                                                               <  2:440674/140‑1 (MQ=255)                                                                        tCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTg                                                               >  1:660106/1‑140 (MQ=255)                                                                        tCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTg                                                               <  2:904127/140‑1 (MQ=255)                                                                            atgatgCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCATCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGc                                                           <  1:824618/140‑1 (MQ=255)                                                                                 gCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTacc                                                      >  1:1011924/1‑140 (MQ=255)                                                                                  cGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTacca                                                     >  2:360738/1‑140 (MQ=255)                                                                                  cGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTacca                                                     <  1:237830/140‑1 (MQ=255)                                                                                    gATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTaccacc                                                   >  2:681604/1‑140 (MQ=255)                                                                                                      cGCAGCAAGCACACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGAt                                 >  2:820238/1‑140 (MQ=255)                                                                                                             agcaCACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTttgtt                          <  2:423153/140‑1 (MQ=255)                                                                                                             agcaCACGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTCCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTttgtt                          >  2:940218/1‑140 (MQ=255)                                                                                                                   cGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCtggt                    <  2:660106/140‑1 (MQ=255)                                                                                                                   cGGCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGACTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCtggt                    <  1:292840/140‑1 (MQ=255)                                                                                                                    ggCTGTGTCTTAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCtggtg                   >  1:56083/1‑140 (MQ=255)                                                                                                                              tAACCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGc         <  2:56083/140‑1 (MQ=255)                                                                                                                                aCCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCaa       >  1:496147/1‑140 (MQ=255)                                                                                                                                aCCCAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCaa       >  1:300180/1‑140 (MQ=255)                                                                                                                                   cAACGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACgg    >  2:975491/1‑140 (MQ=255)                                                                                                                                     aCGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGTCTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTg  >  1:276546/1‑140 (MQ=255)                                                                                                                                     aCGCGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCtggt                    <  1:1003363/122‑1 (MQ=255)                                                                                                                                        |                                                                                                                                         ATGCTTTCTGCACAAAGGCAACACAAAGAGATCATGTGCTTTTGTTCCCATTCCGCCCTGCCCGTGACCAATCAGATGATGCGGATCGTCGGCTGGCATACCGCAGCAAGCACACGGCTGTGTCTTAACCCAACGTGTGTATTTCTCCTTAACCCAGCGGCGACGTTTAGGCAGCTTCATGAAAGATTCCGGAGACTCTGGATCAACGGTGATGCTTACCACCGTCTTTTCCTGTGGTGATTTTTGTTGCTGGTGGGCGTAAGGCAACGGTG  >  NC_000913/1643366‑1643637

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