Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,639,617	T→C	G147G (GGA→GGG)	rrrQ ←	Qin prophage; putative lysozyme

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,639,617	0	T	C	100.0%	201.6 / NA	56	G147G (GGA→GGG)	rrrQ	Qin prophage; putative lysozyme
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base C (29/27);  total (29/27)

ATTTCCGACGATAAATCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGTCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTTGTGGTTCG  >  NC_000913/1639471‑1639764                                                                                                                                                   |                                                                                                                                                    aTTTCCGACGATAAATCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGc                                                                                                                                                 >  2:49353/1‑151 (MQ=255)   ttCCGACGATAAATCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGcgcccgcc                                                                                                                                                >  2:780451/1‑150 (MQ=255)             aaaTCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATcc                                                                                                                                     >  1:76548/1‑151 (MQ=255)              aaTCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATcc                                                                                                                                     >  2:168761/1‑150 (MQ=255)                  cTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATccacc                                                                                                                                  >  1:707867/1‑149 (MQ=255)                         cacGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCg                                                                                                                         <  1:245090/151‑1 (MQ=255)                          acGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGc                                                                                                                        >  2:94606/1‑151 (MQ=255)                           cGAGAAATATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCt                                                                                                                       <  2:853493/151‑1 (MQ=255)                            gagaAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCtt                                                                                                                      >  1:49651/1‑151 (MQ=255)                               aaaaaTCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCg                                                                                                                    <  1:427416/150‑1 (MQ=255)                                aaaaTCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGc                                                                                                                   >  2:549860/1‑150 (MQ=255)                                      aCCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCAcc                                                                                                            <  1:726735/151‑1 (MQ=255)                                      aCCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCAcc                                                                                                            >  2:790598/1‑151 (MQ=255)                                        cATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCtt                                                                                                          <  2:674560/151‑1 (MQ=255)                                          tGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCttt                                                                                                         >  1:329559/1‑150 (MQ=255)                                                 cTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCAc                                                                                                 >  1:953145/1‑151 (MQ=255)                                                    aTTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCAcc                                                                                                >  1:157257/1‑149 (MQ=255)                                                    aTTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCAcc                                                                                                >  1:515567/1‑149 (MQ=255)                                                          tGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTc                                                                                         >  2:284182/1‑150 (MQ=255)                                                          tGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTAAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTc                                                                                         >  1:989048/1‑150 (MQ=255)                                                            ttCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAg                                                                                       <  2:376172/150‑1 (MQ=255)                                                                    cccAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTtat                                                                               >  1:372850/1‑150 (MQ=255)                                                                        gcagGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAAc                                                                          >  1:936249/1‑151 (MQ=255)                                                                         cagGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACg                                                                         <  1:270026/151‑1 (MQ=255)                                                                          agGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGt                                                                        >  2:464244/1‑151 (MQ=255)                                                                               aaTGCACTCACCTGGTCCCGACGAATAGCCTGAACGTAACAGTTATTTGAACGATTGCGGCAATCGCGCCCGCCATCCTTAACCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCgcc                                                                    <  2:276754/150‑3 (MQ=255)                                                                                         ccTGGTCACGACGAATACCCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACtt                                                         <  1:596775/151‑1 (MQ=255)                                                                                           tGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTAc                                                       >  1:270520/1‑151 (MQ=255)                                                                                            ggTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTa                                                        <  2:509808/149‑1 (MQ=255)                                                                                             gTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCg                                                     <  2:522167/151‑1 (MQ=255)                                                                                              tCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCgg                                                    >  1:372734/1‑151 (MQ=255)                                                                                              tCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCgg                                                    >  2:677193/1‑151 (MQ=255)                                                                                              tCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCgg                                                    >  2:220995/1‑151 (MQ=255)                                                                                               cACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCgg                                                    <  1:185208/150‑1 (MQ=255)                                                                                                    cgaATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCa                                               >  1:450660/1‑150 (MQ=255)                                                                                                       aTAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGTCGGGAAACACTTACCGGGGCCAATGt                                           <  2:724210/151‑1 (MQ=255)                                                                                                       aTAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGt                                           <  1:58986/151‑1 (MQ=255)                                                                                                        tAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGtt                                          >  2:870055/1‑151 (MQ=255)                                                                                                         aaCCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGtt                                          <  1:953524/150‑1 (MQ=255)                                                                                                             tGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAgg                                     >  2:66458/1‑151 (MQ=255)                                                                                                              gACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAggg                                    >  1:744147/1‑151 (MQ=255)                                                                                                                 cGTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACa                                 <  2:707787/151‑1 (MQ=255)                                                                                                                  gTAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACa                                 >  1:46633/1‑150 (MQ=255)                                                                                                                   tAACAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCTCTTATAAAACGTCGACGGGAAACTCTTACCGGGGCCAATGTTATATGGActa                                <  2:276076/150‑3 (MQ=255)                                                                                                                    aaCAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACaaa                               >  2:542549/1‑150 (MQ=255)                                                                                                                    aaCAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACaaa                               <  1:670295/150‑1 (MQ=255)                                                                                                                     aCAGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAAt                             <  2:127369/151‑1 (MQ=255)                                                                                                                       aGTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGa                           <  2:911183/151‑1 (MQ=255)                                                                                                                        gTTATTTGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGAc                          <  2:780575/151‑1 (MQ=255)                                                                                                                              tGAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGAt                     <  1:335763/150‑1 (MQ=255)                                                                                                                               gAACGAATGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATa                    <  2:854488/150‑1 (MQ=255)                                                                                                                                    aaTGCGGCAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCt              >  1:673582/1‑151 (MQ=255)                                                                                                                                           cAATCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTtgtg       <  2:890482/151‑1 (MQ=255)                                                                                                                                            aaTCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTTGTgg      <  2:158802/151‑1 (MQ=255)                                                                                                                                              tCGCGCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTTGTGGtt    >  2:180030/1‑151 (MQ=255)                                                                                                                                               cgcgCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTTGTGGTTc   <  1:68459/151‑1 (MQ=255)                                                                                                                                               cgcgCCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTTGTGGTTc   <  1:870135/151‑1 (MQ=255)                                                                                                                                                gcgcCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTTGTGGTTCg  <  1:780532/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    ATTTCCGACGATAAATCCTGACAACACGAGAAAAATCACCATGAATATCTGATTCACTGTTCTATCCCCCAGCAGGTTAATGCGCTCTCCTGGTCACGACGAATAACCTGACCGTAACAGTTATTTGAACGAATGCGGCAATCGCGTCCGCCATCCTTAATCCACCAGCGAATCGCTTCGCATGCACCTTTACGATCACCAGCATTCAGCCGCTTATAAAACGTCGACGGGAAACACTTACCGGGGCCAATGTTATAGGGACAAAATGACGCGATACCCGCTTTTTGTGGTTCG  >  NC_000913/1639471‑1639764

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