Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,305,274	G→T	H442N (CAT→AAT)	yjcS ←	alkyl sulfatase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,305,274	0	G	T	99.0%	331.8 / ‑3.1	95	H442N (CAT→AAT)	yjcS	alkyl sulfatase
Reads supporting (aligned to +/- strand):  ref base G (1/0);  new base T (45/49);  total (46/49)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.35e-01

AATCACCTGTTTCAGCAGTTCTGCCGACCAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCG  >  NC_000913/4305137‑4305407                                                                                                                                          |                                                                                                                                      aaTCACCTGTTTCAGCAGTTCTGCCGACCAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCa                                                                                                                                     >  1:1295965/1‑140 (MQ=255)   tCACCTGTTTCAGCAGTTCTGCCGACCAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTTCCCATCTCCACCTGACCATACGGATTCAgg                                                                                                                                   <  2:1585877/140‑1 (MQ=255)      ccTGTTTCAGCAGTTCTGCCGACCAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTa                                                                                                                                <  2:2113081/140‑1 (MQ=255)       cTGTTTCAGCAGTTCTGCCGACCAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAg                                                                                                                               <  1:986763/140‑1 (MQ=255)        tGTTTCAGCAGTTCTGCCGACCAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGc                                                                                                                              <  2:2112178/140‑1 (MQ=255)             cagcagTTCTGCCGACCAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGt                                                                                                                         <  2:154915/140‑1 (MQ=255)             cagcagTTCTGCCGACCAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGt                                                                                                                         <  2:706045/140‑1 (MQ=255)               gcagTTCTGCCGACCAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGtt                                                                                                                                 <  1:1796615/130‑1 (MQ=255)               gcagTTCTGCCGACCAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGtt                                                                                                                                 >  2:1796615/1‑130 (MQ=255)               gcagTTCTGCCGACCAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTa                                                                                                                       <  1:603431/140‑1 (MQ=255)               gcagTTCTGCCGACCAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTa                                                                                                                       >  2:878434/1‑140 (MQ=255)                cagTTCTGCCGACCAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTAc                                                                                                                      <  1:1281394/140‑1 (MQ=255)                  gTTCTGCCGACCAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACcg                                                                                                                    <  2:1758604/140‑1 (MQ=255)                  gTTCTGCCGACCAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACcg                                                                                                                    <  2:1429628/140‑1 (MQ=255)                   ttCTGCCGACCAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACcgt                                                                                                                   >  1:1889862/1‑140 (MQ=255)                    tCTGCCGACCAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACcgtc                                                                                                                  >  2:332966/1‑140 (MQ=255)                    tCTGCCGACCAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACcgtc                                                                                                                  >  2:78015/1‑140 (MQ=255)                     cTGCCGACCAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACcgtcg                                                                                                                 >  2:931392/1‑140 (MQ=255)                         cGACCAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAAt                                                                                                             <  1:1370408/140‑1 (MQ=255)                         cGACCAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAAt                                                                                                             <  1:332966/140‑1 (MQ=255)                            ccAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGc                                                                                                          <  1:931392/140‑1 (MQ=255)                             cAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGcc                                                                                                         <  2:591782/140‑1 (MQ=255)                             cAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGcc                                                                                                         <  1:1182817/140‑1 (MQ=255)                              aGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCa                                                                                                        <  1:702585/140‑1 (MQ=255)                                  gTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGTGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGAt                                                                                                    <  1:1014832/140‑1 (MQ=255)                                  gTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGc                                                                                                          >  2:1706677/1‑134 (MQ=255)                                  gTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGc                                                                                                          <  1:1706677/134‑1 (MQ=255)                                  gTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGAt                                                                                                    <  1:2153878/140‑1 (MQ=255)                                  gTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGAt                                                                                                    <  1:575366/140‑1 (MQ=255)                                  gTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGAt                                                                                                    <  2:552774/140‑1 (MQ=255)                                    aaTCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAg                                                                                                  <  2:348503/140‑1 (MQ=255)                                    aaTCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAg                                                                                                  <  1:495493/140‑1 (MQ=255)                                    aaTCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAg                                                                                                  <  1:526878/140‑1 (MQ=255)                                    aaTCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAg                                                                                                  <  1:991624/140‑1 (MQ=255)                                    aaTCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAg                                                                                                  <  1:878434/140‑1 (MQ=255)                                    aaTCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAg                                                                                                  <  1:1905517/140‑1 (MQ=255)                                    aaTCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAg                                                                                                  <  2:1226272/140‑1 (MQ=255)                                    aaTCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTACACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAg                                                                                                  <  2:1235257/140‑1 (MQ=255)                                         ccTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGtt                                                                                             >  2:559638/1‑140 (MQ=255)                                          cttgcttgTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGc                                                                                                          >  2:883283/1‑126 (MQ=255)                                          cttgcttgTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGc                                                                                                          <  1:883283/126‑1 (MQ=255)                                          cttgcttgTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAg                                                                                                  <  1:147147/134‑1 (MQ=255)                                          cttgcttgTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAg                                                                                                  >  2:147147/1‑134 (MQ=255)                                          cttgcttgTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTta                                                                                            >  1:1595851/1‑140 (MQ=255)                                          cttgcttgTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTta                                                                                            >  1:1095846/1‑140 (MQ=255)                                          cttgcttgTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTta                                                                                            >  1:1404399/1‑140 (MQ=255)                                             gcttgTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAg                                                                                                  <  2:1156341/131‑1 (MQ=255)                                             gcttgTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAg                                                                                                  >  1:1156341/1‑131 (MQ=255)                                                       ttcttGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCggg                                                                               <  2:1095846/140‑1 (MQ=255)                                                            gcgcCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTtg                                                                          <  1:2097531/140‑1 (MQ=255)                                                            gcgcCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTtg                                                                          <  2:208322/140‑1 (MQ=255)                                                            gcgcCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTtg                                                                          >  1:545434/1‑140 (MQ=255)                                                            gcgcCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTtg                                                                          <  2:1404399/140‑1 (MQ=255)                                                                 aGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCgg                                                                                <  2:115148/129‑1 (MQ=255)                                                                 aGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCgg                                                                                >  1:115148/1‑129 (MQ=255)                                                                   gTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGa                                                                   >  2:919986/1‑140 (MQ=255)                                                                        tgaCACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGaa                                                              >  1:13432/1‑140 (MQ=255)                                                                        tgaCACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGaa                                                              >  1:1971627/1‑140 (MQ=255)                                                                            acGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGt                                                          >  1:1771636/1‑140 (MQ=255)                                                                            acGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGt                                                          >  1:366353/1‑140 (MQ=255)                                                                              gggCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTaa                                                        >  2:1655386/1‑140 (MQ=255)                                                                              gggCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTaa                                                        >  1:372790/1‑140 (MQ=255)                                                                              gggCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTaa                                                        >  1:2091155/1‑140 (MQ=255)                                                                               ggCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAAt                                                       >  1:1164391/1‑140 (MQ=255)                                                                                    aaCCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCcg                                                  >  2:922506/1‑140 (MQ=255)                                                                                      ccGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCcgcg                                                >  1:1874894/1‑140 (MQ=255)                                                                                       cGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCgg                                               <  1:559638/140‑1 (MQ=255)                                                                                        gCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGc                                              <  1:922506/140‑1 (MQ=255)                                                                                        gCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGc                                              <  1:1400037/140‑1 (MQ=255)                                                                                        gCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGAGGCTGACAGAACCGTAATAGCCGCGGc                                              <  1:78015/140‑1 (MQ=255)                                                                                         cccAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCt                                             >  1:2020226/1‑140 (MQ=255)                                                                                             gcgcCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGcc                                         <  1:919986/140‑1 (MQ=255)                                                                                              cgcCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGccc                                        >  2:1068596/1‑140 (MQ=255)                                                                                                ccTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAg                                      >  1:881395/1‑140 (MQ=255)                                                                                                    cACGTAACGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGTCAGAACCGTAATAGCCGCGGCTGGCCCAGTTAt                                  >  1:1977848/1‑140 (MQ=255)                                                                                                          aCGTTTACCCATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCaa                            <  2:1889862/140‑1 (MQ=255)                                                                                                               tACCCATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCa                       >  1:794257/1‑140 (MQ=255)                                                                                                                   cATCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAggcg                   <  2:1164391/140‑1 (MQ=255)                                                                                                                    aTCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAggcgg                  <  2:372790/140‑1 (MQ=255)                                                                                                                    aTCTCCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAggcgg                  <  2:545434/140‑1 (MQ=255)                                                                                                                       tcCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAg               <  2:1874894/140‑1 (MQ=255)                                                                                                                       tcCACCTGACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAg               <  2:2020226/140‑1 (MQ=255)                                                                                                                              gACCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATc        >  2:1212692/1‑140 (MQ=255)                                                                                                                               aCCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCa       >  1:756821/1‑140 (MQ=255)                                                                                                                               aCCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCa       >  2:791039/1‑140 (MQ=255)                                                                                                                               aCCATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCa       >  2:1293761/1‑140 (MQ=255)                                                                                                                                ccATCCGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCAt      >  2:1337790/1‑140 (MQ=255)                                                                                                                                ccATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCAt      >  2:912992/1‑140 (MQ=255)                                                                                                                                ccATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCAt      >  2:359308/1‑140 (MQ=255)                                                                                                                                 cATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATg     >  2:804620/1‑140 (MQ=255)                                                                                                                                 cATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATg     >  2:1406445/1‑140 (MQ=255)                                                                                                                                 cATACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATg     >  1:1405475/1‑140 (MQ=255)                                                                                                                                   tACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTc   >  2:620831/1‑140 (MQ=255)                                                                                                                                   tACGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTc   >  1:772290/1‑140 (MQ=255)                                                                                                                                    aCGGATTCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCg  >  1:818642/1‑140 (MQ=255)                                                                                                                                          |                                                                                                                                      AATCACCTGTTTCAGCAGTTCTGCCGACCAGCGGTAATCACCTTGCTTGTTCGCTTCTTGCGCCAGGTTGATGACACGGGCAGAACCGCCCAGCGCCTGCACGTAACGTTTACCCATCTCCACCTGACCATACGGATGCAGGTTAGCCGGGTTACCGTCGTAATAGCCAAGATAGAAGTTATACACCGCGCGGGCGTTGTGGCTGACAGAACCGTAATAGCCGCGGCTGGCCCAGTTATTGGCAAGTGCAGGCGGCAGCTTAATCATGTCG  >  NC_000913/4305137‑4305407

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