Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	4,630,385	G→A	S46S (TCC→TCT)	ettA ←	energy‑dependent translational throttle protein EttA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,630,385	0	G	A	89.5%	85.9 / 1.2	39	S46S (TCC→TCT)	ettA	energy‑dependent translational throttle protein EttA
Reads supporting (aligned to +/- strand):  ref base G (3/1);  new base A (25/9);  total (28/11)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.41e-01

TCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTT‑AACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGT‑GGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATAATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGCCGGTCGATCACCCTCGCAGT  >  NC_000913/4630137‑4630622                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              cccAGACGCCTTTGTTCAGCGCCCAGTGTGTAAAAACCGGCATCCGTAGCGGCGTACAGCGCGTACTCTTCACCCAGGTGTTGCAGGGCGTA‑ACCCCCTTCTGAAACCGCTTCTTCAAGGGACTCACGCAGGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATACCCAATCTTGATGTCTGGCTGCGGACGCGTTTCCCCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGT‑g                                                                                                                                                                                                                                               <  1:114207/248‑1 (MQ=255)     gACGGCTTTGTGCAGGGGCCAGCTTGTCAAAATCGCCAGCCGGATCGTCTTACAGCGCATACACTTCACCCAGGCGTTTCAGGGCTTT‑AACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCTTGCGGCAGGTACCCATTCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGT‑AGACTTa                                                                                                                                                                                                                                         <  1:260506/251‑1 (MQ=255)                     ggCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTT‑AACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGT‑AGACTTCCCCGCGCCGTTCAGAc                                                                                                                                                                                                                         >  1:301399/1‑251 (MQ=255)                     ggCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGAGCATACACTTCATCCAGGCGTTTCAGGGCGTT‑AACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGGAACCAATCTTGATATCGGGCTGCGGACGCGCTTCACCTTCAGTGTCTTTATCAATTCCCGCCATAATGCACAGCAGGAT‑AGACTTACCCGCGccgttaaga                                                                                                                                                                                                                          >  1:271356/1‑243 (MQ=255)                       ccAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTT‑AACCACTTCAGAAACCGCTTCTTCAATGGCCTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGGCGCGCTTTCCCTTCGTGGTCTTTATCCATGCCCGCCCTAATGCGCAGCAGGGG‑AGACTTTCCCGCCCCGTTCCGAccc                                                                                                                                                                                                                       >  2:169254/1‑251 (MQ=255)                           cTTGTCAAAATCGGCATCCGGATCGGAGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTT‑AACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGGAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCCCCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCCGGGT‑AGACTTACCCGCCCCGTTCAGcccccgg                                                                                                                                                                                                                    >  2:57670/1‑247 (MQ=255)                                         cATCCGGATCGGCGTACAGCGCATACACTTCATCCACACGTTTCAGGGCGTT‑AACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCCCGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCCATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGT‑GGCCTTACCCGCGCCATTCAGACCCAGGACACCACTTTTTGcc                                                                                                                                                                                                     >  1:213217/1‑251 (MQ=255)                                              ggATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTT‑AACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTTCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGT‑AGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAgg                                                                                                                                                                                                >  1:119721/1‑251 (MQ=255)                                              ggATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTT‑AACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGT‑AGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAgg                                                                                                                                                                                                >  1:162678/1‑251 (MQ=255)                                                      gTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTT‑AACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGGAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATTAATGCCCGCCATAATGCGCAGCAGGGT‑AGGCTTACCCGCGCCGTTCAGACCCAGGGCCCCCATTTTTGTCCCAGGGGAGaaa                                                                                                                                                                                         >  2:137423/1‑250 (MQ=255)                                                                    cTTCATCCAGGCGTTTCAGGGCGTT‑AACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGT‑AGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGttt                                                                                                                                                                          >  1:181490/1‑251 (MQ=255)                                                                           ccggtcgttCAGGGCGTTTAACCACTTAAGAAACCCCTTCTTCAAGGCCCCCCCGCAGGTTGTGTCCCGGGTGCAGCTGCGGTTCCTGCGGCAGGTAACCAATTTTGATGTCGGGCGGCGGACGCGCTTCCCCTTGGATGTCTTTACAAATGCCCGCCATAATGCGCAGCAGGGT‑AGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCa                                                                                                                                                                      <  1:82740/242‑1 (MQ=255)                                                                             ggCGTTTCAGGGCGTT‑AACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATTTCTGGCTGCGGACGCGCTTCACCTTCGCTGGCTTTTTCAATGCCCGCCCTAATGCGCAGCAGGGT‑AGACTTAACCGCGCCCTTCCGACCCACGACACCAATTTTTGCCCCCGGGAAGAAACTACGAGGGATGTTTTTccacat                                                                                                                                                                  >  2:49148/1‑245 (MQ=255)                                                                                    cAGGGCGTT‑AACCACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATTCCCGCCATAATGCGCAGCAGGGT‑AGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGtt                                                                                                                                                          >  1:506566/1‑251 (MQ=255)                                                                                     aGGGCGTT‑AACCACTTATGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCTTGCGGCAGATAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGT‑GGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGttt                                                                                                                                                         >  2:14570/1‑251 (MQ=255)                                                                                         gggttaccctcttcagaaactcctccttcAATGGACCCAGCCACGGTTTTTTCCGGGTTCATTGGCGTTTCCTGCGGCAGTTAACCATTCTTGTTGTCTTGCTGCGGCCGCGCTTCACCTTGGTTGTCTTTATCAATCCCCGCCATAATGCGCAGCAGGGT‑AGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTcgg                                                                                                                                                      <  1:348520/226‑1 (MQ=255)                                                                                                 cACTTCAGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCCATCTTGATGTCTGGCTGCGGCCGCGCTTCACCTTCGATGTCTTTATCAATTCCCGCCCTAATGCGCAGCAGGGT‑AGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGacca                                                                                                                                               >  2:270314/1‑247 (MQ=255)                                                                                                             cGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGCAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGT‑AGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGCAGAAAACCGGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTTCCggacc                                                                                                                                  >  2:438606/1‑247 (MQ=255)                                                                                                                            cTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCACGTTCCCAATCTTGTTGTTTGGCTCCGGGCGCGCTTCACCTTCGCTGACTTTATCAATGCCCGCCATACTGCGCAGCAGGGT‑AGACTTACCCGCGACGTTCAAGCCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGAAGTTTTTCAAAATATTACGTTTCGGCGGAACAAACTTGCCGGCACGATGCATGGTATaa                                                                                                                    >  2:46262/1‑250 (MQ=255)                                                                                                                                  cacgGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGGAACCAATCTTGATGTCTGGCTGCGGGCGCGCTTCACCTTCGTTGTCTTTATAAATGCCCGCCATAATGCGCAGCAGGGTAGGA‑TTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTATGAGAGATGTTTTTCAAAAATTGACGTTTCGGCGGAACAAATTTGGCGAAACGAAGCATGGTATAAAcagagt                                                                                                             >  2:305138/1‑246 (MQ=255)                                                                                                                                     ggTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGGTAACCAATCTTGATGGCTGGCTGCGGACGCGCCTCCCCTTCGCTGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGG‑AGACTTCCCCGCGCCGTTCAGCCCCAGGCAACCAATTTTTGCCCCAGGGAAGACACTCAGAGAGATGTTTTTCAAAATATGTCGCTTCGGCGGAAAAACTTTGCCGTCACGCTGCATGGAATAAACCAAtttgg                                                                                                          >  2:183563/1‑248 (MQ=255)                                                                                                                                                   ccgtcTGCGGTTCCTGCGGCAGGTACCCATTCTTGATTTCTGGCTGCGGGCGCTCTTCCCCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGTGT‑AGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACtt                                                                                             <  1:466462/246‑1 (MQ=255)                                                                                                                                                           ggTTCTTGCGGCAGATACCCAATCTTGATGTCTGGCTGCGGACGCGCTTCCCCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGT‑GGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTAt                                                                                    <  1:14570/251‑1 (MQ=255)                                                                                                                                                           ggTTCCTGCGGCAGGTAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGT‑AGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGGATAAACGAATTGATGCACGTTGGACTTCGCCTCTAt                                                                                    >  1:281841/1‑251 (MQ=255)                                                                                                                                                                      caggtAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGT‑AGACTTACCCGCGCCGTTCAGACCCAGGCCACCAATTTTTGCCCCAGGGAAGAAACTCAGCGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACttt                                                                                                                                          >  2:279491/5‑186 (MQ=255)                                                                                                                                                                          tAACCAATCTTGGTGTCTGGCTGCGGACGCGCTTAACCTTCGATGTCTTTATCAATGCCCGCCATAATGACCAGCAGGGC‑AGACTTCCCCGCGCCGTTCAGACCCAGGACACCAAATTTTGCCCCAGGGAAGAAACTCCGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCCACACGATACATGGTATAAACGAATTAAGCCACGTTGGCCTTCGCCTCTCTGTTGATCGTGAtgat                                                                     >  1:320683/1‑248 (MQ=255)                                                                                                                                                                                                        gcTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGT‑AGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCttttt                                       >  1:499609/1‑251 (MQ=255)                                                                                                                                                                                                          ttCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGT‑AGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTcc                                     >  1:442337/1‑251 (MQ=255)                                                                                                                                                                                                            ccccTTCGATGTTTTTATCAATTCCCGCCATAATGCGCAGCAGGGT‑AGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTccc                                    <  1:270314/248‑1 (MQ=255)                                                                                                                                                                                                                 tCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGT‑AGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGACGTGTAGCCTTTTTCCCCGCCt                               >  2:174816/1‑250 (MQ=255)                                                                                                                                                                                                                      gTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGT‑AGACTTACCCGCGCCGTTCAGCCCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCCTGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTTATGATGAGTTTTCCAAGGCGAAATGAAGCCTTTTTCCCCGCCTAATccc                         >  2:187957/1‑251 (MQ=255)                                                                                                                                                                                                                      gTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGT‑AGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATccc                         >  1:183395/1‑251 (MQ=255)                                                                                                                                                                                                                          tcacccagccccGCCAAAATTCTCATCAGGGT‑AGACTTCCCCGCGCCGTTCAGACCCAGGACCCCAATTTTGGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTTCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTTTGTTTATCGTTATGATGAGTTTTCAAAGGCGAAGTGTCGCCTTTTTCCCCGCCTAATCCCAGc                      <  2:67634/241‑1 (MQ=255)                                                                                                                                                                                                                           tatcaagcccctccATAATGCGCAGCAGTGT‑AGATTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTCCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAgtt                     <  1:447576/238‑3 (MQ=255)                                                                                                                                                                                                                              cAATGCCCGCCATAATGCGCAGCAGGGT‑AGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATccccgttggtc                 >  2:156202/1‑243 (MQ=255)                                                                                                                                                                                                                                 cccccGCCATAATGCTCAGCAGGGT‑AGACATACTCGCGCCGTTCAGACCCAGGACACAAATTTTTGCCCCAGGGAAGAAAACCAGAGAGATGTTTTTCAAAATATGACGTGTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGAt              <  2:356292/249‑1 (MQ=255)                                                                                                                                                                                                                                  ccccGCCATAATGCGCAGCAGGGG‑AGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATAGGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTGTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATc             <  1:57670/250‑1 (MQ=255)                                                                                                                                                                                                                                    cagagatcaTGCGCAACAGGGT‑CGATTGACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGATATGTGTTTCAAAATATGACGTTTcggcgg                                                                                                                                                   <  2:36434/107‑1 (MQ=255)                                                                                                                                                                                                                                        cATAATGCGCAGCAGGGT‑AGACTTACTCGCGCCGTTCAGACCCAGCACACCAATTTTTGCCCCAGGGCAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGCACCACTTTGCCGACACGATGCATGGTCTAAACGAATTGAGCCACGTTGGACTTCGCCTCTCTGTTTATCGTGATGATGGGTTTTCACAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATccacgttggtcgttcactttc       >  2:62715/1‑232 (MQ=255)                                                                                                                                                                                                                                             tGCGCAGCAGGGT‑AGACTTACCCGCGCCGTTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATGATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGTTGGTCGATCACCTTCGCagt  >  1:184363/1‑251 (MQ=255)                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                              TCCAGACGGCCTTGTTCAGCGGCCAGCTTGTCAAAATCGGCATCCGGATCGGCGTACAGCGCATACACTTCATCCAGGCGTTTCAGGGCGTT‑AACCACTTCTGAAACCGCTTCTTCAATGGACTCACGCACGGTGTGTTCCGGGTTCAGCTGCGGTTCCTGCGGCAGATAACCAATCTTGATGTCTGGCTGCGGACGCGCTTCACCTTCGATGTCTTTATCAATGCCCGCCATAATGCGCAGCAGGGT‑GGACTTACCCGCGCCATTCAGACCCAGGACACCAATTTTTGCCCCAGGGAAGAAACTCAGAGAGATGTTTTTCAAAATATGACGTTTCGGCGGAACAACTTTGCCGACACGATGCATGGTATAAACGAATTGAGCCACGTTGGACTTCGCCTCTATGTTTATCGTGATAATGAGTTTTCAAAGGCGAAGTGTAGCCTTTTTCCCCGCCTAATCCCAGCCGGTCGATCACCCTCGCAGT  >  NC_000913/4630137‑4630622

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