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%	144.6 / NA	42	S7S (TCT→TCA)	cysZ	sulfate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base T (0/0);  new base A (20/22);  total (20/22)

TTCTCTGTTGTTCTAACACCTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCTGTGCTGTTAGTGTTT  >  NC_000913/2531238‑2531726                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    ttCTCTGTTGTTCTAACACCTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGcccc                                                                                                                                                                                                                                                >  2:33757/1‑251 (MQ=255)                    cTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGACGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTAt                                                                                                                                                                                                                             >  1:492681/1‑251 (MQ=255)                       gCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTAttt                                                                                                                                                                                                                           <  1:681338/250‑1 (MQ=255)                          accacGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTgcg                                                                                                                                                                                                                       >  1:370473/1‑251 (MQ=255)                                     cATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGaa                                                                                                                                                                                                            <  2:522552/251‑1 (MQ=255)                                           cTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGACATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCCGCGGTTTTTCCTATTTTGCGCAAGGCTGGAAGCTCGt                                                                                                                                                                                                      >  2:83377/1‑251 (MQ=255)                                                         ttGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCGCGTCTCGCAACCTggg                                                                                                                                                                                         <  2:896334/250‑1 (MQ=255)                                                           gCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGAt                                                                                                                                                                                       <  1:33757/250‑1 (MQ=255)                                                             cGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCAGTCACATCAGCCGCACGCAGCGGTTTTTACTATTGTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCg                                                                                                                                                                                    <  1:857441/251‑1 (MQ=255)                                                               attacgtcacgagcaagtgacctcttAACTTTCACAGCACAAACATAAATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAACATACGCAGAATTTTCCAGGGATTGGTCAAAAAGAGCTCATCCTGATTATGGTTTCACCATTCACATCAGCCCCACGCGGCGGTTTTTACTATTTTGCGCAAGGCTGGAAACTCGCCTCGCAACCTGGGATTCGGc                                                                                                                                                                                  >  2:943489/25‑251 (MQ=255)                                                                     ccaagtgcaagtgCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTCTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCg                                                                                                                                                                            <  2:735795/250‑1 (MQ=255)                                                                            caattacgcTATTATCTTTCACAGCACAAATATAGATGGAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAAGGTGCCGGTAAAATACGCAGAATTCTCCGGTGATTGTTCAAAAGTAGCTCATCCTGGTTATGGTGTCATCAGTCACAGCAGGCCGACGGGGCGGTTTTTACTATTGTGCGCAAGGCTGGAAGCTCGTCTCGCTGCCTGGGAGTCGGCGTATCGTTAtttt                                                                                                                                                                     <  2:612465/243‑1 (MQ=255)                                                                                        tAACTTTCACAGCACAAAGATCGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCACTTCCCCGGTTAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGt                                                                                                                                                         >  2:878354/1‑251 (MQ=255)                                                                                        tAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTCACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGt                                                                                                                                                         >  1:140202/1‑251 (MQ=255)                                                                                          aCTTTAACAGCACAAAGTGAGATGAAATCGTGCTTTTTGCTGTGTTTTCGAACATATCCTAACTGTCCATTGCGTAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATGTTGCGCAATGCTGGAAGCTCGTCTCGCAACCTGTGATGCGGCGTTGCGTTATTTTACCGCTGCTGGTCa                                                                                                                                                       <  2:650880/251‑1 (MQ=255)                                                                                              tCACAGCACAAAAATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAata                                                                                                                                                    <  2:492681/250‑1 (MQ=255)                                                                                                         agatagatGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATgg                                                                                                                                        >  2:730103/1‑251 (MQ=255)                                                                                                          gatagatGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATggg                                                                                                                                       >  2:20869/1‑251 (MQ=255)                                                                                                                 gAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGgcgc                                                                                                                                 >  2:531189/1‑250 (MQ=255)                                                                                                                    aTCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGAAATTACCCGGTAAAATACGCAGACTTTTCCTGGGATTGGCCAAAAGGAGCTCATCCTGATTATGGATTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTc                                                                                                                             >  1:431900/1‑251 (MQ=255)                                                                                                                                     tttttCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTacacc                                                                                                            >  2:534288/1‑250 (MQ=255)                                                                                                                                       tttCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCAGTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATAGTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAg                                                                                                           <  1:878354/250‑1 (MQ=255)                                                                                                                                            aaCATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGAt                                                                                                     <  1:534288/251‑1 (MQ=255)                                                                                                                                               atatCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTc                                                                                                  <  1:20869/251‑1 (MQ=255)                                                                                                                                                 atCCTAACTGTCCATTGCGCAATTATCCGGTAACATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTg                                                                                                >  2:778769/1‑251 (MQ=255)                                                                                                                                                     tAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGCTTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCATTATTTTGTTGATGGGGGGCGCATTCTGGTAGCTCTTTACACAGCTCGATGTCTGGagc                                                                                            >  1:903481/1‑249 (MQ=255)                                                                                                                                                      aaCTGTCCATTGCGCAAGTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGTGGCGCATGCTGGTGGCTCTTTACACAGCTCGATGTCTGGAt                                                                                             <  2:579184/249‑1 (MQ=255)                                                                                                                                                                            cGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTcc                                                                     >  1:80853/1‑251 (MQ=255)                                                                                                                                                                             ggTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCg                                                                    >  1:700340/1‑251 (MQ=255)                                                                                                                                                                                          gAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGCTGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGCCTGGCTACAAt                                                       >  2:530450/1‑251 (MQ=255)                                                                                                                                                                                               ttCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCt                                                   <  1:904473/250‑1 (MQ=255)                                                                                                                                                                                                      gATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATc                                           <  2:80853/251‑1 (MQ=255)                                                                                                                                                                                                              aaaaGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCAAACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTgg                                   >  1:903486/1‑251 (MQ=255)                                                                                                                                                                                                                 aGGAGCTCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTGTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTGACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCct                                <  2:903486/251‑1 (MQ=255)                                                                                                                                                                                                                       tCATCCTGATTATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGa                                                                                              >  1:125120/1‑183 (MQ=255)                                                                                                                                                                                                                       tCATCCTGATTATGGTTTCAGCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGGCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAAGATTATGTTGAGGGGGGGCGCATGCTGGTGGCTCTTTACGCAGCTCGATGTCTGGa                                                                                              <  2:125120/183‑1 (MQ=255)                                                                                                                                                                                                                                 tATGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTAGCTGTTGTGGCCTCTGGCTGTCATctct                 <  1:730103/250‑1 (MQ=255)                                                                                                                                                                                                                                  aTGGTTTCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCtgt               <  1:531189/251‑1 (MQ=255)                                                                                                                                                                                                                                       ttCATCATTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATctct                 <  1:773792/244‑1 (MQ=255)                                                                                                                                                                                                                                       ttCATCATTCACATCAGCCCCACGAAGCGGTTTTTACTATTTTGCGCAAGGCTGGACGCTCGTCTCGCAACCTGGGCTTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGCCTGGCTACAATGGCTGCGGTATCTGTTGTGGCCTCTGGCTGTCctctct                 >  2:773792/1‑244 (MQ=255)                                                                                                                                                                                                                                         catcatTCACATCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCTGTGCTGTTa        <  2:700340/251‑1 (MQ=255)                                                                                                                                                                                                                                                cacaTCAGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCTGTGCTGTTAGTGttt  <  1:187619/250‑1 (MQ=255)                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                    TTCTCTGTTGTTCTAACACCTTGCCACCACGGCAAACATTTACTCACTAAGAGTATTTGCCGATTACCTCAAGTGCAAGTGCACTATTAACTTTCACAGCACAAAGATAGATGAAATCGTGCTTTTTGCTGTTTTTTCGAACATATCCTAACTGTCCATTGCGCAATTACCCGGTAAAATACGCAGAATTTTCCTGGGATTGGTCAAAAGGAGCTCATCCTGATTATGGTTTCATCATTCACATCTGCCCCACGCAGCGGTTTTTACTATTTTGCGCAAGGCTGGAAGCTCGTCTCGCAACCTGGGATTCGGCGTTTCGTTATTTTACCGCTGCTGGTCAATATTTTGTTGATGGGGGGCGCATTCTGGTGGCTCTTTACACAGCTCGATGTCTGGATCCCGACTCTCATGAGTTACGTTCCGGACTGGCTACAATGGCTGAGTTATCTGTTGTGGCCTCTGGCTGTCATCTCTGTGCTGTTAGTGTTT  >  NC_000913/2531238‑2531726

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: ‑