Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	2,531,483	T→A	S7S (TCT→TCA)	cysZ →	sulfate:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,531,483	0	T	A	100.0%	188.7 / NA	51	S7S (TCT→TCA)	cysZ	sulfate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (24/27);  total (24/27)

TCTCTGTTGTTCTAACACCTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCTGTGCTGTTAGTGTTTGGC  >  NC_000913/2531239‑2531729                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       tctcTGTTGTTCTAACACCTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGcccc                                                                                                                                                                                                                                                   <  2:692389/250‑1 (MQ=255)    ctGTTGTTCTAACACCTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCAc                                                                                                                                                                                                                                                 <  1:157248/249‑1 (MQ=255)       ttgttCTAACACCTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGc                                                                                                                                                                                                                                            >  1:69523/1‑251 (MQ=255)           tCTAACACCTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGt                                                                                                                                                                                                                                         <  1:621740/250‑1 (MQ=255)                     tGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGTGGTTTTTACTAt                                                                                                                                                                                                                                <  1:499837/249‑1 (MQ=255)                              ggCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAgg                                                                                                                                                                                                                     <  2:93548/251‑1 (MQ=255)                                cAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCt                                                                                                                                                                                                                   <  2:103590/251‑1 (MQ=255)                                    aaaagggttaaatccattaacagCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGaa                                                                                                                                                                                                               <  1:278147/229‑1 (MQ=255)                                     aTTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGa                                                                                                                                                                                                                >  2:337135/1‑249 (MQ=255)                                         aCTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTc                                                                                                                                                                                                           <  2:383864/250‑1 (MQ=255)                                                      aTTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTgg                                                                                                                                                                                             >  1:406962/1‑251 (MQ=255)                                                          gCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGAtt                                                                                                                                                                                         <  1:356780/251‑1 (MQ=255)                                                           ccGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTc                                                                                                                                                                                        <  1:592643/251‑1 (MQ=255)                                                              aTTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGc                                                                                                                                                                                     >  1:463214/1‑251 (MQ=255)                                                                 aCCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGtt                                                                                                                                                                                  >  1:397463/1‑251 (MQ=255)                                                                    tCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCg                                                                                                                                                                               >  1:163324/1‑251 (MQ=255)                                                                                     aTTAACTTTCACACCACAAAGATAGATGAAATCGTGCCTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCACTTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGGAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGCTTCGTTATTTTACCGctgctg                                                                                                                                                              >  2:93869/1‑251 (MQ=255)                                                                                       tAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTgg                                                                                                                                                             <  2:397463/250‑1 (MQ=255)                                                                                                    acaAAGAGAGATGAAATCGTGCTTTTTGATGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTtggttg                                                                                                                                               >  2:330958/1‑251 (MQ=255)                                                                                                     caAAGATAGATGAAAGCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTttgttgn                                                                                                                                              <  2:307442/251‑2 (MQ=255)                                                                                                                      gTGCTTTTTGCTGTTTTTTCTAACATATCCTAACTGTCCATTGCGCATTTCCCCGGTAAAATACGCAGAATTTTCCTGGGATTGGACAAAAGGAGCTCATCCTGCTTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCtgg                                                                                                                             >  2:614667/1‑251 (MQ=255)                                                                                                                      gTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCgctgct                                                                                                                                                               <  1:572212/217‑1 (MQ=255)                                                                                                                      gTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCgctgct                                                                                                                                                               >  2:572212/1‑217 (MQ=255)                                                                                                                        gCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCtggtg                                                                                                                           >  2:702878/1‑251 (MQ=255)                                                                                                                                        tCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTc                                                                                                           >  2:108033/1‑251 (MQ=255)                                                                                                                                           aaCATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGAGTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGAt                                                                                                        <  2:652223/251‑1 (MQ=255)                                                                                                                                             cATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTTCTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGt                                                                                                      <  2:270326/251‑1 (MQ=255)                                                                                                                                                atCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTg                                                                                                   <  1:337135/251‑1 (MQ=255)                                                                                                                                                         gTCCATTGCGCAATTTCCCGGTAAAATACGCAGAATTTTCCTGGGAGTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGATCTTTACACAGCt                                                                                                            <  2:453740/233‑1 (MQ=255)                                                                                                                                                         gTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCt                                                                                                            >  1:453740/1‑233 (MQ=255)                                                                                                                                                         gTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGAc                                                                                          >  1:122894/1‑251 (MQ=255)                                                                                                                                                         gTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGAc                                                                                          >  2:415955/1‑251 (MQ=255)                                                                                                                                                           ccATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGAct                                                                                         <  1:389299/250‑1 (MQ=255)                                                                                                                                                                            ggTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCg                                                                       >  2:471088/1‑251 (MQ=255)                                                                                                                                                                               aaaaTACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGa                                                                     <  2:463214/250‑1 (MQ=255)                                                                                                                                                                                aaaTACGCAGAATTTTCCTGGGCTTGCGCAAAAGGAGCTCATCATGATTATGGTTTCATCATTCACATCAGTCCAACGCCGCGGTTTTTACTATTTTGAGCAAGGCTGGCCGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATTTTTTGTTGATGGTGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGCATGGATACCGACTCTCATGAGTTACGTTCCGGACt                                                                   >  2:602989/1‑251 (MQ=255)                                                                                                                                                                                 aaTACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCACCACGCAGCGGTTTTTACTCTTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGCCTCTCATGACTTACGTTCCGGACTg                                                                  >  2:204382/1‑251 (MQ=255)                                                                                                                                                                                       cAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTAc                                                             >  1:201021/1‑250 (MQ=255)                                                                                                                                                                                        aGAATTTTCGTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTAATCATTCACATAAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGCCTCTCATGAGTTACGTTCCGGACTGGCTACa                                                            >  1:425332/1‑250 (MQ=255)                                                                                                                                                                                                      aTTGGTCAAAAGGGGCACATCCTGATTATGGCTTCATCATTCACATCAGCCCCGCGCAGCTGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGGCTCGCATCCTGGGGTTCGGCGTTTCGGTATTGTGCCGCTGCTGGTCAATATTTTGGTGATGGGGGGCGCATTCTTGTGGCTCTTTACACAGCTCGTTGTCGAGATTCCGACTCTCATGAGTTACGTTCCGATTTGGCGACAATGGCTGAGTTAtgt                                             <  2:201868/251‑3 (MQ=255)                                                                                                                                                                                                         ggTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCtgt                                           <  2:593393/250‑1 (MQ=255)                                                                                                                                                                                                          gTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTCTCtgttg                                         >  2:489727/1‑251 (MQ=255)                                                                                                                                                                                                            cAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTTGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTtgtg                                       <  1:330958/251‑1 (MQ=255)                                                                                                                                                                                                                   aGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTg                                >  2:399169/1‑251 (MQ=255)                                                                                                                                                                                                                        aTCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGACGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGt                           <  1:19445/251‑1 (MQ=255)                                                                                                                                                                                                                             gATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCAt                        <  1:613210/249‑1 (MQ=255)                                                                                                                                                                                                                              aTTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCAtctc                     <  2:406962/251‑1 (MQ=255)                                                                                                                                                                                                                                   ggTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTCTTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCAGACTCTCATGGGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCTTCTCTGTGc                >  2:107025/1‑251 (MQ=255)                                                                                                                                                                                                                                           catTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTGGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCTGTGCTGTTAgtg        <  2:163324/251‑1 (MQ=255)                                                                                                                                                                                                                                              tCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCTGTGCTGTTAGTGtt      <  1:610490/250‑1 (MQ=255)                                                                                                                                                                                                                                                 caTCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCTGTGCTGTTAGTGTTTGGc  >  1:90501/1‑251 (MQ=255)                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                       TCTCTGTTGTTCTAACACCTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCTGTGCTGTTAGTGTTTGGC  >  NC_000913/2531239‑2531729

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