Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	731,569	C→T	R663W (CGG→TGG)	rhsC →	rhs element protein RhsC

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	731,569	0	C	T	98.3%	176.6 / NA	59	R663W (CGG→TGG)	rhsC	rhs element protein RhsC
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base T (27/31);  total (28/31)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

GTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACC‑GCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCATATCAGCGAAGGGCACCGGGTGACGGTGCACTATG  >  NC_000913/731332‑731809                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  gTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAgg                                                                                                                                                                                                                                       >  1:269671/1‑250 (MQ=255)  tGAAAGACACGCAGGGCCATGAAACGCGGTATGCATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGa                                                                                                                                                                                                                                      >  1:160960/1‑250 (MQ=255)         cacGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCt                                                                                                                                                                                                                              <  1:181573/251‑1 (MQ=255)           cGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACCACACCTTCCGTTACGATGACCTCGACTGGCTGATACCGGaacacggctct                                                                                                                                                                                                                            >  2:272879/1‑242 (MQ=14)           cGCAGGGCCATGAAACGCGGAATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTGCCGTTACGATGTACTCGACTGGCTGATACAGGAATCCGGCttt                                                                                                                                                                                                                            <  2:207240/251‑1 (MQ=255)                ggCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGGGGGGAAAAGCCATCTGACCCACGCAGGGCGGCTTGAGGCCCGGTATGTAATAGGATGCGGCCGGACGGGTCATCCGCCTGCCCAGAGAAAACGGCAGCAACCCCACCTTCCGTTACGACGTACTCGACTGGCTGATACAGGAAACCGGCTTTGCCgg                                                                                                                                                                                                                       >  1:130582/1‑251 (MQ=14)                      gAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTGCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGca                                                                                                                                                                                                                  <  2:179993/250‑1 (MQ=255)                       aaaCGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGcaca                                                                                                                                                                                                                <  1:364838/251‑1 (MQ=255)                       aaaCGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGcac                                                                                                                                                                                                                 >  1:119162/1‑250 (MQ=255)                        aaCGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGcacac                                                                                                                                                                                                               >  1:112285/1‑251 (MQ=255)                           gcgGTAGGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCGGATACAGGAAACCGGCTTTGACGGCCGCACACAGc                                                                                                                                                                                                            <  2:178680/251‑1 (MQ=255)                                     tACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATcacc                                                                                                                                                                                                   >  1:153354/1‑250 (MQ=255)                                       cAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGCAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGGACCACGCAGGGCGGTCTGACGCGCAGGATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTAGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACg                                                                                                                                                                                                >  1:18247/1‑251 (MQ=255)                                           cccgacgGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCTCACACAGCGTTATCACCACGACCt                                                                                                                                                                                            <  2:293348/246‑1 (MQ=255)                                                  gTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGc                                                                                                                                                                                     >  1:56837/1‑251 (MQ=255)                                                    gacctgccgACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGTTCTTACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACGCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCa                                                                                                                                                                                    <  1:14732/250‑1 (MQ=255)                                                           ccaccGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACAGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATc                                                                                                                                                                            >  2:92579/1‑251 (MQ=255)                                                                gTCAGTGCCCCGGACGGCAGCAGAAACGGGACGCAGGACGATGCGTGGGGAAAAGCCATCTGTTCCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAg                                                                                                                                                                       <  2:78179/251‑1 (MQ=255)                                                                      gCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGa                                                                                                                                                                 <  1:54928/251‑1 (MQ=38)                                                                          cGGACGGCAGCAGAAACGGGACCCAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAg                                                                                                                                                             <  2:255801/251‑1 (MQ=255)                                                                                    cagAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCa                                                                                                                                                   <  2:119162/251‑1 (MQ=255)                                                                                     agAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCAc                                                                                                                                                  <  2:62446/251‑1 (MQ=255)                                                                                           gggACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTg                                                                                                                                            >  1:221274/1‑251 (MQ=255)                                                                                                 caGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTggaaat                                                                                                                                       >  1:15557/1‑246 (MQ=255)                                                                                                    tACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACCCCACCTTCCGTTCCGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCCCACAGCGTTATCACCACGCCCTGACCGGCAAACTTATCCGCAGCGAGGCTGCGGGGCTGGTCACCCACTGGCACtctgc                                                                                                                                   >  2:139100/1‑247 (MQ=14)                                                                                                           cGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGTCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCa                                                                                                                            <  1:56391/251‑1 (MQ=255)                                                                                                                  aaaaGCCATCTGTACCCCGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCCCCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGcc                                                                                                                     <  2:197749/251‑1 (MQ=255)                                                                                                                  aaaaGCCAGCGGGACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGcc                                                                                                                     <  2:36723/251‑1 (MQ=255)                                                                                                                             gTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCTCACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCAcgcac                                                                                                             <  1:344884/248‑1 (MQ=25)                                                                                                                             gTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACATTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGCCCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGAACTATGACGAAGCAGACCGCCTCAcgcacc‑gc                                                                                                         >  1:299215/1‑251 (MQ=255)                                                                                                                                   cGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCCCGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCCCTGGCACTATGCCGAAGCAGACCGCCTCACGCACC‑GCACCGTg                                                                                                   >  2:136538/1‑251 (MQ=255)                                                                                                                                   cGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACC‑GCACCGTg                                                                                                   <  1:355691/251‑1 (MQ=21)                                                                                                                                     cAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCCCGCACC‑GCACCGTGa                                                                                                  >  2:186863/1‑250 (MQ=25)                                                                                                                                        ggCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGCTGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACC‑GCACCGTGAATgg                                                                                              >  2:278379/1‑251 (MQ=25)                                                                                                                                          cGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACC‑GCACCGTGAATGGCg                                                                                            <  1:87509/251‑1 (MQ=37)                                                                                                                                             tCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACC‑GCACCGTGAATGGCGaa                                                                                          >  1:341566/1‑250 (MQ=37)                                                                                                                                                gACGCGCAGTATGGAAGACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACGCCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCGTGGACGGCCGCACACAGCGTTAGCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCGGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCCCC‑GCACCGTGAATGGCGAAACCg                                                                                      <  2:46322/251‑1 (MQ=11)                                                                                                                                                     gcAGGATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACC‑GCACCGTGAATGGCGAAACCGCagag                                                                                 <  2:115810/251‑1 (MQ=37)                                                                                                                                                            agAATACGATGCTGCCGGGCGGGTCATCCGCCTGACCAGGGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACC‑GCACCGTGAATGGCGAATCCGCAGAGCGGTgg                                                                           <  2:28290/249‑1 (MQ=255)                                                                                                                                                             gAATACGATGCTGCCGGACGGGTCCTCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCCGCGAGGATGAGGGGCTGGTCACCCCCTGGCACCATGACGAAGCCGACCGCCTCAC‑CACCCGCACCGTGACTGGCGAAACCGCCGAGCGGTGGcc                                                                         >  2:347785/1‑250 (MQ=11)                                                                                                                                                                        tgcCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACC‑GCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGAAc                                                              >  1:13767/1‑251 (MQ=37)                                                                                                                                                                        tgcCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCAAACCAACTTCCGTTACGCTGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGCCCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACC‑GCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGAAc                                                              >  1:30795/1‑251 (MQ=255)                                                                                                                                                                                gggTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCCCCGTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAACCTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACC‑GCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGAACGTggctgg                                                      <  1:295457/251‑1 (MQ=255)                                                                                                                                                                                  ggcATCCGCCGGACCATTGAAAACGGCTGCCCCACCACCTTCCGTTACGATGTCCTCGACTGGCTGATACAGGAACCCGGCTTTGACGGCCGCACACAGCGTTATCTCCGCGCCCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCGCCCACTGGCACGATGACGAAGCAGACCGCCTCACGCACC‑GCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGAACGCgtctggct                                                    <  2:332638/249‑1 (MQ=14)                                                                                                                                                                                        cGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACC‑GCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGAACGCGGCTGGCTgac                                                 >  1:190330/1‑248 (MQ=39)                                                                                                                                                                                                cAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACC‑GCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGcc                                      <  2:75020/251‑1 (MQ=38)                                                                                                                                                                                                  gTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCCGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGCAGCAGACCGCCTCACGCACC‑GCACCGTGAATGGCGAAACCGCAGAGCGGGGGCAGTATGACGAACGCGGCTGGCTGCCAGACATCAGCCat                                    >  2:135806/1‑251 (MQ=255)                                                                                                                                                                                                   tGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACC‑GCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCat                                    >  1:321930/1‑250 (MQ=39)                                                                                                                                                                                                    gAAAACGGCAGCCACACCACCTTCCTTTACGATGTACTCGACTGGCTGAGACAGGAAACCGGCTTTGACGGCCGCACACGGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCTCGCACC‑GCACCGTGAATGGCGAAACCGCAGAGCTGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCatat                                  <  2:261242/251‑1 (MQ=255)                                                                                                                                                                                                     cacaCGGCAGCCACACCACCTTACGTTACGATGTACTCGACTGGCAGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCCATATTACGAAGCAGACCGCCTCACGCACC‑GCACCGGGAATGGCGAAACCGCAGAGCGGTGGCAGCATGACGACCGCGGCTGGCGGAAAGACATCAGCCaagtc                                 >  2:247577/4‑247 (MQ=14)                                                                                                                                                                                                            cAGCCACCCCACCTGCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGTCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACC‑GCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCATATCAGCGAAg                          <  2:273887/251‑1 (MQ=255)                                                                                                                                                                                                                cacaCCACCTTCCGTTACGCTGTACTCGACTGGCTGATCCAGGAACCCGGCTTTGACGGCCGCCCACCGCGGTATCACCCCGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGAGGCTGGTCACCCACTGGAACTATGACGAAGCAGACCGCCTCACGCCCC‑GCACCGCGAATGACGAAACTGCAGAGCGGTGGCAGGTTGACGAACGCGGCCGGCTGACAGACATCAGccctttcgtcgaaggtcc                      >  2:51310/1‑235 (MQ=14)                                                                                                                                                                                                                cacaCCACCTTCCGTTACGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGGCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACC‑GCAACGTGAAGTGCGACACAGCAGAGCGGAGGCAGTATGACGAACGCGACTGGCTGCCAGACATCAGcccactcaactaagggcc                      >  2:126685/1‑235 (MQ=11)                                                                                                                                                                                                                   accgcctgccGTTACGATGTACTCGACTGTCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCCCCACGACCTGACCGGCAAACTGATCCGCAGCGAGGATGAGGGGCTGGTCTCCCACTGGCACTATGACGAAGCAGACCGCCTCACGCGCC‑GCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCATATCAGCGAAGGGCACCg                   <  2:56837/243‑1 (MQ=14)                                                                                                                                                                                                                   accaccTTCCGTTACGATGTACGCGACTGGCTGATACAGGAACCCTGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTTCCCGGCAAACTTATCCGCAGCGAGGATGAGGTGCTGGTCACCCACTGGCACTAGGACGAAGCAGACCGCCTCACGCACC‑GCTCCGTGAATGGCGGAACCGCAGAGCTGTGTCAGTATGACGAACGCGGCTTGCTGACTGACATCAGCCATATCAGCGAAGGGCACCg                   <  2:23623/251‑1 (MQ=14)                                                                                                                                                                                                                                aCGATGTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACC‑GCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCATATCAGCGAAGGGCACCGGGTGACGGTGCAc      <  1:16310/251‑1 (MQ=255)                                                                                                                                                                                                                                    tGTACTCGACTGGCAAATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTCTGACGAAGCAGACCGCCTCACGCACC‑GCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCATATCAGCGAAGGGCACCGGGTGACGGTGCACTATg  >  1:88261/1‑251 (MQ=255)                                                                                                                                                                                                                                     gTACTCGACTGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACC‑GCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCATATCAGCGAAGGGCACCGGGTGACGGTGCACTATg  <  1:103029/250‑1 (MQ=255)                                                                                                                                                                                                                                     gTACTCGACTGGCTGATACAGGAAACCGGCGTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACC‑GCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCATATCAGCGAAGGGCACCGGGTGACGGTGCACTATg  <  1:139100/250‑1 (MQ=255)                                                                                                                                                                                                                                              |                                                                                                                                                                                                                                                  GTGAAAGACACGCAGGGCCATGAAACGCGGTATGAATACAACGCCGCCGGTGACCTGACCACCGTCATTGCCCCGGACGGCAGCAGAAACGGGACACAGTACGATGCGTGGGGAAAAGCCATCTGTACCACGCAGGGCGGTCTGACGCGCAGTATGGAATACGATGCTGCCGGACGGGTCATCCGCCTGACCAGTGAAAACGGCAGCCACACCACCTTCCGTTACGATGTACTCGACCGGCTGATACAGGAAACCGGCTTTGACGGCCGCACACAGCGTTATCACCACGACCTGACCGGCAAACTTATCCGCAGCGAGGATGAGGGGCTGGTCACCCACTGGCACTATGACGAAGCAGACCGCCTCACGCACC‑GCACCGTGAATGGCGAAACCGCAGAGCGGTGGCAGTATGACGAACGCGGCTGGCTGACAGACATCAGCCATATCAGCGAAGGGCACCGGGTGACGGTGCACTATG  >  NC_000913/731332‑731809

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