Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,377,240	T→G	T61P (ACC→CCC)	sspA ←	stringent starvation protein A

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,377,240	0	T	G	98.5%	200.2 / ‑7.6	66	T61P (ACC→CCC)	sspA	stringent starvation protein A
Reads supporting (aligned to +/- strand):  ref base T (1/0);  new base G (19/46);  total (20/46)
Fisher's exact test for biased strand distribution p-value = 3.03e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGACCGGCGCAATCGCCAGCAGTTCTTCGCGCAGTTGCTTACGTGCGGCATCTGCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTG  >  NC_000913/3376994‑3377485                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       aGACCGGCGCAATCGTCAGCAGTTCTTCGCGCAGTTGCTTACGTGCGGCATGTGCTTCATAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTTTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTTGTGGAGGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGa                                                                                                                                                                                                                                                    <  1:237193/251‑1 (MQ=255)  gACCGGCGCAATCGCCAGCAGTTCTTCGCGCAGTTGCTTACGTGCGGCATCTGCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGaa                                                                                                                                                                                                                                                   <  2:94212/251‑1 (MQ=255)         gcAATCGCCAGCAGTTCTTCGCGAAGTTGATTACGTGCGGCACCTACTTCAGACGCTGAACCGTTGCTGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCAGCCAGACATTACATAATGATGAGGAATTCCCTCAGGGACAGCTCACGAACCACCACGGTCGGCACGCTCTg                                                                                                                                                                                                                                            >  2:91168/1‑251 (MQ=255)             tCGCCAGCAGTTCTTCGCGCAGTTGCTTACGTGCGGCATCTGCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGGCAGCTCCCGATCCACCAAGGGCGGAACGCTCTAATTc                                                                                                                                                                                                                                        >  1:120588/1‑251 (MQ=255)             acGCCAGCAGTTCTTCGCGCAGTTGCTTACGTGCGGCAGCTGCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCAGAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGAttt                                                                                                                                                                                                                                        <  2:23748/250‑2 (MQ=255)              cGCCAGCAGTTCTTCGCGCAGTTGCTTACGTGCGGGATCTGCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTc                                                                                                                                                                                                                                        <  1:102910/250‑1 (MQ=255)                  agcagTTCTTCGCGCAGTTGCTTACGTGCGGCATCTGCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGtt                                                                                                                                                                                                                                   <  1:444856/251‑1 (MQ=255)                             gcgcAGTTGCTTACGTGCGGCATCTGCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGGGATTCCCCCAGGGGCAGCTCACGATCCACCAGGGGCGGAACGCTCTTATTCGGGGTGAGGTGAATc                                                                                                                                                                                                                         >  2:147035/1‑250 (MQ=255)                               gcAGTTGCTTACGTGCGGCATCTGCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTTTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGTCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATGCCCACAGGGTCAGCTCACGAGCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGa                                                                                                                                                                                                                      <  2:320501/251‑1 (MQ=255)                                    tGCTTACGTGCGGCATCTGCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCCTCCAGATATTCCATAATGATGCGCGATTCCCACAGGGTCAGCTCCCGATCCACCAGGGGCGGAACGCTCCGGTTCCGGTTGAGGTCAATCAGATCCTg                                                                                                                                                                                                                 >  2:220074/1‑251 (MQ=255)                                    tGCTTACGTGCGGCATCTGCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCt                                                                                                                                                                                                                  <  1:115239/250‑1 (MQ=255)                                     gCTTACGTGCGGCATCTGCTTCAGAAGCTGAACCGTTGATGTTTGTGTTCATCAGCGTGTCCCAGTCTTTTTCGATGCGATTCATGTACAGACGGGTTTCACCGCTAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATTATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGa                                                                                                                                                                                                                <  1:448919/251‑1 (MQ=255)                                     gCTTACGTGCGGCATCTGCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTg                                                                                                                                                                                                                 <  1:91266/250‑1 (MQ=255)                                       ttACGTGCGGCATCTGCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATATACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGGTGAGGTCAATCAGATCCTGAgg                                                                                                                                                                                                              >  1:293932/1‑251 (MQ=255)                                        tACGTGCGGCATCTGCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAgg                                                                                                                                                                                                              <  1:140547/250‑1 (MQ=255)                                          cGTGCGGCATCTGCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCCCAGGGTCAGCTCACGATCCACCAGGGGCGGCACGCTCTGATTCGGGTTGAGGTCACTCAGATCCTGAGGCgg                                                                                                                                                                                                           >  1:255701/1‑251 (MQ=255)                                           gTGCGGCATCTGCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTTCCACAAGGTCAGCTCCCGATCCACCAGGGGCGGGACGCTCTGATTCGGGTTGAGGGCAATCCGATCCTGAGGCggg                                                                                                                                                                                                          >  2:361185/1‑250 (MQ=255)                                             gCGGCATCTGCTTCAGAAGCTGAACCGTTGATGATGGTTTTCATCGGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACTGCTTTCTCCGCGAGCTACCGTGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGAtt                                                                                                                                                                                                        <  2:373780/251‑1 (MQ=255)                                             gCGGCATCTGCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGAtt                                                                                                                                                                                                        <  1:173639/251‑1 (MQ=255)                                               ggCATCTTCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTGTCGATGCGATGCATGTACAGACGGCTTGCGCCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGAGAGCCTCAGGGTCAGCGCGCGAACCTCCAGGGGCGGATCGCTCTGATTCGGGGTGAGGTCAATCAGATCCTGAGCCGGATT‑Gt                                                                                                                                                                                                     <  2:210403/251‑1 (MQ=255)                                                      gCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCt                                                                                                                                                                                                  <  2:245323/247‑1 (MQ=255)                                                         ccAGAAGCTGAACCGTTGATGATGGTTGTCATCAGCGGGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCAc                                                                                                                                                                                           <  2:142558/250‑1 (MQ=255)                                                          cAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCAc                                                                                                                                                                                           <  1:70635/250‑1 (MQ=255)                                                           aGAAGCTGACCCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGCAAAGCTCCTCCCGATATTCCATAATGATGAGGGATTCCCCCAGGGTCAGCTCACGATCCCCCCGGGGCGGAACGCTCTGATTCGGGGTGAGGTCAATCACATCCCAGGGCGGGAT‑GTTCTTTTCCACgg                                                                                                                                                                                         >  2:228891/1‑250 (MQ=255)                                                              aGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGtt                                                                                                                                                                                      >  1:326340/1‑251 (MQ=255)                                                               gCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTc                                                                                                                                                                                     <  2:380302/251‑1 (MQ=255)                                                                 tGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCACCCAGATATTCCATAATGATGCGAGGTTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCCCTGATTCGGGGTGAGGTCAATCAGATCCTGAGGCGGATTTGTC‑TTTTCCCCGTGTTCgg                                                                                                                                                                                   >  2:171971/1‑250 (MQ=255)                                                                             gaaggTGTTCATTAGCGTGTACCAGTCTTTTGCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCGCATCCAGATATTCCATAAGGATGCGAGATTCCCACAGGGTCAGCTCACGAGCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGAGCCTGAGGAGGATT‑GTCCTTTTCCACGTGTTCGAGCTCGAAACTTa                                                                                                                                                                       <  2:428251/248‑1 (MQ=255)                                                                                      catcaGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACTCCTTtc                                                                                                                                                               <  2:19064/250‑1 (MQ=255)                                                                                          aGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTcagc                                                                                                                                                          <  1:216452/251‑1 (MQ=255)                                                                                               gtACCAGTCTTTTTCGATGCGAAGCGTGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGAGGGTAACGCTCCTCCAGGTATTCCAAAATGATGCGCGATTACCCCAGGGTCAGCTCACTATCCACCAGGGGAGGAACGCTCTGATTCGGGTTTGGGGCACTCAGATCCTGGTGCGGGTA‑GTCCTTTTCTACGTGTTCGATATCGAGACTTACACCTTTCTcacccactc                                                                                                                                                     >  2:413308/1‑248 (MQ=255)                                                                                                    aGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGTATT‑GTCCTTTTCCACTTGTTCGATCTCGAAACTGACACCTTTCTCAGCCAGCACAATg                                                                                                                                                <  2:399678/251‑1 (MQ=255)                                                                                                      tCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCCTAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCCGATCCTGAGGCGGATT‑Gt                                                                                                                                                                                                     >  1:227633/1‑196 (MQ=255)                                                                                                      tCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑Gt                                                                                                                                                                                                     <  2:227633/196‑1 (MQ=255)                                                                                                         ttttCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATGCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACTTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGa                                                                                                                                            <  1:407838/250‑1 (MQ=255)                                                                                                            tCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCt                                                                                                                                         <  2:321434/250‑1 (MQ=255)                                                                                                                 gcgatgcATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATgg                                                                                                                                    <  2:408192/250‑1 (MQ=255)                                                                                                                        aTGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑TTCCTTTTCCACGTGTTTGATCTCGAAACTTACACCTTTCTCAGCCAACACAATGCGGACCTGATGGCTATAggt                                                                                                                            >  1:74921/1‑249 (MQ=255)                                                                                                                             cAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCCTCCAGATATTCCATAATGATGCGAGATTCCCCCAGGGTCAGCTCACGATTCACCAGGGGCGGAACGCTCTGCTTCGGGTTGAGGTCAATCAGAACCTGAGGCGGATT‑GTCCTTTTCCCCGTGTTCGATCTCGAAACTTACACCTTTCTAAGCCATAACAATGCGGCCCCGATGGCTATAAATGTCa                                                                                                                        >  2:456547/1‑250 (MQ=255)                                                                                                                              aGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGAGGTCAGt                                                                                                                      <  2:213100/251‑1 (MQ=255)                                                                                                                                  ggCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGa                                                                                                                  <  1:302370/251‑1 (MQ=255)                                                                                                                                   gCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCCCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCCCGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCCGCACAATGCGGGCCTGATGGCTATAGATGTCAGTAGGAc                                                                                                                 >  2:163213/1‑251 (MQ=255)                                                                                                                                        cACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACAATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGa                                                                                                             <  2:30335/250‑1 (MQ=255)                                                                                                                                               aGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCg                                                                                                     <  1:35273/251‑1 (MQ=255)                                                                                                                                                     cGGGTAAACAGGCATCCGTGGCGGATGCGGGGAACACTCATCCAGATATTTCATAAAGGTGCGAGATTTCCACAGGGGCAGCTCACGATCCACCAGGGGCGGAACGCTCTGAATCGGGGTGAGGGCAATCAGATCCTGCGGCGGATT‑GTCCTTTTCCACGTGTTCGGTCTCGAAACTTACACCTTTCTCAGCCAGCACAATTCGGACCTGATGGCTATAGATGTCCGTAGGGCCGGAAAACAGCGTCAtt                                                                                                >  1:114182/1‑250 (MQ=255)                                                                                                                                                        gTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCg                                                                                            <  2:337539/251‑1 (MQ=255)                                                                                                                                                           aaCAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAAc                                                                                         <  2:147794/251‑1 (MQ=255)                                                                                                                                                                 cATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTttgt                                                                                   <  2:326340/251‑1 (MQ=255)                                                                                                                                                                      gTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCa                                                                              <  1:163213/251‑1 (MQ=255)                                                                                                                                                                           ggATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGAc                                                                         <  1:147035/251‑1 (MQ=255)                                                                                                                                                                           ggATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGAc                                                                         <  2:130062/251‑1 (MQ=255)                                                                                                                                                                                        gCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAAcc                                                            <  1:15035/251‑1 (MQ=255)                                                                                                                                                                                             tCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCta                                                                                                                                 >  1:116163/1‑177 (MQ=255)                                                                                                                                                                                             tCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCta                                                                                                                                 <  2:116163/177‑1 (MQ=255)                                                                                                                                                                                                    atTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGt                                                >  1:376050/1‑251 (MQ=255)                                                                                                                                                                                                        cATAATGAAGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGa                                            <  2:257819/251‑1 (MQ=255)                                                                                                                                                                                                              gatgCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGAt                                                                                                                            >  1:60386/1‑165 (MQ=255)                                                                                                                                                                                                              gatgCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGAt                                                                                                                            <  2:60386/165‑1 (MQ=255)                                                                                                                                                                                                                    agaTTCCCACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCt                                <  2:376050/251‑1 (MQ=255)                                                                                                                                                                                                                         cccACAGGGTCAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCa                            <  1:96230/250‑1 (MQ=255)                                                                                                                                                                                                                          ccACAGGGTCAGCTCTCGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGAGGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCAGTACCGAACGTTTGTTGGCAGCGACAGTCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGc                          <  1:288128/251‑1 (MQ=255)                                                                                                                                                                                                                                   cAGCTCACGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTg                 <  1:361185/251‑1 (MQ=255)                                                                                                                                                                                                                                    aGCTCACGATCCCCCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCGCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCATTAGGACCGGAAAACAGCGTCATGACCGAACGGTTGTTGGCAGCGACAGCCATGAAAACCTCCAGTTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTgg                <  2:167876/251‑1 (MQ=255)                                                                                                                                                                                                                                         aCGATCCACCAGGGGCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGt           >  1:146088/1‑251 (MQ=255)                                                                                                                                                                                                                                          cGATCCACCCGGGGCGGAACGCGCTGATTCGGGTTGAGGTCAATCAAATCCTGAGGCGGATT‑GTCCTTTACCCCGTGTTCGATCTCGAACCTTACACATTTCTCAACCAGCACAACGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCATCATTACCGACCGCTTGTTCGCAGCGCCAGCCACGAAAACCCCCAGGTATAGTCAGAACATTTACGGCTTCCAGCCACCAGATGGCCAGTc          <  2:453213/251‑1 (MQ=255)                                                                                                                                                                                                                                                  caGGGGCGGAACGCGCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGGGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTtg  <  1:451546/251‑1 (MQ=255)                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                       AGACCGGCGCAATCGCCAGCAGTTCTTCGCGCAGTTGCTTACGTGCGGCATCTGCTTCAGAAGCTGAACCGTTGATGATGGTGTTCATCAGCGTGTACCAGTCTTTTTCGATGCGATGCATGTACAGACGGCTTTCACCGCGAGCTACCGGGTAAACAGGCATCAGTGGCGGATGCGGGAAACGCTCATCCAGATATTCCATAATGATGCGAGATTCCCACAGGGTCAGCTCACGATCCACCAGGGTCGGAACGCTCTGATTCGGGTTGAGGTCAATCAGATCCTGAGGCGGATT‑GTCCTTTTCCACGTGTTCGATCTCGAAACTTACACCTTTCTCAGCCAGCACAATGCGGACCTGATGGCTATAGATGTCAGTAGGACCGGAAAACAGCGTCATTACCGAACGTTTGTTGGCAGCGACAGCCATGAAAACCTCCAGGTATAGTCAGAATTTTTACTGCTACCAGCCACCAGGTGGCCAGTCAGAAGTTG  >  NC_000913/3376994‑3377485

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