Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,529,203	C→A	pseudogene (1302/2037 nt)	rhsE →	pseudogene, Rhs family

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,529,203	0	C	A	100.0%	92.6 / NA	41	pseudogene (1302/2037 nt)	rhsE	pseudogene, Rhs family
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base A (17/24);  total (17/24)

GGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGA  >  NC_000913/1529014‑1529490                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                ggtgTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTa                                                                                                                                                                                   >  1:598226/1‑300 (MQ=14)       tcgggCTGAACTGGTGGGGCTGTTGGACGGGCTGAAGGAAGAAATCCGGGCAGCCCGCGTGGGCAGTGAGAGCCGGGCGTGGCTTGCGCAGTGCGGCCTGACGGTTGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGGGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGGGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTg                                                                                                                                                                            <  1:214896/298‑1 (MQ=11)        ccGGCTGAACAGGTGCGGTTGCTGCACAGTCTGGAGGAAGAAATCCGCGCAGACCGCGTGAGCAGTGAAAGCCGGGAGTGGCTTGCGCAGCGCGGGCTGACGGTCGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTCACCACTGCGACCATCTGGGACTGCTGCCTGCGCTTTTCAGCGAATACGCTAATACGGCTTGGTGCGCGGAGTATGATGAAAGGCGCAACCAGCCTCATTAGGCGAACACGCATCACCTCCACCAGCCGTACAGTCTg                                                                                                                                                                            >  2:213183/1‑300 (MQ=11)          ggCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGTAAAGCTCATCTTTATCAATGCGACCACCGGGGACGGCCGACGGCGCTTGTCCGCGGAGACCGAAATAAGGCGCGGCGCGGGGAGTATGATGAATGGGGCGACCATCATAAAGAGGAGGACTGGCGTCACATCCACCAACCGTTACGTccgagg                                                                                                                                                                         >  2:214894/1‑296 (MQ=11)          ggCTGAACGGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAGAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGGGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGcc                                                                                                                                                                          <  2:346850/300‑1 (MQ=9)              gAAGTGGTGCGGCGGCGGAGCAGGCTAGAGGAAGAAAGCCGGGCAGCCGGCGTGGGCAGTGAAAGGCGGGCGTGGGTTACGGAGGACGGACGGACGATGGAGCAAATGGACAGAGAGGTGGAACCGGGATACAAGCCGGGGGGAAAAGCTCATCTTTATCACGGCGACCACCGGGGACTGCGGCTGGCGCTGATCAGCGGAGACGGCAATACGGCGTGGCGCGGGGAGTAGGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAggg                                                                                                                                                                      <  2:385722/300‑1 (MQ=11)               aaCTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgc                                                                                                                                                                     <  1:507847/300‑1 (MQ=14)                               ggACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGTCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATTATTGCGACCACTGGGGATTGCCGCCGGCGATTATCAGAGAAGACGGCAATACGGCGGGGCGCGGGGAGGACGATGAAGGGGGCCACCAGCTTAATGAGGAGAACCAGCAACACCTGAACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATgaggag                                                                                                                                                    >  2:212333/1‑301 (MQ=11)                                gcaaGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAACGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATgaggag                                                                                                                                                    <  2:323183/297‑1 (MQ=12)                                     gCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCgggg                                                                                                                                              >  1:432742/1‑301 (MQ=18)                                     gCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCggg                                                                                                                                               <  1:212335/300‑1 (MQ=18)                                      cTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTTAAAGCCGCGCGTGGCTTGCGCAGTGCGGGCTGACGGTGAAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGTCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGc                                                                                                                                             <  1:570099/301‑1 (MQ=255)                                       tGGGGGAAGAAATCCGGCCAGACCGCGTGAGCAGTGAACGCCGTGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGCCAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGc                                                                                                                                             <  2:598220/300‑1 (MQ=255)                                       gggaggaAGAAATCCGGGCAGACAGAGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGGGCGGGCTGACGGTGGAGCAACGGGCCAGAAAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGc                                                                                                                                             <  2:43207/299‑1 (MQ=11)                                        ggaggaAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGTAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCgggg                                                                                                                                              >  2:128894/1‑298 (MQ=16)                                                aaaTCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACtata                                                                                                                                    <  1:110088/300‑1 (MQ=18)                                                         cAGACCGCGTTGGCAGTGAAAGCCGGGAGTGGCTTGGGCAGTGCGGGCTGACGGTGGAGAAACTGGCCAGACAGGTGGAGCCGGAATACACGCCGGCGCGAAAGGCTCAGCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATAAGGAGTCGGGGCTGTACTATAACCGTCAcc                                                                                                                           <  2:109119/300‑1 (MQ=11)                                                            ctcttccgatctgTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCa                                                                                                                       <  1:512650/289‑1 (MQ=255)                                                                 gTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTCCCTTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTTACCGGCACTAc                                                                                                                   >  1:115402/1‑300 (MQ=255)                                                                 gTGAGCAGGGAAAGCCGGGCGGGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTAc                                                                                                                   <  2:33673/300‑1 (MQ=255)                                                                   gAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGGACTATAACCGTCACCGGCACTACGaa                                                                                                                >  1:194590/1‑300 (MQ=255)                                                                   ctctttcccaaaacgacgcggttccgagctgTGCGGGAGGAGGGTGGAGCAACTGGCCAGGCAGGGGGGGCCGGAATACACACCGGCGCGAAAAGCTCATATTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGAt                                                                                                                <  2:570483/271‑1 (MQ=11)                                                                       aGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCGGACAGGTGGAGCCGGAATACACACCGGCTCGAAAAGCTCTTCTTTATCCCTGAGAACACCGGGGACTGCTGCTGGCGCTTACCAGCGAAGAGGGCAATACGGCGCGGCGCGGGGAGTACGATGTATGGGGCAAACAGCTTAATGGGGAGAATCCGCACATCCTGCACCAGTCGCACAGTCCTCCAGGGCAGCAGCATGACGAGGAGTAGGGGACTTATTATAACTGCCACCGGGAcatcgagacgg                                                                                                            >  2:96939/1‑291 (MQ=2)                                                                        gTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCCGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCAATatcggaagag                                                                                                           >  2:512645/1‑292 (MQ=255)                                                                          gAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACCGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCCCGAAAAGCCCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATTAGCGAAGACGGTATTACGGCCTGGCGCGGGGAGTATGACGAACGAGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGCACGGTCTGCCAGGCCAGCCTCTTGATGAGGAGTCGGGGGCGTACTGTCACCGTCACCGGTACTCCGATGCGCTg                                                                                                          >  2:21789/1‑300 (MQ=11)                                                                                    tacccgacgctcttccgatcttGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTTCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGta                                                                                               <  2:224654/280‑1 (MQ=11)                                                                                     gTGGCGTGCGCAGTGCGGGCGGACGGTTGAGCAACTGGCCAGACAGGTGGAGCCGGAAGACACACCGGCGCGAAAAGCTCATCTTTATCGCTGCGACCACCGGGGACTGCCGATGTCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGTAGTATGATGAATGGGTCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCACCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGta                                                                                               <  1:636271/300‑1 (MQ=11)                                                                                                 gTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGAtagatcggaagagcatactgctgaactccat                                                                                  >  1:570488/1‑271 (MQ=11)                                                                                                    cGGGGTGACGGTAGAGCAACGGGCCAGAAAGGGTGAGCAGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCTCCGGGGACTGCCGCGGGCGCGTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATAATGAATGGGGCAACCAGCTTAATGTGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGAccc                                                                                <  2:115400/300‑1 (MQ=14)                                                                                                         tGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTTCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGtaagatcggaagagcacacgt                                                                            >  1:224656/1‑280 (MQ=11)                                                                                                                  aGCAACAGACCAGTCATGTGTGGCCGGAAGACACACCGGCGCGAAAAGCTCATCTTTATCACGGCGACCACCGGGGACGGCCCCTGGCGCTTAACAGCTAAGACGGCAGTACGGCGTTGCGCGGGGAGTATGATGAATGGGGCAACCAGCGTAATGAGGAGAACCCGCAGCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGGCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTGGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGagag                                                                  <  2:307978/300‑1 (MQ=14)                                                                                                                    cAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAgg                                                                 <  1:128896/299‑1 (MQ=255)                                                                                                                         ggCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCGCTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGtggatg                                                           <  1:271280/300‑1 (MQ=255)                                                                                                                            aaGAAAGGTGGGGCCGGAATCCACACCGGCGCGAAAAGCTCATCTGTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGCAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAat                                                       <  2:194588/300‑1 (MQ=255)                                                                                                                                                      cGGCGCGAAAAGCTCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCAt                             >  1:228678/1‑301 (MQ=255)                                                                                                                                                                 gctCATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAATTGTCACCGGCACTACGACCCGTTGCAGGGGCGGTATATCACCTCGGACCCGATTGGGTTGCGGGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTgttcg                  >  2:53801/3‑297 (MQ=17)                                                                                                                                                                    cATCTTTATCACTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGAc                >  1:23285/1‑300 (MQ=255)                                                                                                                                                                         ttATCACTGCGACCACCGGGGACGGCCGCTGGCGCTTATCAGCGAAGACGGCAAGACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGAACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGGATATGTATCAGTATCCGTTGAAACCCATACAAGTGATAGACCCAATg          <  1:21789/301‑1 (MQ=255)                                                                                                                                                                            ccgcTGCGACCACCGTGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAAAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGGTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATgggg       <  2:506780/298‑1 (MQ=255)                                                                                                                                                                               cTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCCCTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGCTTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGACTCCCATCCAAGTGATAGACGCAATggggaga    >  1:506785/1‑298 (MQ=255)                                                                                                                                                                                  cGACCCCCGGGGACTGGCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGa  <  1:617974/300‑1 (MQ=255)                                                                                                                                                                                              |                                                                                                                                                                                                                                                                                                GGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGA  >  NC_000913/1529014‑1529490

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