Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,067,555	C→T	E574E (GAG→GAA)	yihQ ←	sulfoquinovosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,067,555	0	C	T	96.2%	176.7 / NA	54	E574E (GAG→GAA)	yihQ	sulfoquinovosidase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (23/28);  minor base A (2/0);  total (26/28)
Fisher's exact test for biased strand distribution p-value = 2.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCGATGGGCGCATTAACGGTAACTTCCCCGCCCCGGAACGCTTCACCCGTCCAGGCGTGGACCCAGTTATCCTCCGGCAGATAGAGCGTCCAGTCGCTACGGCCTTCTTCATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCCTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGG  >  NC_000913/4067308‑4067790                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            tGCCGATGGGCGCATTAACGGTAACTTCCCCGCCCCGGACCGCTTCACCCGTCCAGGCGTGGACCCAGTTATCCTCCGGCAGATAGAGCGTCCAGTCGCTACGGCCTTCTTCATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCttct                                                                                                                                                                                                                                          >  2:224324/1‑251 (MQ=255)   ccGATGGGCGCATTAACGGTAACTTCCCCGCCCCGGAACGCTTCACCCGTCCAGGCGTGGACCCAGTTATCCTCCGGCAGATAGAGCGTCCAGTCGCTACGGCCTTCTTCATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCttt                                                                                                                                                                                                                                        <  2:116184/251‑1 (MQ=255)      aTGGGCGCATTAACGGTAACTTCCCCGCCCCGGAACGCTTCACCCGTCCAGGCGTGGACCCAGTTATCCTCCGGCAGATAGAGCGTCCAGTCGCTACGGCCTTCTTCATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTc                                                                                                                                                                                                                                       >  1:180789/1‑249 (MQ=255)        gggCGCATTAACGGTAACTTCCCCGCCCCGGAACGCTTCACCCGTCCAGGCGTGGACCCAGTTATCCTCCGGCAGATAGAGCGTCCAGTCGCTACGGCCTTCTTCATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTATTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCATCGCGCCGGCTTCTTTCAggg                                                                                                                                                                                                                                   >  2:140139/1‑250 (MQ=255)         ggCGCATTAACGGTAACTTCCCCGCCCCGGAACGCTTCACCCGTCCAGGCGTGGCCCCAGTTATCCTCCGGCAGATAGAGCGTCCAGTCGCTACGGCCTTCTTCATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGGGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGt                                                                                                                                                                                                                                   >  2:115657/1‑250 (MQ=255)                       aaCTTCCCCGCCCCGGCACGCTTCACCCGTCCAGGCGTGGACCCAGTTATCCTCCGGCAGATAGAGCGTCCAGTCGCTACGGCCTTCTTCATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGAGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCGTCTTTCAGGTAAGGTTTCAgggggg                                                                                                                                                                                                                    >  1:200235/1‑248 (MQ=255)                                    cgGAACGCTTCACCCGTCCAGGCGTGGACCCAGTTATCCTCCGGCAGATAGAGCGTCCAGTAGCTACGGCCTTCTTCATGCACCGGGGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGGGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCATGGTGGTGACGACggtggt                                                                                                                                                                                                       >  2:204883/1‑251 (MQ=255)                                       aaCGCTTCACCCTGCCAGGCGTTGACCCATTTATCCTCCGGCAGATAGAGCGTCCAGTCGCCACGACCTTCTTCATGCACCTGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGTGTGTAAGTGTGCGCATCGTCTGCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGTACTTCGCATTCAGCGCGACGGCTTCTTTCATGTAAGGTTTCAGGGTGGTGAAGATGGTTGTCAt                                                                                                                                                                                                    <  2:85140/251‑1 (MQ=255)                                       aaCGCTTCACAAGTCCAGGCGTTGACCCATTTATCCTCCGGCAGCTAGAGCGTCCAGTCGCTACGGCCTTCTTCATGCACCGGAGCGACCAGCATGTCGCGACCTAACAGGTACTGATATTTCAGGCTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAAACGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGAGGAGAAGCCGGTGGTCAt                                                                                                                                                                                                    >  2:230653/1‑251 (MQ=255)                                            ttCACGCGTCCAGGTGTGCACCCCGTTATCCTCAGGCAGATTGAGCATCCCGTCGCTCCGGCCTTCTACATGCCCCGGAGCGACCCGATTGTCGCCACCTAACCGGTACTGATACTTCCGGGACTGAGTGTGACTATCGTCTCCGTAATGCAGGGACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATAGAGCGCGACGGCATATTTCAGATAAGGCTTCAGGGCGCTGAAGACGGGCGTCATaacgg                                                                                                                                                                                               >  2:41877/1‑247 (MQ=255)                                               ccccGTCCAGGCGTGGACCCAGTTATCCTCCGGCAGATAGAGCGTCCAGTCGCTACGGCCTGCTTCATGCCCCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGCCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATGCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGa                                                                                                                                                                                            <  2:249176/250‑1 (MQ=255)                                                cccGTCCAGGCGTGGACCCAGTTATCCTCCGGCAGATAGAGCGTCCAGTCGCTACGGCCTTCTTCATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGaa                                                                                                                                                                                           <  1:29704/251‑1 (MQ=255)                                                 ccGTCCAGGCGTGGACCCAGTTATCCTCCGGCAGATAGAGCGTCCAGTCGCTACGGCCTTCTTCATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGaaa                                                                                                                                                                                          <  1:247295/251‑1 (MQ=255)                                                                    gTTATCCTCCGGCAGATAGAGCGTCCAGTCGCTACGGCCTTCTTCATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGt                                                                                                                                                                       <  1:178112/251‑1 (MQ=255)                                                                    gTTATCCTCCGGCAGATAGAGCGTCCAGTCGCTACGGCCTTCTTCATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGt                                                                                                                                                                       >  1:38280/1‑251 (MQ=255)                                                                      tATCCTCCGGCAGATAGAGCGTCCAGTCGCTACGGCCTTCTTCATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTc                                                                                                                                                                      <  2:89363/250‑1 (MQ=255)                                                                          ctccGGCAGATAGAGCGTCCAGTCGCTACGGCCTTCTTCATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGcagatcggaaagagcacacgt                                                                                                                                                                 >  1:142079/1‑231 (MQ=255)                                                                          ctccGGCAGATAGAGCGTCCAGTCGCTACGGCCTTCTTCATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGc                                                                                                                                                                                     <  2:142079/231‑1 (MQ=255)                                                                            ccGGCAGATAGAGCGTCCAGTCGCTACGGCCTTCTTCATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTc                                                                                                                                                                <  2:146421/250‑1 (MQ=255)                                                                               gCAGATAGAGCGTCCAGTCGCTACGGCCTTCTTCATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCAGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAAc                                                                                                                                                            >  1:455/1‑251 (MQ=255)                                                                                cAGATAGAGCGTCCAGTCGCTACGGCCTTCTTCATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACt                                                                                                                                                           <  1:67330/251‑1 (MQ=255)                                                                                     agagCGGCCAGTCGCTACGGCCTTCTTCATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCa                                                                                                                                                       <  1:120692/250‑1 (MQ=255)                                                                                      gagCGTCCAGTCGCTACGGCCTTCTTCATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGt                                                                                                                                                     <  1:75893/251‑1 (MQ=255)                                                                                             cAGTCGCTACGGCCTTCTTCATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGcc                                                                                                                                              <  2:78654/251‑1 (MQ=255)                                                                                                  gCTACGGCCTTCTTCATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGAc                                                                                                                                         >  2:128960/1‑251 (MQ=255)                                                                                                  gCTACGGCCTTCTTCATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGAc                                                                                                                                         >  1:165551/1‑251 (MQ=255)                                                                                                       ggCCGTCTTCATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTGTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTa                                                                                                                                    <  2:259/251‑1 (MQ=255)                                                                                                         ccTTCTTCATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACGTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTGGTCGCCAGGACGGTTAcc                                                                                                                                  <  1:230653/251‑1 (MQ=255)                                                                                                           ttcttcATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGCAATGGGCGATGGTTTCTGCGTCGCCGTCAACCTGCCAGTTGTCGCCAGGCCGGTTACCtt                                                                                                                                >  1:263421/1‑251 (MQ=255)                                                                                                                cATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTgg                                                                                                                           <  1:47872/251‑1 (MQ=255)                                                                                                                          aGCGACCAGAATGTCGCGACCTAACAGGAACTGATCTTTCAGGGTGTAAGAGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTAGCCGTCAAACTGCCAGTTGTCGCAAGGACGGTTACCTTCGGGGGTGCGcaac                                                                                                                  >  2:244034/1‑248 (MQ=255)                                                                                                                           gcgccgaGAATGTCGCGACCTAACAGGTACGGAGATTTCTGGTTGGAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAGCCGGCAGTCCGGGCGTCGCAGTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCAGACGGGCGAAATGGGCGATGGGTTCGGCTTCGCGGGCAAACTGCCAGTTGTCGCCAGGATGGTCAGCTTTGTGGGTGGGcatcat                                                                                                                <  2:157915/245‑1 (MQ=255)                                                                                                                               ccAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGc                                                                                                            <  1:157349/251‑1 (MQ=255)                                                                                                                                  gAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTg                                                                                                         >  2:109919/1‑251 (MQ=255)                                                                                                                                   aaTGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGGGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGGTGTCGCCAGGCCGGTTACCTTCGTGGGTGCGCATCATCGGCGTg                                                                                                         >  2:268/1‑250 (MQ=255)                                                                                                                                            aCCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTg                                                                                                >  1:106966/1‑250 (MQ=255)                                                                                                                                            aCCTAACAGGAACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTAGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGCTGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGATTACCTTCGTGGTGGCGCATCAGCGGCGTGAAGGCGCTg                                                                                                >  2:230875/1‑250 (MQ=255)                                                                                                                                                aaCAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATc                                                                                           <  1:115657/251‑1 (MQ=255)                                                                                                                                                      tACTGATATTTAAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCCGGTAGGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGATGGCGCTGAACTCggccc                                                                                      >  2:169158/1‑246 (MQ=255)                                                                                                                                                                       gtAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGGGCGCATCATCGGCGTGAAGGCGCTGCAATCGCACCAGCGCAGCAGCAGCTCtt                                                                    >  2:203436/1‑251 (MQ=255)                                                                                                                                                                       gtAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGGGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCtt                                                                    <  1:169158/251‑1 (MQ=255)                                                                                                                                                                         aaGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAAACGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCCGCGCAGCAGCAGCTCTTtg                                                                  >  1:12630/1‑251 (MQ=255)                                                                                                                                                                               gcgcATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTgcg                                                             <  2:46118/250‑1 (MQ=255)                                                                                                                                                                               gcgcATCGTCGTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTgcgc                                                            <  2:200235/251‑1 (MQ=255)                                                                                                                                                                                cgcATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCt                                                           <  1:63727/251‑1 (MQ=255)                                                                                                                                                                                            gTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACa                                               <  1:224324/251‑1 (MQ=255)                                                                                                                                                                                                tGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGggt                                           >  2:252623/1‑251 (MQ=255)                                                                                                                                                                                                 gCAGGAACAGCGGGCGCATAATCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCGTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGTCGATGTTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGTACTGTTACCTTCGTGGGTGCGCATCAGCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGggtg                                          <  1:121661/251‑1 (MQ=255)                                                                                                                                                                                                       aCAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAACCAGGGGGGTGTa                                     >  2:22217/1‑250 (MQ=255)                                                                                                                                                                                                       aCAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTaa                                    <  1:95508/251‑1 (MQ=255)                                                                                                                                                                                                            gggCGCAAAACCGGCAGGCCGGACTTCGCATCCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGcc                               <  2:38280/251‑1 (MQ=255)                                                                                                                                                                                                             ggCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAgcg                                                                                  <  1:130945/199‑1 (MQ=255)                                                                                                                                                                                                             ggCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGCAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAgcg                                                                                  >  2:130945/1‑199 (MQ=255)                                                                                                                                                                                                                        ggCAGGCTGGACTTCGCATTCAGCGCGACGGCTTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCtgtg                   >  2:160643/1‑251 (MQ=255)                                                                                                                                                                                                                                          ttCAGAGCGACGGCATCCTTAAGGTAAGGCTTCAGGGTGGTGAAGTCGGTGGTCATACGGGCGAAATGGGCGATGGTTTCAGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGCGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCTGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCgg  >  2:248804/1‑250 (MQ=255)                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                            TGCCGATGGGCGCATTAACGGTAACTTCCCCGCCCCGGAACGCTTCACCCGTCCAGGCGTGGACCCAGTTATCCTCCGGCAGATAGAGCGTCCAGTCGCTACGGCCTTCTTCATGCACCGGAGCGACCAGAATGTCGCGACCTAACAGGTACTGATATTTCAGGGTGTAAGTGTGCGCATCGTCTTCGTAATGCAGGAACAGCGGGCGCATAACCGGCAGGCCGGACTTCGCATTCAGCGCGACGGCCTCTTTCAGGTAAGGTTTCAGGGTGGTGAAGACGGTGGTCATACGGGCGAAATGGGCGATGGTTTCTGCGTCGCCGTCAAACTGCCAGTTGTCGCCAGGACGGTTACCTTCGTGGGTGCGCATCATCGGCGTGAAGGCGCTGAAATCGCACCAGCGCAGCAGCAGCTCTTTGCTGCGCTTCATCTCAAACAGGGTGGTGTAACCGCCAATGTCGCTGTGGTGCAGGCCATGTCCGG  >  NC_000913/4067308‑4067790

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