Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,449,865	G→A	P7S (CCT→TCT)	rpsC ←	30S ribosomal subunit protein S3

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,449,865	0	G	A	100.0%	137.2 / NA	40	P7S (CCT→TCT)	rpsC	30S ribosomal subunit protein S3
Reads supporting (aligned to +/- strand):  ref base G (0/0);  new base A (14/26);  total (14/26)

CAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGTTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTT  >  NC_000913/3449623‑3450109                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     cAGTTTTTCTACGTCTTCACCTTTTTTGCCGATAACGATACCCGTGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACTGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCGCTGTCCAGGTTGTCAGTGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACtt                                                                                                                                                                                                                                              <  1:383712/251‑1 (MQ=255)      tttCTACGTCTTCACCTTTTTTCCCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTg                                                                                                                                                                                                                                          <  2:221556/250‑1 (MQ=255)        tCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGAc                                                                                                                                                                                                                                        <  1:362937/250‑1 (MQ=255)          tACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGAcccc                                                                                                                                                                                                                                     >  2:211582/1‑250 (MQ=255)                 tCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTa                                                                                                                                                                                                                               <  1:288676/250‑1 (MQ=255)                  cACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGt                                                                                                                                                                                                                             <  1:76990/251‑1 (MQ=255)                              cGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGtctc                                                                                                                                                                                                                 >  1:36447/1‑251 (MQ=255)                               gATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGtctc                                                                                                                                                                                                                 >  1:65336/1‑250 (MQ=255)                                  aaCGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGc                                                                                                                                                                                                              <  1:116899/250‑1 (MQ=255)                                     gATACCCGGGCGAGCGGTGTGAATGGTTACACGGATGCTCTTAGCCCGACGCTCGATAACGATACTAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCGCTGTCCAGCTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATTGTTTTACAATTCCCGGGCGAATACCATTAGAATGTACTTTCTGCCCCATTGCTAGTCTCCAGAGTCTCAGCGATc                                                                                                                                                                                                          <  2:64373/251‑1 (MQ=255)                                      aTACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCg                                                                                                                                                                                                         >  2:242654/1‑251 (MQ=255)                                                   cAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACcac                                                                                                                                                                                              >  1:120991/1‑249 (MQ=255)                                                                  cacGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTgc                                                                                                                                                                             <  2:221702/251‑1 (MQ=255)                                                                        tgcaattAGACGGACGCTCGATAACGATACGAGATACGGACGCGTTAGCCAGTTCCTTAGGCAGGTACTGACGTACTTTAAAAGCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGTAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGACGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTc                                                                                                                                                                        <  2:112632/245‑1 (MQ=255)                                                                              tAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTATAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGTCTGGTGCGCTTCAgta                                                                                                                                                                    <  1:88840/248‑3 (MQ=255)                                                                                               aCGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCtttt                                                                                                                                                >  1:136306/1‑251 (MQ=255)                                                                                                 gATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGGACTGACGTACTTTAAAATCGATGTCCCGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTAGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGc                                                                                                                                              >  1:235151/1‑251 (MQ=255)                                                                                                      gagaTTCGTACGCTTTATCCAGTTCCTTAGTCATGTACTGACGTACTTTAAAATCGCTGTCCAGTTTTTCAGCGAATTCTTTGTTGTTCGCAAACCAGGTAGAGTTCCATTGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCATAGTCTCAGCGATCGGACACAACCACAGTGATGTTGCTGGTGCGCTTCAGTATGCGATCTGCACGCCCTTTTGCAcgcg                                                                                                                                         <  1:33030/251‑1 (MQ=255)                                                                                                        gaTACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGc                                                                                                                                       <  2:120991/251‑1 (MQ=255)                                                                                                            cGGACGCTTTAGCCAGTTCCTTAGTCATGTACTGACGTACTTTAAAATCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCAGATTGATGTGGCGGTTGCTCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATaa                                                                                                                                   <  1:18588/251‑1 (MQ=255)                                                                                                                    ttAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTc                                                                                                                           >  1:284581/1‑251 (MQ=255)                                                                                                                     tAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCTCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGTTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTc                                                                                                                           <  1:47949/250‑1 (MQ=255)                                                                                                                                   gTCAGGTACTGACGTACTTTAAAATCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGAACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCg                                                                                                            >  1:295387/1‑251 (MQ=255)                                                                                                                                       ggTACTGACGTACTTTAAAATCGCTTTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTa                                                                                                        <  1:415580/251‑1 (MQ=255)                                                                                                                                                 tACTTTAAAATCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGGGCGCTTCAGGATGCGAACTGCACGACCTTTTGCACGCGGCATAAAGCGCTTCATGCTCGGGCCTTCGTCTACGAAAAtttt                                                                                              >  2:13657/1‑251 (MQ=255)                                                                                                                                                             gCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATCCCATTAGAATGTACTTTCTGACCTATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAg                                                                                  <  1:151418/251‑1 (MQ=255)                                                                                                                                                                gTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATc                                                                               <  1:24976/251‑1 (MQ=255)                                                                                                                                                                  ccAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGTTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGt                                                                             <  2:175016/251‑1 (MQ=255)                                                                                                                                                                           cAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCa                                                                     <  2:239468/250‑1 (MQ=255)                                                                                                                                                                                         gtgtTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTgtgt                                                      <  2:53647/251‑1 (MQ=255)                                                                                                                                                                                                   aCCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATTCCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTGCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCAGCGTTGTGTTCAGCGTTAg                                            <  2:59583/251‑1 (MQ=255)                                                                                                                                                                                                     cAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGTTCAGCGTTAGCa                                          <  1:199476/251‑1 (MQ=255)                                                                                                                                                                                                      aGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGTTCAGCGTTAGCa                                          >  1:295565/1‑250 (MQ=255)                                                                                                                                                                                                               tCCATGGTTTTACAATTCCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACTCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCTATCGTTGTGTTCAGCGTTAGCAATGGCAGAt                                 <  1:173020/250‑1 (MQ=255)                                                                                                                                                                                                                ccATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGTTCAGCGTTAGCAATGGCAGAtt                                <  1:315831/250‑1 (MQ=255)                                                                                                                                                                                                                 cATGGTTTTACAATACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCCTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGTTCAGCGTTAGCAATGGCAGATTcc                              >  2:393460/1‑251 (MQ=255)                                                                                                                                                                                                                  aTGGTTTTACAATATCCAGGCGAATATCATTATAAGGTACTTTCTGACCCATTGCTAGGCTCCAGAGTCTCAGAGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGGCCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGTTCAGCGTTAGCAATGGCAGATTCCa                             <  2:284581/251‑1 (MQ=255)                                                                                                                                                                                                                             aTACCCAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGTTCAGCGTTAGCAATGGAAGATTCCAGAACTTTCTTg                  >  1:60391/1‑251 (MQ=255)                                                                                                                                                                                                                                  cAGGCGAATACCATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGTTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAg             <  1:179542/251‑1 (MQ=255)                                                                                                                                                                                                                                             cATTAGAATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGTTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGATCAGTACAGCCGCtt  >  2:112214/1‑251 (MQ=255)                                                                                                                                                                                                                                                   |                                                                                                                                                                                                                                                     CAGTTTTTCTACGTCTTCACCTTTTTTACCGATAACGATACCCGGGCGAGCAGTGTGAATGGTTACACGGATGCTCTTAGCCGGACGCTCGATAACGATACGAGATACGGACGCTTTAGCCAGTTCCTTAGTCAGGTACTGACGTACTTTAAAATCGCTGTCCAGGTTGTCAGCGAATTCTTTGGTGTTCGCAAACCAGGTAGAGTTCCATGGTTTTACAATACCCAGGCGAATACCATTAGGATGTACTTTCTGACCCATTGCTAGTCTCCAGAGTCTCAGCGATCGGACACAACCACAGTGATGTGGCTGGTGCGCTTCAGGATGCGATCTGCACGACCTTTTGCACGCGGCATAATGCGCTTCATGCTCGGGCCTTCGTCTACGAAAATTTTCGTAACTTTCAGATCGTCAATGTCAGCGCCATCGTTGTGTTCAGCGTTAGCAATGGCAGATTCCAGAACTTTCTTGACCAGTACAGCCGCTT  >  NC_000913/3449623‑3450109

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