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.5%	308.0 / ‑2.1	68	coding (215/993 nt)	rbsR	DNA‑binding transcriptional dual regulator RbsR
Reads supporting (aligned to +/- strand):  ref base C (0/1);  new base . (33/34);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

GACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCGACA  >  NC_000913/3938197‑3938686                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       gACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGTTGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCCATGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTAt                                                                                                                                                                                                                                                 >  2:78965/1‑251 (MQ=255)     cATTTTTAGACAGGCAGAGGAGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTc                                                                                                                                                                                                                                               >  1:237247/1‑249 (MQ=255)                  gCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTgg                                                                                                                                                                                                                                >  2:80556/1‑251 (MQ=255)                    agagGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCg                                                                                                                                                                                                                              >  2:504830/1‑251 (MQ=255)                         tGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGa                                                                                                                                                                                                                          >  1:334841/1‑250 (MQ=255)                                 ggccaaaatgtgagattgtGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTGTCTCACGTTATCAATAAAGAGCGCTTAGTCAGTGAAGCGATTCCCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATCCCATTGGCATGTTGATCACTGCCAGTTCCA‑TTCTTTCTATGCAGAACTGGGGCGTGGCGTTGAACGCAgct                                                                                                                                                                                                                  <  2:494832/232‑1 (MQ=255)                                 ggCTACAATGAAAGATGTGGCCCGCCGGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGGTGATCACTGCCAGTACCAATC‑TTTCTATGCAGAACTGGTGCGTGGCGTTGAACGCAgct                                                                                                                                                                                                                  <  2:511155/250‑1 (MQ=255)                                     aCAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCAGCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGTCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCtt                                                                                                                                                                                                              <  1:408584/250‑1 (MQ=255)                                       aaTGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCg                                                                                                                                                                                                            >  1:175480/1‑250 (MQ=255)                                         cgAAAGATGTTGCCCGCCTGGCAGGCGTTTCTACCTCATCAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTTGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTAGAAc                                                                                                                                                                                                         >  2:22273/2‑251 (MQ=255)                                              gATGTTGCACGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGACAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACAATTGGCATGTTGATAACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGACCGCACCTGCTTCGAACACGGt                                                                                                                                                                                                    >  2:383172/1‑251 (MQ=255)                                                    gCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAg                                                                                                                                                                                               <  1:504830/250‑1 (MQ=255)                                                      ccGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATCAAGATCGCTTCGTCAGTGAAGCGATTACCGCCACAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTTGCGCGTAGCCTACAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTCCCAATC‑TTTCAATTCAGTACTGGTGCGTGGCGTTGAACGCAGCTGCTTAGAACGCGGTTATAGtcc                                                                                                                                                                                            >  2:133420/1‑250 (MQ=255)                                                       cGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGtctc                                                                                                                                                                                           <  2:175480/251‑1 (MQ=255)                                                       cGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTAAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTCATAGtctc                                                                                                                                                                                           >  1:84131/1‑251 (MQ=255)                                                           tGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGGGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTcc                                                                                                                                                                                       >  2:104136/1‑251 (MQ=255)                                                            ggcgggcgTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTcc                                                                                                                                                                                       <  1:75355/250‑1 (MQ=255)                                                            ggcgggcgTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTc                                                                                                                                                                                        <  1:114421/249‑1 (MQ=255)                                                            ggcgggcgTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGACTCGTCCt                                                                                                                                                                                      <  2:271844/251‑1 (MQ=255)                                                               gggcgTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTg                                                                                                                                                                                   <  1:389532/251‑1 (MQ=255)                                                                             cAACATTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGGGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGc                                                                                                                                                                      <  2:334841/250‑1 (MQ=255)                                                                              aaCAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCg                                                                                                                                                                     >  1:93351/1‑250 (MQ=255)                                                                                 aGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCGTCGAACACGGTTATAGTCTCCTCCTTTGCAATACCGAAGGCGATGa                                                                                                                                                                 <  2:88289/251‑1 (MQ=255)                                                                                         aCGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCgg                                                                                                                                                          >  1:255733/1‑250 (MQ=255)                                                                                                cAATAAAGATCGCTTCGACAGGGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCTTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCg                                                                                                                                                 <  2:569175/251‑1 (MQ=255)                                                                                                    aaaGATCGCTTCGTCAGTGAAGCGATTACCGCCACAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTAGCGCGTAGCATCAAACTCAATCAAACACATACCATTGGCATGTTGATCCCTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGCACGCGGTTATAGTCTCATCCTTTGCAATACCGAAGGCGATTAACAGCGGATG‑AATCGCAAt                                                                                                                                             >  2:430874/1‑251 (MQ=255)                                                                                                    aaaGAGCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAAGCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCaag                                                                                                                                             <  2:462598/251‑2 (MQ=255)                                                                                                       gATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGAAATCTg                                                                                                                                          >  1:542402/1‑251 (MQ=255)                                                                                                        aTCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTgg                                                                                                                                         >  1:297345/1‑251 (MQ=255)                                                                                                               cGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCt                                                                                                                                  >  1:136446/1‑251 (MQ=255)                                                                                                                gTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AAGCGCAATCTGGAAACGCTg                                                                                                                                 <  1:440093/251‑1 (MQ=255)                                                                                                                         gataTTACCGCCAAATTTGAAGCGGCTATTAAAGAACTCAATTACCCGCCATCATCTCTGGCGCGTAGCCTCAAACTCAACCAAACACATACCACTGGCATGTTGATCACTTCCAGTACCAATC‑TTTCTATTCATAACGGGTGCGTGGCGTTGAACGCAGCTTCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGG‑TGTAATCGCAATCTGGAAACGCTGATGCaaaaa                                                                                                                        <  1:479761/248‑1 (MQ=255)                                                                                                                          cGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACGGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAAc                                                                                                                       <  1:80556/251‑1 (MQ=255)                                                                                                                           gATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAAcg                                                                                                                      <  2:189479/251‑1 (MQ=255)                                                                                                                           gATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATAACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCACTACCGCAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGCTGCAAAAAcg                                                                                                                      >  1:547087/1‑251 (MQ=255)                                                                                                                            aTTACCGCCAACGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAAcgc                                                                                                                     >  1:114618/1‑251 (MQ=255)                                                                                                                                    cAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATgg                                                                                                             >  2:160373/1‑251 (MQ=255)                                                                                                                                           gAAGCTGCGATTAAAGAACACAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGGGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAACCGCGTTGATGGCTTGCTg                                                                                                      >  2:435590/1‑251 (MQ=255)                                                                                                                                               cggcgATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGAGCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTTCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTAc                                                                                                  <  1:264774/251‑1 (MQ=255)                                                                                                                                                ggcgATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACt                                                                                                 >  2:192581/1‑251 (MQ=255)                                                                                                                                                       aaaGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCAcc                                                                                          <  2:78688/251‑1 (MQ=255)                                                                                                                                                                    aCGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCaa                                                                              <  1:9122/250‑1 (MQ=255)                                                                                                                                                                    aCGCGCCATCAGCTCTGGCGCGCAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGATATAGTCTCGTCCTTTGCAATACCGCAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAACACGCGTTGTTGGCTTGCTGTTACTGTGCAACGAAACGCATCAAc                                                                             >  1:94791/1‑251 (MQ=255)                                                                                                                                                                     cgcgCCATCAGCTCTGGCGCGTAGCCTCAAACTCACTCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAAcc                                                                            >  2:351534/1‑251 (MQ=255)                                                                                                                                                                        gcCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTc                                                                         <  1:124923/251‑1 (MQ=255)                                                                                                                                                                        gcCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTc                                                                         <  1:456225/251‑1 (MQ=255)                                                                                                                                                                           aTCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTcgcg                                                                      >  2:259909/1‑251 (MQ=255)                                                                                                                                                                                 tctGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAAt                                                                >  1:441973/1‑251 (MQ=255)                                                                                                                                                                                        cgTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCaa                                                         <  1:472724/251‑1 (MQ=255)                                                                                                                                                                                         gTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAAc                                                        <  2:169798/251‑1 (MQ=255)                                                                                                                                                                                               ccaaccTCAATCAAACCCATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑GTTCTACTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATTCCGAAGTCGATGAACAGCGGATG‑AATCTCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTTAAATCATGCAACGTTATc                                                  <  1:112059/246‑1 (MQ=255)                                                                                                                                                                                                  aaCTCAATCAAACACATACCAGTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGa                                               <  2:110193/251‑1 (MQ=255)                                                                                                                                                                                                   aCTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCATAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAACACGCGTTGATGGCTTGCTGTTACTGTGCACCGCAACACATCAACCTTCGCGTGAAATCATGCAACGTTATCCGAc                                              >  2:192051/1‑251 (MQ=255)                                                                                                                                                                                                   aCTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGAc                                              >  1:95650/1‑251 (MQ=255)                                                                                                                                                                                                      caatcaaACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGAAAGt                                           >  1:399840/1‑251 (MQ=255)                                                                                                                                                                                                       aatcaaACACATACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTg                                          >  1:204170/1‑251 (MQ=255)                                                                                                                                                                                                              cacaTACCATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGATTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGAGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAAAGTTATCCGACAGTGCCTACtg                                   >  1:246156/1‑251 (MQ=255)                                                                                                                                                                                                                     cATTGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTgatgat                            <  1:153219/251‑1 (MQ=255)                                                                                                                                                                                                                        tGGCATGTTGATCACTGCCAGTACCAATC‑TTTCTAGTCAGAACTGGTGCGTGGCGTTGAACGCAGCTTCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCGTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATgg                          <  1:383172/250‑1 (MQ=255)                                                                                                                                                                                                                            aTGTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTgg                     <  1:296473/251‑1 (MQ=255)                                                                                                                                                                                                                              gTTGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTggg                    <  1:370264/250‑1 (MQ=255)                                                                                                                                                                                                                               ttGATCACTGCCAGTACCAATC‑TTTCTATTCAGCACTGGTGCGTGGCGTTGCACGCAGCTGCTTCGAAAGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGACCAGCGGATG‑AATCGCAATCTGGAAACGCTGATg                                                                                                                              >  2:35345/1‑143 (MQ=255)                                                                                                                                                                                                                               ttGATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATg                                                                                                                              <  1:35345/143‑1 (MQ=255)                                                                                                                                                                                                                                 gagcACTGCCAGTACCAATC‑TTTCTAGTCAGAACTGGTTCGTGGCGTTGAACGCAGCTGCTTCGATCGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACTGCGGATG‑AATCTCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCTCCTAAACGCATCAACCTTCGCGTGAAATCATGCATCGTTATCCGACAGTGCCTACTGTGATGATGGCCTGGGCTcc                <  2:343342/248‑1 (MQ=255)                                                                                                                                                                                                                                 gATCACTGCCAGTACCAATC‑TTTCTATTCAGAACTGGGGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGACTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AAGCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGTCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTcc                <  1:369526/251‑1 (MQ=255)                                                                                                                                                                                                                                    cACTGCCAGTACCAATC‑TTTCGATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCGGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTTCGtt             <  2:441973/251‑1 (MQ=255)                                                                                                                                                                                                                                       tGCCAGTACCAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACTCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGACACGCATCAACCTTCGCGTGAATTCATGCAACGTTATCCGACAGTGCCTACTGTGGTGATGGCCTGGGCTCCGTTCg           >  2:161045/1‑250 (MQ=255)                                                                                                                                                                                                                                               ccAATC‑TTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTGATAGGCTCGTCCTTTGCAAGACCGAAGGCGATGAACAGCGGATG‑AAGCGCAGTCTGGAAACGCTGATGCAAAGACTCGGTGATGGCTTGCGGTTACTGTTCACCGAAGCGCATCATCCTTCGCGTGAAATCATGCAACGTGAGCGGATGGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCGACa  <  1:336561/251‑1 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       GACGCATTTTTAGACAGGCAGAGGTGACGCTTGGCTACAATGAAAGATGTTGCCCGCCTGGCGGGCGTTTCTACCTCAACAGTTTCTCACGTTATCAATAAAGATCGCTTCGTCAGTGAAGCGATTACCGCCAAAGTTGAAGCGGCGATTAAAGAACTCAATTACGCGCCATCAGCTCTGGCGCGTAGCCTCAAACTCAATCAAACACATACCATTGGCATGTTGATCACTGCCAGTACCAATCCTTTCTATTCAGAACTGGTGCGTGGCGTTGAACGCAGCTGCTTCGAACGCGGTTATAGTCTCGTCCTTTGCAATACCGAAGGCGATGAACAGCGGATG‑AATCGCAATCTGGAAACGCTGATGCAAAAACGCGTTGATGGCTTGCTGTTACTGTGCACCGAAACGCATCAACCTTCGCGTGAAATCATGCAACGTTATCCGACAGTGCCTACTGTGATGATGGACTGGGCTCCGTTCGATGGCGACA  >  NC_000913/3938197‑3938686

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: ‑