Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,529,218	C→T	pseudogene (1317/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,218	0	C	T	100.0%	150.6 / NA	50	pseudogene (1317/2037 nt)	rhsE	pseudogene, Rhs family
Reads supporting (aligned to +/- strand):  ref base C (0/0);  new base T (26/24);  total (26/24)

GTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAA‑GGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATAT  >  NC_000913/1528922‑1529505                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                 gTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAA‑GGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACAGGGGCCTGCCGCTGGCGCttatc                                                                                                                                                                                                                                                                                               <  2:551358/300‑1 (MQ=17)   cGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAA‑GGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAg                                                                                                                                                                                                                                                                                             <  2:440341/300‑1 (MQ=17)           agaACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAA‑GGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAg                                                                                                                                                                                                                                                                                        <  2:23377/297‑1 (MQ=17)                 gCGAGCGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAA‑GGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATa                                                                                                                                                                                                                                                                               <  2:506020/300‑1 (MQ=17)                             aaacccgagggcgCAGCCT‑GCCGAGCGGCTCAAGGGAGAAGGGGAGGGGGAACGACCGAGGGGGGGGGTTCCTGGCGGAACGGGTGCGGCTGATAGAGAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGAAGTAAAAGGCGGGGGTGGCGTGCGCAGTGCGGGGTGACGGTGGAGAAACTGGCCAGACGGGTGGAGCCGGAAAACACGCCGGCGCGAAAAGTTCATTTTTACCCCTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCgg                                                                                                                                                                                                                                                                  <  1:205056/292‑1 (MQ=2)                                cgaacgggcgatgCCTGGCAGAGCAGATCCAGCAGGAA‑GGGAGGGAGAGCGGCCAAGGCGTGGTGTTCCCGGCTGCACGGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCgggg                                                                                                                                                                                                                                                                <  1:278897/288‑1 (MQ=11)                                      ggcgCAGCCTGGCAGAGACGCTCCAGCAGGAA‑GGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTat‑gg                                                                                                                                                                                                                                                        >  2:511327/1‑300 (MQ=25)                                        cgcgGCCTGGCAGGGACGCTCCAGAAGGAA‑GGGAGTGAGAACGGCCACGGCGGGGTGTTCCCGGCGGAACAGGGGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGAAGTGAAGGCCTGTCGTGGCTTGCGCAGTGCGGGCTGAGGGTGAAGCAACTGGCAAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCACTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTat‑gatg                                                                                                                                                                                                                                                      <  2:268344/301‑1 (MQ=11)                                         gcAGCCTGGCAGAGACGCTCCAGCAGGAA‑GGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTat‑gatg                                                                                                                                                                                                                                                      <  2:378444/300‑1 (MQ=25)                                                                aagGAA‑AGGAGTGAGAACGGCCACGGCGTGGTGTGCCCGGCTGAACTGGGGCGGCTGCTAGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATgagg                                                                                                                                                                                                                              <  1:27856/300‑1 (MQ=255)                                                                     a‑GGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACCCCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCCATACGGCGTGGCCGGGGGGGTAG‑GGTGAAGGGGGCAACCAGCGTAATgaggagga                                                                                                                                                                                                                          >  1:264549/1‑298 (MQ=11)                                                                                 aCGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATTTTTAACACTGCGACCACCGGGGCCTGCCGCCGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GAGGAATGGGGCCACCACCGTAGTGAGGCGCACACGCAGCAcc                                                                                                                                                                                                               >  1:448020/1‑300 (MQ=11)                                                                                  cGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCAc                                                                                                                                                                                                                >  2:507513/1‑298 (MQ=37)                                                                                      cttccgatcttgtTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCAcc                                                                                                                                                                                                         <  2:397736/291‑1 (MQ=255)                                                                                                tgtTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCCGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCCACCAGCTTAATGAGGAGAACCCGCATCACCTGCCCCAgaccggcag                                                                                                                                                                                               >  1:397268/1‑293 (MQ=255)                                                                                                 gtTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAACGTTCATTTTTACCACTGCGCCCACCGGGGCCTGCCGCTGGCGCTTATCAGCGTACGCGTCAATACCGCGTGGCTCGGGGAGTATCGAT‑CATGTGTCCACCAGCGTTAGGACGACAACCCGCCACACTCGCCCCAACCGTACCg                                                                                                                                                                                               >  1:427028/1‑300 (MQ=11)                                                                                                 gtTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGt                                                                                                                                                                                              >  2:417684/1‑301 (MQ=255)                                                                                                      cccgacgctcttccgatctgctgGACAGGCGGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGGGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGcc                                                                                                                                                                                         <  1:112411/284‑1 (MQ=255)                                                                                                            aaCGGGGGCGGCTGATGGACAGGCGGGAGGAAGAAATCCGGGCAGACCGCGTGGGCAGTGAAAGCCGGGCGTGGCGTGGGTAGGGCGGGCTGACGGTGGAGCAACGGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgc                                                                                                                                                                                    <  1:417189/300‑1 (MQ=14)                                                                                                             aCTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgcag                                                                                                                                                                                  <  2:264862/301‑1 (MQ=255)                                                                                                                         gctgGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAgatcggaagagcacacgt                                                                                                                                                                      >  2:112551/1‑284 (MQ=255)                                                                                                                         gctgGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATgag                                                                                                                                                                      >  2:17017/1‑301 (MQ=255)                                                                                                                          ctgGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTTACGGTGGAGCAACTGGCCCGACAGGTGGAGCCGGAATACACACCGGCGCGACAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCTAAGACGGCAATACGGCGCGCCCCGGGTAGTCT‑GATTAATTGGGCAGCCAGCCTAATGAGCAGAACCCGCTACAACGTCACCAGCCGTCCCGTCTGACACGGCAGCAGCATGCTgagg                                                                                                                                                                     >  1:121165/1‑301 (MQ=14)                                                                                                                                 ggCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTACTGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCCGCCAGGGCCGCAGCATGAATCGGGGTCggg                                                                                                                                                              >  1:418686/1‑301 (MQ=255)                                                                                                                                  gCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCgggg                                                                                                                                                             >  2:515398/1‑301 (MQ=255)                                                                                                                                   cTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATCATGAGGAGTCggggg                                                                                                                                                            >  1:88167/1‑300 (MQ=255)                                                                                                                                        ggaAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTa                                                                                                                                                        >  2:196053/1‑300 (MQ=255)                                                                                                                                                ggcgggaggacgaggggaggagTGAAAGCCGTGACAGGCTTGCGCAGGGCGGGCGGACGGTGGGGAAACTGGCCAGAAAGGTGGAGCCGGAAGACGCAGCGGCGCGAAAAGTTCAGTATTAGCACTGCGACCACCGGGGCCTGCCGCGGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAATCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAAcc                                                                                                                                                <  1:106718/284‑1 (MQ=14)                                                                                                                                                 ccGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGt                                                                                                                                              >  2:493723/1‑301 (MQ=255)                                                                                                                                                  cGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGt                                                                                                                                              >  2:117627/1‑300 (MQ=255)                                                                                                                                                          ccGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCa                                                                                                                                      >  2:164272/1‑300 (MQ=255)                                                                                                                                                                     gTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGAATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGAACTACGAGCCGtg                                                                                                                          >  1:96205/1‑300 (MQ=255)                                                                                                                                                                      ggAAAGCCGGGGGGGGCTGGCGAAGTGCGGGCTGACGGGGGAGCAACTGGCCAGACAGGTGGGGCCGGAATACACACCGGCGAGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATc                                                                                                                              <  1:510693/295‑1 (MQ=255)                                                                                                                                                                                                 cGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGAccc                                                                                               <  2:88274/300‑1 (MQ=255)                                                                                                                                                                                                   ggCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGAt                                                                                            >  2:159600/1‑301 (MQ=255)                                                                                                                                                                                                         cGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGtt                                                                                      >  2:398438/1‑301 (MQ=255)                                                                                                                                                                                                              aaGCAACTGACCAGACAGGTGAGACCGGAATACACACCGGCGCGAAAAGTTCATTGTTACCACTGCGACAACCGGGGCCGGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTgaga                                                                                  <  1:159403/299‑1 (MQ=255)                                                                                                                                                                                                                    cTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCCCCCCGGACCCGATTGGGTTGGGGGGtgggt                                                                           >  1:239155/1‑299 (MQ=255)                                                                                                                                                                                                                                  ggAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGTGAAGACGGCAATACGGCGTGGCGCGGGGCGTAT‑GATGAAGGGGGCAACCAGCTTAATGAGGAGAACCAGCATCACCTGCCCCAGCCGTACCGTCTGCCACGGCCGCAGCTGCATTAGGACGCGGGGCTGTCCTATAACCTTCACCGGCACTCCGATCCGTTTCAGGGGCGGTATATCACCCCCGCCCCGATTTGGTTTAGGGGTGGAGGGAatagttagtcg                                                             >  1:473092/1‑293 (MQ=255)                                                                                                                                                                                                                                    aacaggAATAGACAGCGGCGCGAAAAGTTCAGTTTTACCACTGCGACCACCGGGGCCTGCGGGTGGCGGTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGGAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGGATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCTGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTa                                                           <  1:16994/297‑1 (MQ=17)                                                                                                                                                                                                                                     gCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGAACCGTTGCAGGGGCGGTTTTTCACCCCGGACCCGACTGGGTTGAGAGGTGGATGGCATATGTATCAGTAt                                                          >  2:516166/1‑301 (MQ=255)                                                                                                                                                                                                                                      ccGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATc                                                         <  2:239451/301‑1 (MQ=255)                                                                                                                                                                                                                                          aataccacccGGCGCGAAAAGTCCATTTTTACCACGGGGACCACCGGGGCCGGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGtt                                                     <  1:493112/293‑1 (MQ=255)                                                                                                                                                                                                                                                cacCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCCTGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGCTTGGGTTGAGAGGTGGATGGAATCTGTATCATTATCTGTTGAATc                                                >  1:245580/1‑300 (MQ=255)                                                                                                                                                                                                                                                      cgcgAAAAGTTCATTTTTACCACGGCGACCACCGGGGCCTGCCGCTGGCGCTTAGCAGCGAGGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTAGATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATAc                                          <  2:473686/300‑1 (MQ=255)                                                                                                                                                                                                                                                           aaaGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGa                                    >  2:170593/1‑301 (MQ=255)                                                                                                                                                                                                                                                                        cccACTGCGACCACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATggg                       <  2:96332/300‑1 (MQ=255)                                                                                                                                                                                                                                                                          cccgacggtcttccgatcttGCCGCGGGCGCTTATCAGCGAAAACGGCAAGACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGt                     <  1:504940/282‑1 (MQ=255)                                                                                                                                                                                                                                                                               cGACAACCGGGGCCTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGa                 <  2:448567/300‑1 (MQ=255)                                                                                                                                                                                                                                                                               cGAACACCGGGGCCTGCCGCGGGAGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGa                 <  1:506879/300‑1 (MQ=255)                                                                                                                                                                                                                                                                                             tGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGtagatcggaagagcacacgt  >  2:505574/1‑282 (MQ=255)                                                                                                                                                                                                                                                                                                          |                                                                                                                                                                                                                                                                                                 GTCGAGACAGAGAACGGCGAGCGGGAAAAAGCGCAGCGGCGCAGCCTGGCAGAGACGCTCCAGCAGGAA‑GGGAGTGAGAACGGCCACGGCGTGGTGTTCCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCACTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTAT‑GATGAATGGGGCAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAATGGGGTTAGATGCGATTGAGAATAT  >  NC_000913/1528922‑1529505

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