Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,621,090	C→T	R633R (CGC→CGT)	rhsB →	rhs element protein RhsB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,621,090	0	C	T	97.8%	81.6 / NA	45	R633R (CGC→CGT)	rhsB	rhs element protein RhsB
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (25/19);  total (26/19)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTCAGTTGCTGAGCTTCACCGACTGTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCGGTGGCAGTATGACGAACGTGGCTGGCTGACAG  >  NC_000913/3620801‑3621373                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                            gTCAGTTGCTGAGCTTCACCGACTGTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAAt                                                                                                                                                                                                                                                                                   >  1:513761/1‑300 (MQ=14)      ttGCTGAGCTTCACCGACTGTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGNTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACTGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGAt                                                                                                                                                                                                                                                                              <  2:652595/300‑1 (MQ=14)           gAGCTTCACTGACTGTTCCGGTTATGTAACCCGTTATGATCATGACCGCTTCGGGCAGATGACTGCGATGCACTACGATGAAGTGCTGAGTCAGTACCGCGCATACGACAGACGTGGACAGTGAATTGCCGGGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACTGGACATAGTACGATGCGTGGGGAAAGGCCGTCCGTAACACGCAGGGCGGGCTGACGCGTAGTATGGAATACGAtgctgc                                                                                                                                                                                                                                                                         <  1:677294/300‑1 (MQ=11)           gAGCTTCACCGACTGTTCCGGTTATGTAACCCGTTATGGCCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGCATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCAGTCCGTACCACGAAGGGCGGGCTGATGCGTAGTATGGAATACGAtgct                                                                                                                                                                                                                                                                           >  2:212056/1‑298 (MQ=9)           gAGCTTCACCGACTGTTCCGGTTATGTAACCCGTTATGGCAATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGATATGAATACAACATCGCCGGTGACCTGACCGCCGTCCTTGCCCCGGACGGCAGCAGAAACGGGACACAGCACGATGCGTGGGGAAAGGCCGTCCGTAACACGCAGGGCGGGCTGACGCGTAGTATGGAATACGAtgctgc                                                                                                                                                                                                                                                                         >  2:212552/1‑300 (MQ=255)                         gTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCACGCATACGAAAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATcc                                                                                                                                                                                                                                                           >  1:53368/1‑300 (MQ=14)                             cGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCACCGTCATTGCACCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGACGTCCGTACCAAGCAGGGCGGGCTGACGCGTAGTATGGCATACGATGCTGCCGGACGGGTCATCCGCCt                                                                                                                                                                                                                                                       >  2:140249/1‑300 (MQ=9)                                           gtgatgaccatgaccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTNAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACg                                                                                                                                                                                                                                         <  2:671035/297‑1 (MQ=14)                                           gTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAAc                                                                                                                                                                                                                                          >  1:671035/1‑299 (MQ=14)                                                gaccatgaccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGAACAGTGAAAACGGCAGc                                                                                                                                                                                                                                    >  1:680181/1‑300 (MQ=14)                                                gaccatgaccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCAAGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAAAGGCAGc                                                                                                                                                                                                                                    >  2:280506/1‑300 (MQ=14)                                                     tgaccGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCAGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGNCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGAAGGGTAATCCGCCTGACCAGTGAAAACGACAGCCacac                                                                                                                                                                                                                               >  2:659872/1‑300 (MQ=9)                                                               gggCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGt                                                                                                                                                                                                                     >  1:682859/1‑300 (MQ=14)                                                                                        aGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATAc                                                                                                                                                                                            >  1:467287/1‑300 (MQ=14)                                                                                                   aGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCAGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGc                                                                                                                                                                                 >  2:756743/1‑300 (MQ=14)                                                                                                             gcgcATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATCCAGGAAACCGGCTTTGACGGcc                                                                                                                                                                       >  1:344370/1‑300 (MQ=14)                                                                                                                 aTACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGTTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGcac                                                                                                                                                                   <  2:170751/300‑1 (MQ=14)                                                                                                                 aTACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGcac                                                                                                                                                                   <  1:395732/300‑1 (MQ=14)                                                                                                                   aCGACAGCTGTGGACAGTTAATTGCCGTTAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGcacac                                                                                                                                                                 <  2:147485/300‑1 (MQ=14)                                                                                                                   aCGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCNGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGcacac                                                                                                                                                                 <  2:618052/300‑1 (MQ=14)                                                                                                                         gCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGtt                                                                                                                                                           >  1:798123/1‑300 (MQ=14)                                                                                                                         gCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGAACGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGtt                                                                                                                                                           >  1:249418/1‑300 (MQ=14)                                                                                                                            gTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATc                                                                                                                                                        <  1:627896/300‑1 (MQ=14)                                                                                                                                   gTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAACCGTCCGTACCCCGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGACTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATccccac                                                                                                                                                   >  1:238952/1‑298 (MQ=9)                                                                                                                                                           gggCCATGAAGCGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAAAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATAAGATACTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCAACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGATTTGACGGCCGCACCTAGCGTTACCACCACGAACTGa                                                                                                                                            >  2:112727/1‑281 (MQ=11)                                                                                                                                                           gggCAATGAAACGCGGTATGAATACAACATTGCCGGTGACCTGACCGCCGTCATTGCACCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACAGGCTTTGACGGACGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATcagaa                                                                                                                          >  2:107811/1‑295 (MQ=255)                                                                                                                                                                         ggTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGt                                                                                                            <  1:312367/299‑1 (MQ=14)                                                                                                                                                                          gTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGNAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCt                                                                                                          >  1:169841/1‑300 (MQ=14)                                                                                                                                                                             tGAATACAACATCGCCGGTGACCTGACCGCCTTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCTTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGt                                                                                                       <  2:344370/300‑1 (MQ=14)                                                                                                                                                                              gAATACAACATCGCCGGTGACCTGACTGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGGTGCGTGGGGAAAGGCCGTCCGTACCAAGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGAACAGTGAAAACGGCAGCCACACCAACTTCCGTTACGATGTACTTGACCGGTTGCTACAGGAAACCGGCTTTGACGGCCGCACATAGCGTTATCACCACGCCCTGACCGGCAAACTCATCCGCCGCGAGAATGAGGGTCTGGTc                                                                                                      >  2:246061/1‑300 (MQ=11)                                                                                                                                                                                                            gTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGACTGACCAGTGAAAACGGCGGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGAcc                                                                          <  1:792428/298‑1 (MQ=14)                                                                                                                                                                                                            gTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGACTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGAcc                                                                          <  1:792486/298‑1 (MQ=14)                                                                                                                                                                                                                       ggACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCAcgc                                                                 <  1:816251/296‑1 (MQ=14)                                                                                                                                                                                                                              agcagAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTTACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGt                                                       <  1:212552/299‑1 (MQ=14)                                                                                                                                                                                                                              agcagAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGt                                                       <  1:212056/299‑1 (MQ=14)                                                                                                                                                                                                                                      cGGGGCACAGTACGATGCGTGGGGAAAGACCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCATGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGGGGATGAGGGTCTGGTCACCCACTGGCACTATGACGACGCAGACCGCCTCACGCCCCGCACCATGAAGGGt                                                >  1:537590/1‑298 (MQ=11)                                                                                                                                                                                                                                      cGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGa                                              >  1:537430/1‑300 (MQ=14)                                                                                                                                                                                                                                         gACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAAc                                           <  1:483972/300‑1 (MQ=14)                                                                                                                                                                                                                                             caGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCa                                       >  1:559554/1‑300 (MQ=14)                                                                                                                                                                                                                                             caGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCAACACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCa                                       >  2:763090/1‑300 (MQ=14)                                                                                                                                                                                                                                                      gCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACTGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCTTGAAGGGTGAAACCGCAGAGCGGTgg                              <  2:526049/300‑1 (MQ=14)                                                                                                                                                                                                                                                              aaaGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCGGTGGCAGTATGa                      <  2:775457/300‑1 (MQ=14)                                                                                                                                                                                                                                                                    cgtccgtACCACGCAGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCGGTGGCAGTATGACGAAc                 >  1:626799/1‑299 (MQ=14)                                                                                                                                                                                                                                                                                  aGGGCGGGCTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACAACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCGGTGGCAGTATGACGAACGTGGCTGGCTgacag  <  1:454122/300‑1 (MQ=14)                                                                                                                                                                                                                                                                                   gggcgggcTGACGCGTAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGAACAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGAACGGCTGATACAGGAAACCGGCTTTGACGGCAGCACACAGCGTTATCACCACGACCTGACCGGCCAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACACACTGGCACTCTGACGAAGCAGACCGCCCCACGCACCGCACCGTGAAGGGTGAAAACGCAGAGCGGTGGAAGTATGACGCACGTggctggctcca    >  2:85586/1‑294 (MQ=11)                                                                                                                                                                                                                                                                                                  |                                                                                                                                                                                                                                                                                            GTCAGTTGCTGAGCTTCACCGACTGTTCCGGTTATGTAACCCGTTATGACCATGACCGCTTCGGGCAGATGACGGCGGTGCACCGCGAGGAAGGGCTGAGTCAGTACCGCGCATACGACAGCCGTGGACAGTTAATTGCCGTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACATCGCCGGTGACCTGACCGCCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAGGCCGTCCGTACCACGCAGGGCGGGCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTTGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGTCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACCGCACCGTGAAGGGTGAAACCGCAGAGCGGTGGCAGTATGACGAACGTGGCTGGCTGACAG  >  NC_000913/3620801‑3621373

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