Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,400,978	C→T	G266S (GGT→AGT)	abgT ←	p‑aminobenzoyl glutamate:H(+) symporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,400,978	0	C	T	100.0%	155.5 / NA	44	G266S (GGT→AGT)	abgT	p‑aminobenzoyl glutamate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (23/21);  total (23/21)

CTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACCTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATTTGCGGATTGAACGCAGCTGCCGCTTCCGTGCTGATCCCCGACAGCAACACGTCGGTTGTGAC  >  NC_000913/1400750‑1401221                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    ctTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAAc                                                                                                                                                                                                                                 <  2:112281/249‑1 (MQ=255)          gTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCATAAGCAAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCg                                                                                                                                                                                                                      >  2:421890/1‑251 (MQ=255)            tCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCATAAGCCAGCGCGACAACAAAGACAAAAAGACTGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGCCACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGt                                                                                                                                                                                                                    >  1:432440/1‑251 (MQ=255)                              aaaTCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCATAAGCCAGCGAGGCAACAAAGAAAAAAAGCATGATCAGTGGCACGCTCCCTTTAATAAAGGGTGATGGCATCACGGTGTGATTACTCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCCATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTtctc                                                                                                                                                                                                  >  1:284693/1‑251 (MQ=255)                                     ccTGACGTCGAATTGTGCGGGTAGCGATGCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTg                                                                                                                                                                                           <  2:233150/251‑1 (MQ=255)                                             cGAATTGTGCGGGTAGCGATGCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTg                                                                                                                                                                                   <  1:344274/251‑1 (MQ=255)                                                                       aGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCgatg                                                                                                                                                          >  1:32762/1‑250 (MQ=255)                                                                         ccAGCGAGACAACAAAGAAAAAAAGAAGGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGAtt                                                                                                                                                       <  1:188423/251‑1 (MQ=255)                                                                          cAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGAttt                                                                                                                                                      >  1:307203/1‑251 (MQ=255)                                                                                    acaaAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTAt                                                                                                                                            >  2:308560/1‑251 (MQ=255)                                                                                    acaaAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTAt                                                                                                                                            <  1:239170/251‑1 (MQ=255)                                                                                      aaaGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATc                                                                                                                                           >  2:291315/1‑250 (MQ=255)                                                                                         gAAAAAAAGAATGATCAGTGGCACGATACCTTTAATACAGGGTGATGGCATAACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGgc                                                                                                                                       >  2:407486/1‑251 (MQ=255)                                                                                           aaaaaaGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGgccg                                                                                                                                     <  1:371526/251‑1 (MQ=255)                                                                                            aaaaaGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGgccgc                                                                                                                                    >  2:250037/1‑251 (MQ=255)                                                                                             aaaaGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTTACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGgccgcc                                                                                                                                   >  2:240699/1‑251 (MQ=255)                                                                                                       tCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTc                                                                                                                         >  1:381030/1‑251 (MQ=255)                                                                                                                cGATTCCTTTAATAAAGGGTGATGGCAACACGGTGTGATTAATCGGATCGCGCAATATCCCGTTGTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTTTGCGTAAACCAAAACGCTGATTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTTCCATGGACCTAACCGGGGCTCGATGATTTTGGCGGTTATCAGGCCGCCAACAATCGTCAGTACGACt                                                                                                                <  2:386662/251‑1 (MQ=255)                                                                                                                  aTACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTAc                                                                                                              >  1:198316/1‑251 (MQ=255)                                                                                                                          aataaAGGGGTAGGGCATCACGGTGTGATTAATCGGGTCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGAAATCGCAGCAATAAAAAGTAGCGATACGCCACTTGCTATGCGTAAACCAAAACGCTGACGTTCGGTCAATGTCTGCGGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGGTCGATGATTGTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTgg                                                                                                      <  2:167252/251‑1 (MQ=255)                                                                                                                           ataaAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAAGATCCCGTTGTGCGGGATCTCCATCAGCGCAATCGCAGCAATAAAAATGAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCGGCAGTTTCTCATCGCTGTTTCCCTGCCATTGCCCTAACCGTGGCTCGATGATGTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGTCTACGGAGCTGGc                                                                                                     <  2:404814/251‑1 (MQ=255)                                                                                                                             aaaGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCa                                                                                                   <  1:379233/251‑1 (MQ=255)                                                                                                                                 ggTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATaa                                                                                                >  2:151243/1‑250 (MQ=255)                                                                                                                                 ggTGATGGCATCACGGTGTGATGAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATTCGTAAACCAAAACGCTGACTTTCGGTCATTGTCGGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTATCCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATaaa                                                                                               <  2:252959/251‑1 (MQ=255)                                                                                                                                          aTCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGtt                                                                                      >  1:107674/1‑251 (MQ=255)                                                                                                                                                  gtgATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAAtt                                                                               >  2:84843/1‑250 (MQ=255)                                                                                                                                                   tgATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGgccg                                                                                                                                     >  1:178908/1‑195 (MQ=255)                                                                                                                                                   tgATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTac                                                                             >  1:345783/1‑251 (MQ=255)                                                                                                                                                   tgATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCAGCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGgccg                                                                                                                                     <  2:178908/195‑1 (MQ=255)                                                                                                                                                       tAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTg                                                                         >  1:329368/1‑251 (MQ=255)                                                                                                                                                         aTCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGAc                                                                       <  1:407486/251‑1 (MQ=255)                                                                                                                                                                cgcgCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGttt                                                                                                                                                                                        <  2:119704/131‑1 (MQ=255)                                                                                                                                                                cgcgCAATATCCCGTTTTGCGGGATCACCATAAGCGCAATCGCAGCACTAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGttt                                                                                                                                                                                        >  1:119704/1‑131 (MQ=255)                                                                                                                                                                    cAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAACTACCAGTTATCAATTACACTGACGTGCATTTGCg                                                            >  2:48556/1‑251 (MQ=255)                                                                                                                                                                        atCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATc                                                                                                                                                                                               <  1:207536/116‑1 (MQ=255)                                                                                                                                                                        atCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATc                                                                                                                                                                                               >  2:207536/1‑116 (MQ=255)                                                                                                                                                                            cGTTGTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCAGTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAAGTACACTGACGTGCATTTGCGGATTGAAc                                                    <  2:271776/251‑1 (MQ=255)                                                                                                                                                                                ttGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGTCCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGGTAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATTTGCGGATTGAACGCa                                                 <  1:277463/250‑1 (MQ=255)                                                                                                                                                                                                        cGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAGCGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATTTGCGGATTGAACGCAGCTGCCGCTTCCGTGCTGATcccc                         <  2:141932/250‑1 (MQ=255)                                                                                                                                                                                                        cGCAGCAATAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATTTGCGGATTGAACGCAGCTGCCGCTTCCGTGCTGATCCCCg                        <  2:329368/251‑1 (MQ=255)                                                                                                                                                                                                               aTAAAAAGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATTTGCGGATTGAACGCAGCTGCCGCTTCCGTGCTGATCCCCGACAGCaa                 <  1:240699/251‑1 (MQ=255)                                                                                                                                                                                                                   aaaGTAGCGATACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATTTGCGGATTGAACGCAGCTGCCGCTTCCGTGCTGATCCCCGACAGCAACACg             >  1:288786/1‑251 (MQ=255)                                                                                                                                                                                                                            aTACGACACTTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATTTGCGGATTGAACGCAGCTGCCGCTTCCGTGCTGATCCCCGACAGCAACACGTCGGTtgtg    >  1:389051/1‑251 (MQ=255)                                                                                                                                                                                                                              aCGACACTTGCTATGCGTAAACCAAAACGCTGACGTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGGGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATTTGCGGATTGAACGCAGCTGCCGCTTCCGTGCTGATCCCCGTCAGCAACACGTCGGTTGTGAc  <  2:432440/251‑1 (MQ=255)                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                    CTTTCATCGGTTCAATCATTAAATGCGGTAAATCCGCCTGACGTCGAATTGTGCGGGTAGCGATGCCATAAGCCAGCGAGACAACAAAGAAAAAAAGAATGATCAGTGGCACGATACCTTTAATAAAGGGTGATGGCATCACGGTGTGATTAATCGGATCGCGCAATATCCCGTTTTGCGGGATCACCATCAGCGCAATCGCAGCAATAAAAAGTAGCGATACGACACCTGCTATGCGTAAACCAAAACGCTGACTTTCGGTCAATGTCTGCAGTTTCTCATCGCTGTTTCCCTGCCATTGACCTAACCGTGGCTCGATGATTTTGTCGGTTATCAGGCCGCCAACAATCGTCAGTACGACTACGGAGCTGGCCATAAAATACCAGTTATCAATTACACTGACGTGCATTTGCGGATTGAACGCAGCTGCCGCTTCCGTGCTGATCCCCGACAGCAACACGTCGGTTGTGAC  >  NC_000913/1400750‑1401221

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