Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,938,441	Δ1 bp	coding (215/993 nt)	rbsR →	DNA‑binding transcriptional dual regulator RbsR

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,938,441	0	C	.	98.4%	289.2 / ‑2.4	63	coding (215/993 nt)	rbsR	DNA‑binding transcriptional dual regulator RbsR
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (33/29);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 4.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.29e-01

GACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCGACAGCG  >  NC_000913/3938197‑3938689                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         ggcGCATTTTTAGACAGGCAGAGGTGATGCTTGGCTACAGTGAAAGATGTTGCCAGCAGGGCTGGCGTGTCTACCTCAACAGTTTCTCACTTTATCAAGAAAGATCGCTTCGGCAGTGAAGAGATTCCCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCCCAACCTCAATCAAACACAGACCATTGGCATGTTGATCACTGCCAGTACCAA‑CCTTTcgaa                                                                                                                                                                                                                                                   <  2:105187/249‑4 (MQ=255)      aTTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTc                                                                                                                                                                                                                                                 >  1:37674/1‑248 (MQ=255)           tAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTg                                                                                                                                                                                                                                          <  1:214585/250‑1 (MQ=255)            aGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTgg                                                                                                                                                                                                                                         <  1:274999/250‑1 (MQ=255)              aCAGGCAAAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTG‑CGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCg                                                                                                                                                                                                                                     <  1:84627/251‑1 (MQ=255)                    agagGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCg                                                                                                                                                                                                                                >  1:146996/1‑251 (MQ=255)                        gTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTg                                                                                                                                                                                                                             >  1:229836/1‑250 (MQ=255)                             gCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGc                                                                                                                                                                                                                        >  1:278938/1‑250 (MQ=255)                               ttGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTAACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGAATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGCACGCAgc                                                                                                                                                                                                                     >  1:80366/1‑251 (MQ=255)                                tGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCAAGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCACTTACGCGCCATCAGCTCTGGCGCGTAGCCTCCAACTCAATCAAACACATACCCTTGGCCTGTTGATCACTGCCCGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAgct                                                                                                                                                                                                                    >  2:184223/1‑251 (MQ=255)                                         tGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGCTCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGGGCGTGGCGTTGAACGCAGCTGCTTCGAAc                                                                                                                                                                                                           >  2:37026/1‑251 (MQ=255)                                         tGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAAc                                                                                                                                                                                                           >  1:311506/1‑251 (MQ=255)                                          gAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAAcg                                                                                                                                                                                                          <  2:313136/251‑1 (MQ=255)                                             aGATGGTGCCCGCCTGGCGTGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCTCGTAGCCTCAAACTCAATCAAACACATACCATGGGCATGTGGATCACGGCCAGTACCAATC‑TTTCGATTCAGAACGGGTGCGTGGCGTTGAACGCAGCTGTTTCGAACGCgg                                                                                                                                                                                                       <  2:345208/251‑1 (MQ=255)                                               aTGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGtt                                                                                                                                                                                                     <  1:90515/251‑1 (MQ=255)                                                               gggcgTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTg                                                                                                                                                                                     >  1:279579/1‑251 (MQ=255)                                                               gggcgTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACGCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGGCCTTTg                                                                                                                                                                                     <  2:383533/251‑1 (MQ=255)                                                                             cAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTAACGCCAAAGTTGAAGCGGCGATTAACGAACTCAATTACGCGACATCAGCTCTGGCGCGCAGCCTCAAACTCAATCAAACACATAACATTGGCATGTTGATCACTGCCAGAACCAATC‑TTTCTATTCAGAACTGGTGCGTGACGTTGAACGCAGCTGCTTCGCACGCGGTTATAGTCTCGTCCTTTACAATACCGAAGGCg                                                                                                                                                                       >  1:401351/1‑251 (MQ=255)                                                                                        gccggtagcAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTGGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGc                                                                                                                                                              <  1:328148/241‑1 (MQ=255)                                                                                        cACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGAGGAGCAGCgg                                                                                                                                                            <  1:308588/251‑1 (MQ=255)                                                                                             tATCAATAAAGATCGCTTCGTCACTGAAGCGCTTAACGCCAAAGTTGAAGCGGCGATTAAAGAACTCACTTACGCGCCCTCAGCACTGACGCGTAGCCTCAAACTCACTCAAACACATACCATTCGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGCACTGGTGCGTGGCGTTGACCGCAGCTACTTCGCACGCGGTTATAGTCTCGTCCTTTGCAATACCGCAGGCGATGAACAGCGGATGa                                                                                                                                                        >  2:147053/1‑250 (MQ=255)                                                                                                 aaTAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGCACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGc                                                                                                                                                   >  1:101819/1‑251 (MQ=255)                                                                                                    aaaGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTACAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGCTCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAAt                                                                                                                                                >  1:37900/1‑251 (MQ=255)                                                                                                              tCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGCACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGc                                                                                                                                      >  1:46167/1‑251 (MQ=255)                                                                                                               cGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTTAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGc                                                                                                                                      <  1:302092/250‑1 (MQ=255)                                                                                                               cGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCt                                                                                                                                     <  2:106430/251‑1 (MQ=255)                                                                                                                   aggtaaGCGATTACCGCCAAAGTTGAAGCGGCGAGTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACAGACCATTGGCATGTTGATCACTGCCAGGACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCGGCGAACTCGGTTATAGTCTCGTCCGTTTCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATg                                                                                                                                 <  1:30091/247‑1 (MQ=255)                                                                                                                   aGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGTTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGGCTCGTc                                                                                                                                                                                          <  2:156617/194‑1 (MQ=255)                                                                                                                   aGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTc                                                                                                                                                                                          >  1:156617/1‑194 (MQ=255)                                                                                                                   aGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATg                                                                                                                                 <  1:37026/251‑1 (MQ=255)                                                                                                                      gaagcgaaaaacgacAAAGTTGAAGCGGCGAATAAAGAACAAAATTACGAGCCATCAGCTCTGGCGCGTAGACTAAAACACAATCAAACACATACCATTGGCATGTAGAACACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTAGAGTCTCGTCCTTTGCAATACAGAAGGCGAAGAACAGCGGATGAATCGCAATGTGGAAACGCTGATGCa                                                                                                                               >  2:279295/15‑250 (MQ=255)                                                                                                                       aaGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCaa                                                                                                                              <  1:180895/250‑1 (MQ=255)                                                                                                                        aGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCaaa                                                                                                                             >  1:88318/1‑250 (MQ=255)                                                                                                                           gATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGCACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGTTGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCCAAAAcg                                                                                                                         >  2:299946/1‑251 (MQ=255)                                                                                                                              tACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAAcgcg                                                                                                                       <  2:136318/250‑1 (MQ=255)                                                                                                                                ccGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGGGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCGTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTg                                                                                                                    <  1:233182/251‑1 (MQ=255)                                                                                                                                            aagcggccaatcaaGAACTCAATTAAGCGCCATCAGCTCTGCCGCGTAGCCTCAAACTCAATCAAACACACACCATTGGCATGATGATCCCTGCCAGTACCAATC‑TTTCTATTCAGACATGGTGCGTAGCGTTGAACGCAGCTGCTTAGCACGCGGTTATAGCCTCCTCCTTTGCAATACCGAAGGCGATGTACAGCGGATGAATCGCAATCTGGCAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGt                                                                                                        >  2:164075/13‑251 (MQ=255)                                                                                                                                                    aTTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGc                                                                                                <  1:147060/251‑1 (MQ=255)                                                                                                                                                    aTTAAAGAACTACAGTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTAAGAACTGGTGCGTAGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGCTGAATCGCAAGCTGGAAACGCTGATGCACACACGCGTTGATGGCTTGCTGATACtgtg                                                                                                 >  2:148369/1‑250 (MQ=255)                                                                                                                                                         aGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCg                                                                                            >  1:64662/1‑250 (MQ=255)                                                                                                                                                            aCTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAAc                                                                                        >  1:110750/1‑251 (MQ=255)                                                                                                                                                             cTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACg                                                                                       >  1:17479/1‑251 (MQ=255)                                                                                                                                                                     cgcgCCATCAGCTCTGGCGAGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGAGGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAAcc                                                                               <  2:9043/251‑1 (MQ=255)                                                                                                                                                                         ccATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCGATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTGGATGGCTTGCTGTTACTGTGCACCTATACGCATCTACCTTcg                                                                           <  2:267261/251‑1 (MQ=255)                                                                                                                                                                                    ggCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATc                                                                  <  1:69794/249‑1 (MQ=255)                                                                                                                                                                                        cgTAGCCTCAAACGCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCaa                                                            <  2:278938/251‑1 (MQ=255)                                                                                                                                                                                          tAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACg                                                          >  1:43526/1‑251 (MQ=255)                                                                                                                                                                                                  aaCTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCg                                                   >  1:25729/1‑250 (MQ=255)                                                                                                                                                                                                     tcaatcaaACACATACCCTTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATTCGAc                                                 >  2:132241/1‑249 (MQ=255)                                                                                                                                                                                                      caatcaaACACATACCATTGGCATGTTGATCCCTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCACCGTTATCCGACAGt                                              >  2:367521/1‑251 (MQ=255)                                                                                                                                                                                                      caatcaaACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGt                                              <  1:23656/251‑1 (MQ=255)                                                                                                                                                                                                           aaaCACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTg                                                                                                         <  1:113971/187‑1 (MQ=255)                                                                                                                                                                                                           aaaCACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTg                                                                                                         >  2:113971/1‑187 (MQ=255)                                                                                                                                                                                                            aaCACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTAc                                        >  2:141526/1‑251 (MQ=255)                                                                                                                                                                                                             acacATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACt                                       >  2:295311/1‑251 (MQ=255)                                                                                                                                                                                                                             tGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACt                          >  1:351619/1‑248 (MQ=255)                                                                                                                                                                                                                             tGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACt                          <  2:351619/248‑1 (MQ=255)                                                                                                                                                                                                                                tGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAAGCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCt                     <  2:127126/250‑1 (MQ=255)                                                                                                                                                                                                                                     aCTGCCAGTACCAATC‑TTTCTATTCAGAACTGGAGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGCACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGCAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTc               >  2:364954/1‑251 (MQ=255)                                                                                                                                                                                                                                     aCTGCCAGTACCAATC‑GTGCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTc               <  2:284333/251‑1 (MQ=255)                                                                                                                                                                                                                                        gCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCg              <  1:184306/249‑1 (MQ=255)                                                                                                                                                                                                                                           aGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGTCTCCGTTCGATGGc         >  1:286229/1‑251 (MQ=255)                                                                                                                                                                                                                                            gTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCg        <  2:37674/251‑1 (MQ=255)                                                                                                                                                                                                                                                 caatcTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCGACAg    <  1:260440/247‑1 (MQ=255)                                                                                                                                                                                                                                                   atcTTTCTATTCAGAAGGGGTGCGTGGCGTTGCACGCAGCTGCTTCGTACGCGGTTATAGTCTCGTCCTTTGAAATACCGACGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGCTGCTGGCTTGCTGTTACTGTGCACCGAAACGCAACAACCTTCGCGTGAAATCATGCCACGTTCTCCGACCGCGCCTACTGTGATGATGGAATGGGCGCCGTTCGTTGGCGCCAGCg  >  2:247528/3‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                         GACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATGAATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCGACAGCG  >  NC_000913/3938197‑3938689

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