Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,403,010	C→T	R325C (CGT→TGT)	hflK →	regulator of FtsH protease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,403,010	0	C	T	96.1%	160.8 / NA	52	R325C (CGT→TGT)	hflK	regulator of FtsH protease
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (25/24);  minor base A (2/0);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 4.92e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCG  >  NC_000913/4402765‑4403256                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       cgccgGAAGAAGTAAAAGAGGCGTTTGACGATGCGATTGCCGCGCGTGAACACGCACAGCAATACATTCGTGCAGCAGAAGCGTCTACCAACGAAGTTCAGCCGCGTGCGAACGGTAAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACCAGGCCCAGACCATCCTGGAAGCTCAGGGAGAAGTGGCGCGCTTTGCTAACCTTCTGCCGGAATATAAAGCCGCGCCGGCAATTACTCGCGAGTGTCTg                                                                                                                                                                                                                                                   >  2:380271/1‑251 (MQ=255)          aagTAAAAGCGTCGTTTGACGATGCTATTGACGCGCGTTAAAACGAACAGCAATACATTCGTGAAGCTGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGAACATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAATGCCGCGCCGTAAATTACTCGAGAGTGTCTGTATATCgag                                                                                                                                                                                                                                          <  2:351473/251‑1 (MQ=255)            gTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCTAGTGTCTGTATATCGAGAc                                                                                                                                                                                                                                        >  1:140890/1‑251 (MQ=255)              aaaaGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGCATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACg                                                                                                                                                                                                                                       >  1:77707/1‑250 (MQ=255)                    tgcgTTTGACGATGCGATTGCCGCGCGTGAAAAAGACCAGCAATACATTCGTGAAGCAGACGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGTCGCGCTTTGCTAAACTTCGGCCGGCATATAAAGCCGCGCCGGAACTTACTCGCGAGGGTCTGTATATTGAGACGATggca                                                                                                                                                                                                                                 >  2:82279/2‑248 (MQ=255)                       gTTTGAAGATGCGATTGCCGCGCGTGAAACCGACCAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGCGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGCAATTACTCGCGAGTGTCTGTATATCGAGACGATGGaaaa                                                                                                                                                                                                                               >  2:319562/1‑249 (MQ=255)                          tGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAgtg                                                                                                                                                                                                                           <  2:401813/250‑1 (MQ=255)                          tGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAgtg                                                                                                                                                                                                                           <  1:193214/250‑1 (MQ=255)                           gACGATGCGATTGCCGCTCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGtt                                                                                                                                                                                                                         <  2:120585/251‑1 (MQ=255)                                         ctcgcgTGACAACGAACAGCAATAAATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGATCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATCAACCCGCGCCGGAACTTACTCGCGAGAGTCTGTATATCGAGACGATGGAAAAAGTGTTAGGTAACACCCga                                                                                                                                                                                                            >  2:215250/3‑249 (MQ=255)                                             cgTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGAGTTGGGTAACACCCGCAAAGt                                                                                                                                                                                                       >  2:345144/1‑251 (MQ=255)                                                 aaaaCGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGAGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGAGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGAAGATGGAAAAAGTGGTGGGTAACACCCGCAAAGTGCTg                                                                                                                                                                                                   >  2:9979/1‑251 (MQ=255)                                                       aaCAGCAATACCTTCGTGAAGCAGAAGCGAATACCAACGAAGTTCAGCCACGTGCGAACGGTCAGGCGCAACGCATCCTCGAAGAGGCGCGTGCGTACAAAGACCAGACCATCCTGGAAGCTCAGGGTGACGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCAGCGCCGGAAATTACTAGCGCGTGTCTGAATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGGTAAc                                                                                                                                                                                             >  2:355444/1‑251 (MQ=255)                                                          aGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGCGACGATGGAAAAAGTGTTGGGTACCACCCGCAAAGTGCTGGTTAACGAt                                                                                                                                                                                          >  2:323711/1‑251 (MQ=255)                                                                 cATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACTATAAAGGTg                                                                                                                                                                                   <  1:125727/251‑1 (MQ=255)                                                                        gAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACCTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCt                                                                                                                                                                            <  1:388076/251‑1 (MQ=255)                                                                             aGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATgg                                                                                                                                                                       >  1:200428/1‑251 (MQ=255)                                                                              gAAGCGCATACCAACGAAGTTCAGCCGCGTGCGAACGGACAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGCCGATGGAAAAAGTGTTGGGTAACACACGAAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGt                                                                                                                                                                      >  2:163047/1‑251 (MQ=255)                                                                                aGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTc                                                                                                                                                                    >  1:141409/1‑251 (MQ=255)                                                                                     agacccacgaaatgaaGCCGCGTGCGAACGGTAAGGCGCGACGTATCCTCGAAGAGGCGCGTGCCTACAAGGCCCAATCCATCCTGGAAGGTCAGGGTGAAGTGGCGCGCTTTGCTAAGCTTCGGCCGGAATATAAAGCCGCGCATGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGTGTAACACCCGCAAAGTGCTGGTTAACGATAAAAGTGGCAACCTGATTGTTCTGcc                                                                                                                                                                <  2:5503/235‑1 (MQ=255)                                                                                        ccaaccacgTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCATCGAAGCGGCGCGTGCGTACAAGGCCCAGACCCTCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTCACCTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGATGGTTCACGATCAAGGGGGCAACCTGATGGTTCTGCCGtt                                                                                                                                                             >  2:401441/9‑250 (MQ=255)                                                                                        ccAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTAGCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAAGATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTa                                                                                                                                                            <  2:342271/251‑1 (MQ=255)                                                                                         cAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGtt                                                                                                                                                             >  1:30211/1‑249 (MQ=255)                                                                                         cAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGCGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGACTGTATATCGAGACGATGGAAAAAGAGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGtt                                                                                                                                                             >  2:92957/1‑249 (MQ=255)                                                                                           aCGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGAc                                                                                                                                                         <  1:243395/251‑1 (MQ=255)                                                                                              aaGTTCAGCCGCGTGCGAACGGTCAGGCGCAACTGATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAg                                                                                                                                                      <  2:323972/251‑1 (MQ=255)                                                                                                                 cGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAAATTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGAGTCTGTCTATCGAGACGATGGAAAAAGAGTTGGGTAACACCCGCAAAGTGCTGGGTAACGATAAAGGTGGCCACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACg                                                                                                                                   >  1:231085/1‑251 (MQ=255)                                                                                                                  ggTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGc                                                                                                                                  >  2:99075/1‑251 (MQ=255)                                                                                                                   gTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGCTAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGcc                                                                                                                                 >  1:244112/1‑251 (MQ=255)                                                                                                                     cAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAATGTGGTAACGcccc                                                                                                                               <  1:323711/251‑1 (MQ=255)                                                                                                                         cgcAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTgcg                                                                                                                           <  1:31577/251‑1 (MQ=255)                                                                                                                            aaCGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTgcggcg                                                                                                                        <  1:92957/251‑1 (MQ=255)                                                                                                                                       gAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCggta                                                                                                              >  2:190469/1‑247 (MQ=255)                                                                                                                                           agGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGt                                                                                                         >  1:199044/1‑251 (MQ=255)                                                                                                                                            ggCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTg                                                                                                        <  1:196806/251‑1 (MQ=255)                                                                                                                                                cgtgcgtACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAg                                                                                                    <  1:259783/251‑1 (MQ=255)                                                                                                                                                    cgtACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGTCGCGCTTTTCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTGCTCGCGAGTGTCTGTATATCGAGACGATGGAAAAATGGTTGGGTAACACCCGGAATGTGCTGGTTAACGATAAAGTTGGCAGCCGGATGGTTCTGTTGTTATACCAGGTGCATAAAGGTGTGAACTCCCCTGCGGCTAAGAGCCATAACGGGTCCAGCAAt                                                                                                <  2:203570/251‑1 (MQ=255)                                                                                                                                                       aCAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATctg                                                                                             >  1:355766/1‑251 (MQ=255)                                                                                                                                                                       ccTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCa                                                                              <  1:345144/250‑1 (MQ=255)                                                                                                                                                                         tGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGACGGAATATAAAGCCGCGCCGGAAATTACTCGAGAGTGTCTGTATAGCGAGACGATGGAAAAAGTGTTGGGTAACCCCCGCAAAGTGATGGTTAACGATAAAGGGGGCAACCTGATGATTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAAAGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATATGCTGCGTCTGCCGCCAGcc                                                                           >  1:393127/1‑251 (MQ=255)                                                                                                                                                                               gtCAGGGGGAAGGGGGGCGCGTTGCTAAGCGTCTGCCGGAATATAAAGCCGCGGCGGAAATTACGCGCGAGTGTCTGTATAGCGAGACGATGGAAAAAGGGTGGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGGTGGTTCTGCCGTTAGACCAGAGGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCGGCGTCGGCCGCCAGCCACTtcc                                                                     <  2:355639/250‑1 (MQ=255)                                                                                                                                                                                        aaGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCCGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAAcc                                                            <  1:380271/251‑1 (MQ=255)                                                                                                                                                                                           tGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAAc                                                                                                            >  2:81024/1‑200 (MQ=255)                                                                                                                                                                                           tGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAAc                                                                                                            <  1:81024/200‑1 (MQ=255)                                                                                                                                                                                                       cTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGAAAGCAac                                             >  1:325481/1‑251 (MQ=255)                                                                                                                                                                                                             ttCTGCCGGAATATAAAGCAGCGCCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGTTACAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCACCAcggcg                                       >  1:173183/1‑248 (MQ=255)                                                                                                                                                                                                                           aaaGCCGCGCCGGAAATTACTCGCGAGTGTCGGTATATCGAGACGAGGGAAAAAGTGTTGGGTAACACCCGCAAAGGGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGGGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCCAGCAACACGTCGTCCACCAGTCAggg                         <  1:242152/251‑1 (MQ=255)                                                                                                                                                                                                                                 gcgcCGGAAATTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGGGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGatat                   <  2:329800/251‑1 (MQ=255)                                                                                                                                                                                                                                        aaaTTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGACAGCCTCTTCCTCCACCACCAGTGGAGCAAGCAACACGTCGTCCACCGGTCAGGGCGATATTATGGAc            >  1:249005/1‑251 (MQ=255)                                                                                                                                                                                                                                         aaTTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGAcc           >  1:71943/1‑251 (MQ=255)                                                                                                                                                                                                                                          aTTACTCGCGAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCa          <  1:252707/251‑1 (MQ=255)                                                                                                                                                                                                                                                  cgAGTGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGGGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCcgcg  <  2:26762/251‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                       CGCCGGAAGAAGTAAAAGCGGCGTTTGACGATGCGATTGCCGCGCGTGAAAACGAACAGCAATACATTCGTGAAGCAGAAGCGTATACCAACGAAGTTCAGCCGCGTGCGAACGGTCAGGCGCAACGTATCCTCGAAGAGGCGCGTGCGTACAAGGCCCAGACCATCCTGGAAGCTCAGGGTGAAGTGGCGCGCTTTGCTAAACTTCTGCCGGAATATAAAGCCGCGCCGGAAATTACTCGCGAGCGTCTGTATATCGAGACGATGGAAAAAGTGTTGGGTAACACCCGCAAAGTGCTGGTTAACGATAAAGGTGGCAACCTGATGGTTCTGCCGTTAGACCAGATGCTGAAAGGTGGTAACGCCCCTGCGGCGAAGAGCGATAACGGTGCCAGCAATCTGCTGCGTCTGCCGCCAGCCTCTTCCTCCACAACCAGTGGAGCAAGCAACACGTCGTCCACCAGTCAGGGCGATATTATGGACCAACGCCGCG  >  NC_000913/4402765‑4403256

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