Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,239,627	G→A	intergenic (+82/‑317)	alx → / → sstT	putative membrane‑bound redox modulator Alx/serine/threonine:Na(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,239,627	0	G	A	97.8%	144.4 / ‑6.5	46	intergenic (+82/‑317)	alx/sstT	putative membrane‑bound redox modulator Alx/serine/threonine:Na(+) symporter
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (22/23);  minor base T (1/0);  total (23/23)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

CTGGCGGTGATTCTGGTGTTTATCGGTATCAAGATGCTGATTGTCGACTTCTACCATATTCCAATCGCCGTCTCGCTGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGTTTCAAAATATGACCTCTCTTTAAAATCCAGCATTTTTCGCTTCCCGAAGCTGTAACTTTCCTTATACTCGACCTTGCAAACACT  >  NC_000913/3239384‑3239871                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     cTGGCGGTGATTCTGGTGTTTATCGGTATCAAGATGCTGATTGTCGACTTCTACCATATTCCAATCGCCGTCTCGCTGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGGGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCCCCGTTTTCTCCCTCCCCCTTCCAGGGAGAGGGTCGGGATGCGGGt                                                                                                                                                                                                                                               >  2:161737/1‑251 (MQ=255)             cTGGTGTTTATCGGTATCAAGATGCTGATTGTCGACTTCTACCATATTCCAATCGCCGTCTCGCTGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCATGGAGAGGGTCGGGATGAGGGTAATTTTTCGCAc                                                                                                                                                                                                                                   >  2:347470/1‑251 (MQ=255)                            aTCAAGATGCTGATTGTCGACTTCTACCATATTCCAATCGCCGTCTCGCTGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGtt                                                                                                                                                                                                                    <  2:191868/251‑1 (MQ=255)                               aaGATGCTGATTGTCGACTTCTACCATATTCCAATCGCCGTCTCGCTGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTccc                                                                                                                                                                                                                 >  2:631013/1‑251 (MQ=255)                                 gATGCTGATTGTCGACTTCTACCATATTCCAATCGCCGTCTAGCTGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCACCGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCNCCt                                                                                                                                                                                                               >  2:593687/1‑251 (MQ=255)                                       gATTGTCGACTTCTACCATATTCCAATCGCCGTCTCGCTGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCAccc                                                                                                                                                                                                          <  1:104806/250‑1 (MQ=255)                                                     aCCATATTCCAATCGCCGTCTCGCTGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCaa                                                                                                                                                                                           <  1:13442/251‑1 (MQ=255)                                                       cATATTCCAATCGCCGTCTAGCTGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCaaa                                                                                                                                                                                          >  1:570399/1‑250 (MQ=255)                                                        atatTCCAATCGCCGTCTCGCTGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGTTGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACg                                                                                                                                                                                        >  2:132246/1‑251 (MQ=255)                                                           ttCCAATCGCCGTCTCGCTGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTTCCCCCAAACgggg                                                                                                                                                                                     >  1:455022/1‑251 (MQ=255)                                                                       tctcGCTGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTg                                                                                                                                                                          >  1:500443/1‑250 (MQ=255)                                                                              gggCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGTGTGAATTATCTGCATGATAAGCAGCGTGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACTTGGCGATGGGACTGACCGAg                                                                                                                                                                    <  2:228577/249‑1 (MQ=255)                                                                                      tgtTTGGCATTCGGGTGATGACGTTTATTATCAACGCCTGGGGGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGGAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACACTCGCCCCAAACAGGGCGGGGGGGCTGACCGATGCCTTTTTTg                                                                                                                                                          >  2:1052/1‑251 (MQ=255)                                                                                       gtTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGa                                                                                                                                                         <  2:312521/251‑1 (MQ=255)                                                                                         ttGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGAtt                                                                                                                                                       >  2:409863/1‑250 (MQ=255)                                                                                                tCTGGTGATGACGTTTATTATCAACGCCTGGGCGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGtcat                                                                                                                                                >  1:550792/1‑251 (MQ=255)                                                                                                      gatgaCGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCt                                                                                                                                          <  1:462820/251‑1 (MQ=255)                                                                                                          aCGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGa                                                                                                                                       >  1:470364/1‑250 (MQ=255)                                                                                                                     cAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTg                                                                                                                           <  1:494220/251‑1 (MQ=255)                                                                                                                         gCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCtt                                                                                                                       <  2:401186/251‑1 (MQ=255)                                                                                                                              ggTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTAttt                                                                                                                  <  2:339986/251‑1 (MQ=255)                                                                                                                                 gAATTATCGGCATGATAAGCGTCTGGGTGGATAATTTTTAATCTGCCGAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGTGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCCCCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGATTCCTTTTTTGATGTTGTTATCAGTCTGGAAGCCGCGCGTTGGCTTTAttttt                                                                                                                <  2:626297/250‑1 (MQ=255)                                                                                                                                    ttATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGt                                                                                                            >  1:137084/1‑251 (MQ=255)                                                                                                                                                     agcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAcc                                                                                                                                                                                                     <  2:555836/145‑1 (MQ=255)                                                                                                                                                     agcGGGGAGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAAcc                                                                                                                                                                                                     >  1:555836/1‑145 (MQ=255)                                                                                                                                                      gcGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCAtt                                                                                          <  2:343004/251‑1 (MQ=255)                                                                                                                                                                 aaTTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGTTTCaaa                                                                                >  2:238658/1‑250 (MQ=255)                                                                                                                                                                    ttttAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGTTTCAAAata                                                                             <  1:631013/250‑1 (MQ=255)                                                                                                                                                                         aTCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTcccc                                                                                                                                                                                                                >  1:106514/1‑114 (MQ=255)                                                                                                                                                                         aTCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTcccc                                                                                                                                                                                                                <  2:106514/114‑1 (MQ=255)                                                                                                                                                                         aTCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGTTTCAAAATATGACct                                                                       >  1:128715/1‑251 (MQ=255)                                                                                                                                                                                 aaGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTATCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGTTTCAAAATATGACCTCTCTTTa                                                                >  1:536812/1‑250 (MQ=255)                                                                                                                                                                                 aaGCCGTGTACCCTGTCATTAACATGAGCACAGTTTTCTCCCTCTCCCGTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGTTTCAAAATATGACCTCTCTTTa                                                                <  2:552705/250‑1 (MQ=255)                                                                                                                                                                                                            gCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGCTTCAAAATATGACCTCTCTTTAAAATCCAGCATTTTTCGCTTCCCGAAGc                                    >  2:571578/1‑251 (MQ=255)                                                                                                                                                                                                              aCCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCATCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGATGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGTTTCAAAATATGACCTCTCTTTAAAATCCAGCATTTTTCGCTTCCCGAAGCTg                                  >  2:370697/1‑251 (MQ=255)                                                                                                                                                                                                                  ttttCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGTGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCAAGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTCAGTTTCAACATATGACCTCTCTTTAAAATCCAGCATTTTTCGCTTCCCGAAGCTGTAAc                              >  1:569745/1‑251 (MQ=255)                                                                                                                                                                                                                   tttCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGTTTCAAAATATGACCTCTCTTTAAAATCCAGCATTTTTCGCTTCCCGAAGCTGTAAc                              <  2:341587/250‑1 (MQ=255)                                                                                                                                                                                                                    ttCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGTTTCAAAATATGACCTCTCTTTAAAATCCAGCATTTTTCGCTTCCCGAAGCTGTAACt                             <  2:458653/250‑1 (MQ=255)                                                                                                                                                                                                                              ccctTCAAGGGAGAGGGTCGGGATGAGGGAAATTTTTCGCACCGATGCTGGACTGTTCCCCTCAACCTAACCCTCTCCCCAAACGGGGCGGGGGGACTGCCCGAGTGCTTTTTTGATGTTGTCATCAGTCTGGAAGACGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAAACATAAAGTTTCAAAATATGAGCTCTCTTTAAAATCCAGCATTTTTCGATTCCCGAACCTGTAAATTTCCTTATACt                  >  2:494814/1‑251 (MQ=255)                                                                                                                                                                                                                                 ttCCAGGGAGAGGGACGTGTTGAGTGTAATTTTTCGCACAGGTGATGGCCTGGTCCCATCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGTATGCCCGCGGCCTTTTTTTATTTTTTAATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGGTTCAAAATATGACCTCTCTTTAAACTCCAGCATTTTTCGCTTCCCGAATCTGTAACTTTCCTTATACTca                >  2:654369/1‑249 (MQ=255)                                                                                                                                                                                                                                    cAGGGAGAGGGTCGGGATGAGGGTAATTTGTCGCACCGATGCTGGCCTGTTCCCCGCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGTTTCAAAATATGACCTCTCTTTAAAATCCAGCAGTTTTCGCTTCCCGAAGCTGTACCTTTCCTTATGCTCGatct            <  2:116474/251‑4 (MQ=255)                                                                                                                                                                                                                                     aGGGAGAGGGTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGTTTCAAAATATGACCTCTCTTTAAAATCCAGCATTTTTCGCTTCCCGAAGCTGTAACTTTCCTTATACTCGACCtt           <  1:286157/251‑1 (MQ=255)                                                                                                                                                                                                                                            gggTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGTTTCAAAATATGACCTCTCTTTAAAATCCAGCATTTTTCGCTTCCCGAAGCTGTAAATTTCCTTATACTCGACCTTGCAAaca    <  2:475900/251‑1 (MQ=255)                                                                                                                                                                                                                                             ggTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGTTTCAAAATATGACCTCTCTTTAAAATCCAGCATTTTTCGCTTCCCGAAGCTGTAACTTTCCTTATACTCGACCTTGCAAacac   <  1:394905/251‑1 (MQ=255)                                                                                                                                                                                                                                             ggTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCACCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGTTTCAAAATATGACCTCTCTTTAAAATCCAGCATTTTTCGCTTCCCGAAGCTGTAACTTTCCTTATACTCGACCTTGCAAacac   <  1:161737/251‑1 (MQ=255)                                                                                                                                                                                                                                              gTCGGGATGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGTTTCAAAATATGACCTCTCTTTAAAATCCAGCATTTTTCGCTTCCCGAAGCTGTAACTTTCCTTATACTCGACCTTGCAAACACt  <  1:428306/251‑1 (MQ=255)                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                     CTGGCGGTGATTCTGGTGTTTATCGGTATCAAGATGCTGATTGTCGACTTCTACCATATTCCAATCGCCGTCTCGCTGGGCGTGGTGTTTGGCATTCTGGTGATGACGTTTATTATCAACGCCTGGGTGAATTATCGGCATGATAAGCAGCGGGGTGGATAATTTTTAATCTGCCTAAGCCGTGTACCCTGTCATTAACATGAGCACCGTTTTCTCCCTCTCCCTTCCAGGGAGAGGGTCGGGGTGAGGGTAATTTTTCGCACCGATGCTGGCCTGTTCCCCTCACCCTAACCCTCTCCCCAAACGGGGCGAGGGGACTGACCGAGTCCTTTTTTGATGTTGTCATCAGTCTGGAAGCCGCACGTTGGCTTTATTTTTATGTCAAAGAAATGTAACCATTAAGTTTCAAAATATGACCTCTCTTTAAAATCCAGCATTTTTCGCTTCCCGAAGCTGTAACTTTCCTTATACTCGACCTTGCAAACACT  >  NC_000913/3239384‑3239871

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