Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,614,208	A→G	20.4%	E243G (GAG→GGG)	yjjV →	putative DNase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,614,208	0	A	G	20.4%	96.4 / 2.5	52	E243G (GAG→GGG)	yjjV	putative DNase
Reads supporting (aligned to +/- strand):  ref base A (20/19);  new base G (7/5);  total (28/24)
Fisher's exact test for biased strand distribution p-value = 7.49e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CGAACGGTTTGTACAGCTGGGCTACAAAATT‑GGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑CGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATC  >  NC_000913/4613978‑4614434                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                    cGAACGGTTTGTACAGCGGGGCAACAAAATT‑GGCGTAGGCGGTACTATCACCTATCCAAGCGCCAGTAAACCCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTgcg                                                                                                                                                                                                                              <  1:901176/238‑1 (MQ=255)         ttGTACAGCTGGGCTACAAAATT‑GGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGGGATTGCGCAAGCGTTg                                                                                                                                                                                                                     >  1:145257/1‑239 (MQ=255)               cgCTGGGCTAAAAAAGT‑GGCGTTGGCGTTCCAATCCCCTTTCCACGGGCCAGTAAAACCCGGGATGTCATCGGAAAATTGCCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCCGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTtaat                                                                                                                                                                                                               <  2:154499/238‑1 (MQ=255)                gCTGGGCTACAAAATT‑GGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGGGTTGCGCCGGGAACCGGCGGATGGGGTTTCGCAAGCGTTGTTtaata                                                                                                                                                                                                              >  1:877008/1‑239 (MQ=255)                gCTGGGCTACAAAATT‑GGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCCTGTGTTCGCCGGGCTTTGCGGGTTGCCCCGGGGACCGGCGGTCGAAATTGCGGAAACGGTGGTtaaaa                                                                                                                                                                                                              >  2:89728/1‑237 (MQ=255)                  ttggcgaaaaaaTT‑GGCTGTGGCCGTCATACCACCTCACCCCGCGCCAGTAAAACCCGCGATGGCATCGAAAAATTACCGCGGGCATCGTTATGGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGGTTTCCGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCGGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTtaataa                                                                                                                                                                                                             <  2:245869/230‑1 (MQ=255)                  ggggCTACAAAATT‑GGCGTGGGCTGTCCTACCCCGCACCCACGTGGCGGGAAAACCCGGGATTTCAGCGCAAAATTACCGCGGGCAGCGTTATTGGTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTGCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTtaataa                                                                                                                                                                                                             <  2:270788/237‑1 (MQ=255)                    ggCTACAAAATT‑GGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGAGGAGATTGCGCAAGGGTGGCTTAATAaca‑c                                                                                                                                                                                                         >  1:601536/1‑239 (MQ=255)                      cTAAAAAATT‑GGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCg                                                                                                                                                                                                                        <  1:185187/223‑1 (MQ=255)                          aaaaTT‑GGCTTAGGCGGTACTTTCCCCCACCCCGGCGCCATTAAACCCCGCGATGGCATCGCAAAATTCCCGCTGGCATCGTTATTGCTGGAAACCGGCGCGCCGGATATGCCTCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑CGTATAc                                                                                                                                                                                                   <  2:305476/239‑1 (MQ=255)                          aaaaTG‑GGGGTAGGCGGTACTATCCCCGATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGTATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑CGTATAc                                                                                                                                                                                                   <  1:140009/239‑1 (MQ=255)                           aaaTTGGGCGTAGGCGGTATTATCCCCTTCCCCCGCGCAAGTAAAACCCGGGATGTCATCGAAAAATTACGGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑CGTATAc                                                                                                                                                                                                   <  2:200564/239‑1 (MQ=255)                                 ggCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGGTTTCGCCGGGCTTTGCGAGTTGCCCCGGGAACCGGGGGGTGAGATTGCGCACGCGTTGGTtaaaa                                                                                                                                                                                                              >  2:185187/1‑221 (MQ=255)                                            tACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGGTGCTTAATAACA‑CGTATACGTTGTTTAACTTGCCGt                                                                                                                                                                                  >  1:74976/1‑239 (MQ=255)                                            tACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGGGCTTTGCGAGTTGGGCCGGGAACCGGCGGGTGGGATTGGGAAAGCGTTTGTTTATAACA‑CGGATACGTTGTTTAACGTGCCgg                                                                                                                                                                                  >  1:54810/1‑238 (MQ=255)                                               caaccccTACCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑CGTATACGTTGTTTAACGTGCCGTAgg                                                                                                                                                                               <  1:552505/234‑1 (MQ=255)                                                aTCCCCTACCCCCGCGCCAGTAAACCCCGCGATGTCCTCGCAAAATTCCCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑CGTATACGTTGTTTAACGTGCCGTAGGc                                                                                                                                                                              <  2:721635/239‑1 (MQ=255)                                                      tATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGAGgaga                                                                                                                                                                                                                                   >  1:485743/1‑181 (MQ=255)                                                      tATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTACCCGTGCTTTGCGAGTTGCCCCGGGAACCGGCGGGTGCGGTTGGCCAAGCGTTGCTTAATAACC‑CGGTACCGTTGTTTAACGTGCCGTAGGCCgggaa                                                                                                                                                                        >  2:845375/1‑236 (MQ=255)                                                      tATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGGGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGGTGGGGTTCGGCAAGCGTTGCTTACACACC‑CGTAAACGGTTTTTAACGTTCCAGAGGCCGGaaa                                                                                                                                                                        >  2:676232/1‑237 (MQ=255)                                                       aTCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGGTGGGATTGCGCAAGCGTTGCTTAATAAAA‑CGAATACGTTGGTTAACGTGCCGGAAGCCGGATaa                                                                                                                                                                       >  2:588896/1‑239 (MQ=255)                                                        tCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACACGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTTCGCAAGCGTTGCTTAAAAACC‑CAGATACGTTTTTTAAAGGGCCCTAGGCCGGAGAAg                                                                                                                                                                      >  2:357671/1‑239 (MQ=255)                                                             gcgcCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGCTGGGCTTGGGCAAGCGTTGCTTAATAACC‑CGTATACGTTGTTTAACGACCCGTAGGCCGGATAAGGGGtt                                                                                                                                                                 >  2:340360/1‑239 (MQ=255)                                                             gcgcCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAACCGTTGCTtaata                                                                                                                                                                                                              >  1:770039/1‑195 (MQ=255)                                                                         cccGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑CGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGGGTTCACGCCCCATcc                                                                                                                                                     >  1:837932/1‑239 (MQ=255)                                                                                   cATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑GGTATACGTTGTTTAACGTGCCGTAGGCCGGTTAAGGGGTTCACGCCGCATCCGGCCGTTgg                                                                                                                                            >  1:890167/1‑238 (MQ=255)                                                                                    aTCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑CGTATACGTTGTTTAACGTGCCGTGGCCCGGATAAGGCGTTCACGCCGCATCCGTAAGTTggg                                                                                                                                           >  1:847609/1‑237 (MQ=255)                                                                                                   cTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGGGATTGCGCCAGCGTTGCTTAATAACA‑CGTATACGTTGTTTAACGTGCCGTAGGCCGGCAAAAGAGTTTACGCCGCATCCGGCAGTTGGCCCAAAAAGCCCGagcg                                                                                                                           >  2:593809/1‑236 (MQ=255)                                                                                                       cATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGGTCGCCGTGCTTTGCGAGTTCCGCCGGGAACCGGCGGCTGAGATTGCGCAAACGTTGCTTAATAACA‑CGACTACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCCCCCCGCATCCGGCAGTGTGCGGACAAGCCCTGAGGCGAc                                                                                                                        >  2:712216/1‑238 (MQ=255)                                                                                                             tATTGCTGGAAACCGACGCGCCGGATAGGCCGCTCAACGGTTTGCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑CGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTaa                                                                                                                  <  1:483514/238‑1 (MQ=255)                                                                                                              aTTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACC‑CGTATACGTTGTTTAACTTGCCGTAGGCCGGGTAAGGCGTTCACGCCGCATCCGGCCGTTGGCGCCCAAAACCTGATGCGACCCTTaaac                                                                                                                >  2:705373/1‑237 (MQ=255)                                                                                                                       acccgccgCGCCGGATAGGCCGCTCAACGGTTTTCGGGGGCCGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑CGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATc                                                                                                       <  1:89728/237‑1 (MQ=255)                                                                                                                         ccccgcggcgtgATATGCCGCTCAAGGGTTTTCGGGGGGAGCCTAACCGCCCGGGGCAGGCTGGCGGTGTGTTCGCCGTGCTTTGCGAGTGGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑CGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCAt                                                                                                     <  2:345603/228‑1 (MQ=38)                                                                                                                            aCGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑CGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCCCAATGCCTGATGCGACGCTTAACGCGTCTTATCATGc                                                                                                   >  1:648434/1‑238 (MQ=255)                                                                                                                              gcccaggttatGCCCCCCAACGGTTTTCGGGGGCAGCCAAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑CGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTa                                                                                                <  1:705373/231‑1 (MQ=255)                                                                                                                                          cccccacagggTTTTCGGGGGCAGCCTAACCGCCCGGGGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑CGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGcc                                                                                    <  1:263262/230‑1 (MQ=255)                                                                                                                                              tCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑CGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCCCAATGCCCGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGaa                                                                                 >  1:824552/1‑238 (MQ=255)                                                                                                                                               cAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑CGTATACGTTGTTTAACGTGCCGTAGGCCGGCTAAGGCGTTCAGGCCGCATCCGGCAGTTGCCGCCAAATGCCTGATGCGACGCTTAACGCGTCTTATCATCCCTACAGGTTTTGGCCGAAcc                                                                               >  1:758501/1‑239 (MQ=255)                                                                                                                                                   gtgtttCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCGGTGTGTTCCCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑CGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAg                                                                           <  1:593809/236‑1 (MQ=255)                                                                                                                                                                      aCCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCCTTGCTTAATAACA‑CGTATACGTTGGTTAACCTGCCGTAGGCCGGGTAAGGCGTTCACGCCGCCTCCCGCAGTTTGCGGCCCATGCCTGATTGGACGCTTAACGCGCCTTATCCTTCCTCACGGTTTTTTCCCAACCGGACGGCGGGTAAGCCGGTCAc                                                         >  2:345057/1‑238 (MQ=255)                                                                                                                                                                            cGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑CGTATACGTTGTTTAACGTGCCGTAGGCCGGCTAAGGCGCTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTCGGCAGGATAAGGGGTGCCAGCCGCa                                                   >  1:342217/1‑238 (MQ=37)                                                                                                                                                                            cGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑CGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGtttcagcccca                                                   >  1:854930/1‑229 (MQ=255)                                                                                                                                                                                  gggCTGCCCGGGTGTTGGCGGTGTTTTGCGGGTGGCGCGGGGACCCGGCGGATGGGATTGCGCAAGCGTTGCTTAATAACC‑CGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGc                                             <  2:648434/237‑1 (MQ=21)                                                                                                                                                                                       ctcggggtgtTGCCCGTTTTTTGCGGGTTGCGGCGGGAACCGGGGGGTGGGATTGCCGAAGCGTGGGTTAA‑AACACCGGATGCGTGGTTTAACGGGCGGTGGGCCGGATAAGGGGTTCCCGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTgg                                       <  2:145257/233‑1 (MQ=17)                                                                                                                                                                                                   gCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑CGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcg                                                        <  1:249045/210‑1 (MQ=255)                                                                                                                                                                                                   gCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑CGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCAcg                                                        >  2:249045/1‑210 (MQ=255)                                                                                                                                                                                                            tGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑CGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGCCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCCCAATGCCTGATGcacc                   >  1:784392/1‑235 (MQ=255)                                                                                                                                                                                                                  ttggcccGGGACCCGGGGGGTGGGATTGGGCAAGGGTTTGTTAATAACC‑GGTATCCTGTGTTTAACGTGCCGGAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCAAAAGGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTc            <  2:877008/234‑1 (MQ=17)                                                                                                                                                                                                                  ttGGGCCGGGAACCGGGGGATGGGATTGCGCAAGCGTTGCTTAATAACA‑CGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGt             <  1:340360/238‑1 (MQ=255)                                                                                                                                                                                                                   tGCGCCGGGAACCGGGGGATGAGATTGCGCAAGCGTTGCTTAATAACC‑CGTATCCGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTc            <  1:845375/238‑1 (MQ=255)                                                                                                                                                                                                                   tGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑CGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCttttgc           >  1:89562/1‑235 (MQ=255)                                                                                                                                                                                                                             aCCGGCGGATGGAATGCCGTAAGCGGTGTTTAATACCA‑CGTAAACGTTGTTTAACGTGCCGTAGGCTGGATAAGTCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCTTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATc  <  1:345057/238‑1 (MQ=17)                                                                                                                                                                                                                             aCCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑CGTATACGGTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATc  <  1:227303/238‑1 (MQ=255)                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                    CGAACGGTTTGTACAGCTGGGCTACAAAATT‑GGCGTAGGCGGTACTATCACCTATCCACGCGCCAGTAAAACCCGCGATGTCATCGCAAAATTACCGCTGGCATCGTTATTGCTGGAAACCGACGCGCCGGATATGCCGCTCAACGGTTTTCAGGGGCAGCCTAACCGCCCGGAGCAGGCTGCCCGTGTGTTCGCCGTGCTTTGCGAGTTGCGCCGGGAACCGGCGGATGAGATTGCGCAAGCGTTGCTTAATAACA‑CGTATACGTTGTTTAACGTGCCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTAACGCGTCTTATCATGCCTACAGGTTTGTGCCGAACCGTAGGCCGGATAAGGCGTTCACGCCGCATCCGGCAGTTGGCGCACAATGCCTGATGCGACGCTTGTCGCGTCTTATC  >  NC_000913/4613978‑4614434

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