Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	268,643	C→T	A106A (GCG→GCA)	ykfA ←	CP4‑6 prophage; putative GTP‑binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	268,643	0	C	T	100.0%	247.0 / NA	63	A106A (GCG→GCA)	ykfA	CP4‑6 prophage; putative GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (27/36);  total (27/36)

CGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGCGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCGACTTGCAGGCGAAAGCGCAGTGGCTCACGTGTACAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGC  >  NC_000913/268496‑268789                                                                                                                                                    |                                                                                                                                                   cGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAgtgc                                                                                                                                                  >  1:782958/1‑150 (MQ=255)  gACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTgcag                                                                                                                                                <  2:532332/151‑1 (MQ=255)   aCTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTgcagc                                                                                                                                               <  1:843540/151‑1 (MQ=255)    cTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTgcagc                                                                                                                                               >  1:870080/1‑150 (MQ=255)      ggCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCAt                                                                                                                                             >  1:893829/1‑150 (MQ=255)         tGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTc                                                                                                                                         >  1:445295/1‑151 (MQ=255)           aTAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCgg                                                                                                                                       >  1:458153/1‑151 (MQ=255)              acaaacaGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTAt                                                                                                                                    >  1:808391/1‑151 (MQ=255)               caaacaGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATc                                                                                                                                   <  1:376562/151‑1 (MQ=255)                 aacaGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATcgc                                                                                                                                 <  1:240100/151‑1 (MQ=255)                    aGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAg                                                                                                                              <  1:931513/151‑1 (MQ=255)                       aCCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAgcgc                                                                                                                           >  1:859038/1‑151 (MQ=255)                         cTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGcc                                                                                                                          >  2:319893/1‑150 (MQ=255)                          ttATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCAc                                                                                                                        >  2:588574/1‑151 (MQ=255)                            aTGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACt                                                                                                                       <  1:233003/150‑1 (MQ=255)                                cGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGc                                                                                                                  <  2:197887/151‑1 (MQ=255)                                    ataCCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCAc                                                                                                              <  1:38306/151‑1 (MQ=255)                                          ctcCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCggg                                                                                                        >  1:729985/1‑151 (MQ=255)                                            ccAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCa                                                                                                      <  2:125671/151‑1 (MQ=255)                                            ccAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCa                                                                                                      <  2:604055/151‑1 (MQ=255)                                             cAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGc                                                                                                       >  2:445483/1‑149 (MQ=255)                                             cAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAg                                                                                                     <  2:1033238/151‑1 (MQ=255)                                               aTCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGc                                                                                                       >  1:875398/1‑147 (MQ=255)                                               aTCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGAt                                                                                                   <  2:342720/151‑1 (MQ=255)                                                 cACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATc                                                                                                  <  2:323236/150‑1 (MQ=255)                                                  aCCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATcc                                                                                                 >  1:668119/1‑150 (MQ=255)                                                  aCCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATcc                                                                                                 >  1:591721/1‑150 (MQ=255)                                                   ccTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCAt                                                                                               <  1:692373/151‑1 (MQ=255)                                                      gatgatAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCa                                                                                             <  1:870370/150‑1 (MQ=255)                                                      gatgatAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAg                                                                                            <  2:976809/151‑1 (MQ=255)                                                       atgatAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGc                                                                                           >  2:862135/1‑151 (MQ=255)                                                       atgatAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACCTTCGCCCACGCCGGGCAGATCCATCAGc                                                                                           >  1:1047720/1‑151 (MQ=255)                                                         gatAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGt                                                                                         <  1:212887/151‑1 (MQ=255)                                                             aaaaTGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTtat                                                                                      >  2:3318/1‑150 (MQ=255)                                                                aTGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTtatat                                                                                    >  2:1065719/1‑149 (MQ=255)                                                                   ttCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCgg                                                                               >  2:312072/1‑151 (MQ=255)                                                                      aTCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCg                                                                            <  2:57826/151‑1 (MQ=255)                                                                       tCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGc                                                                           >  1:465727/1‑151 (MQ=255)                                                                        ccACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGTGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGc                                                                           <  1:711824/150‑1 (MQ=255)                                                                         cACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCg                                                                         >  2:25014/1‑151 (MQ=255)                                                                           ccGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCGAc                                                                       >  2:330042/1‑151 (MQ=255)                                                                           ccGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCGAc                                                                       <  1:830042/151‑1 (MQ=255)                                                                             gTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCGACtt                                                                     >  1:730881/1‑151 (MQ=255)                                                                                 gcgcCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCGACTTGCAg                                                                 <  2:412779/151‑1 (MQ=255)                                                                                    cccGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCGACTTGCAGGCg                                                              <  2:507669/151‑1 (MQ=255)                                                                                          cATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCGACTTGCAGGCGAAAgc                                                         <  2:465716/150‑1 (MQ=255)                                                                                                       cAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCGACTTGCAGGCGAAAGCGCAGTGGCTCACg                                            <  1:123242/150‑1 (MQ=255)                                                                                                         gCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCGACTTGCAGGCGAAAGCGCAGTGGCTCACgtgt                                         <  2:468841/151‑1 (MQ=255)                                                                                                           cacaGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCGACTTGCAGGCGAAAGCGCAGTGGCTCACGTGTAc                                       >  2:771632/1‑151 (MQ=255)                                                                                                            acaGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCGACTTGCAGGCGAAAGCGCAGTGGCTCACGTGTACa                                      <  2:101619/151‑1 (MQ=255)                                                                                                             caGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCGACTTGCAGGCGAAAGCGCAGTGGCTCACGTGTAc                                       >  1:133591/1‑149 (MQ=255)                                                                                                              aGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCGACTTGCAGGCGAAAGCGCAGTGGCTCACGTGTACAgg                                    <  1:418706/151‑1 (MQ=255)                                                                                                                      gTCGAGGCGAGGAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCGACTTGCAGGCGAAAGCGCAGTGGCTCACGTGTACAGGCCGCcac                             >  2:334259/1‑150 (MQ=255)                                                                                                                        cgaggcgaggAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCGACTTGCAGGCGAAAGCGCAGTGGCTCACGTGTACAGGCCGCCACATc                          <  2:947209/151‑1 (MQ=255)                                                                                                                         gaggcgaggAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTGCGCCCACGCCGGGCAGATCCAGCAGCGTTATATAGCGGTCGCCGACTTGCAGGCGAAAGCGCAGTGGCTCACGTGTACAGGCCGCCACATCg                         <  2:87314/151‑1 (MQ=255)                                                                                                                           ggcgaggAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCGACTTGCAGGCGAAAGCGCAGTGGCTCACGTGTACAGGCCGCCACATCGCt                       <  1:918466/151‑1 (MQ=255)                                                                                                                                ggAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCGACTTGCAGGCGAAAGCGCAGTGGCTCACGTGTACAGGCCGCCACATCGCTGACCg                  <  1:59064/151‑1 (MQ=255)                                                                                                                                 gAAGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCGACTTGCAGGCGAAAGCGCAGTGGCTCACGTGTACAGGCCGCCACATCGCTGACCgg                 <  2:147511/151‑1 (MQ=255)                                                                                                                                   aGCTGTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCGACTTGCAGGCGAAAGCGCAGTGGCTCACGTGTACAGGCCGCCACATCGCTGACCGGTg               <  1:1047629/151‑1 (MQ=255)                                                                                                                                       gTTCGCGGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCGACTTGCAGGCGAAAGCGCAGTGGCTCACGTGTACAGGCCGCCACATCGCTGACCGGTGat             <  2:805236/149‑1 (MQ=255)                                                                                                                                           gcgGTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCGACTTGCAGGCGAAAGCGCAGTGGCTCACGTGTACAGGCCGCCACATCGCTGACCGGTGATATAtct       <  1:230315/151‑1 (MQ=255)                                                                                                                                             ggTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCGACTTGCAGGCGAAAGCGCAGTGGCTCACGTGTACAGGCCGCCACATCGCTGACCGGTGATATATCTcc     >  1:859834/1‑151 (MQ=255)                                                                                                                                              gTACAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCGACTTGCAGGCGAAAGCGCAGTGGCTCACGTGTACAGGCCGCCACATCGCTGACCGGTGATATATCTCCg    <  2:423071/151‑1 (MQ=255)                                                                                                                                                aCAGTGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCGACTTGCAGGCGAAAGCGCAGTGGCTCACGTGTACAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGc  >  2:402846/1‑151 (MQ=255)                                                                                                                                                    |                                                                                                                                                   CGACTGGCTGATAACAAACAGCACCTTATGCCGGTATACCTCTCCAATCACCTGATGATAAAAATGTTCATCCACCGTCAGCGCCCGGTCATCGGCCTTAATCAGCCACAGTACCAGGTCGAGGCGAGGAAGCTGTTCGCGGTACAGCGCAGCATACTCGGTATCGCGAGCGCCACTTTCGCCCACGCCGGGCAGATCCATCAGCGTTATATAGCGGTCGCCGACTTGCAGGCGAAAGCGCAGTGGCTCACGTGTACAGGCCGCCACATCGCTGACCGGTGATATATCTCCGGC  >  NC_000913/268496‑268789

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