Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,529,180	T→C	pseudogene (1279/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,180	0	T	C	93.8%	76.9 / ‑5.4	34	pseudogene (1279/2037 nt)	rhsE	pseudogene, Rhs family
Reads supporting (aligned to +/- strand):  ref base T (0/2);  new base C (13/17);  total (14/20)
Fisher's exact test for biased strand distribution p-value = 5.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

CCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCA‑CTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTT  >  NC_000913/1529020‑1529454                                                                                                                                                                 |                                                                                                                                                                                                                                                                                     cccGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATATTTATCA‑CTGCCGCCACCGGGGACCGCCGCTGGCCCTCATCAGCCAAGCCGGGAATCCGGCCTGGCGCGGGGCGTATGGTGCATGGGGACAACC‑GGTTCACGCGGAGCACCTGCATCACCTggcgca                                                                                                                                                     >  2:183379/1‑283 (MQ=11)   cGGCTGAACTGGTGCGGCGGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGcc                                                                                                                                                  <  1:433306/289‑1 (MQ=12)              gTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGCACCGTccg                                                                                                                                        >  1:469214/1‑286 (MQ=12)               tGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCA‑CTGCGACCACCGGGGCCTGCCCCTGGCGCTTATCAGCGAAGACGGCAAGACGGCGAGGGCCGGGGGGGAGGATGACTGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTAcacgctgc                                                                                                                                       >  2:491396/1‑281 (MQ=11)                     tgcggGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgca                                                                                                                                <  1:117011/289‑1 (MQ=12)                        tgGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgcagca                                                                                                                             <  1:542105/289‑1 (MQ=12)                        gtgccggccgggagaaATAATTCCGGGAGCACCGCGGGGTCAGTCGCAGCCGGCGGTGGCTTGCGCAGTGGCGGCAGACGGAGGCGCAACTTGCCCGACAGGTGGGGCCGGAATACCCACCGGCGCGCAACGCTCATTGTTGGCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGCCGGCCATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGgcagc                                                                                                                              <  2:24331/278‑1 (MQ=11)                               gCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACGCACCGGCGCGACAAGCTCATCTTTATCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGCGGCGCGTGGAGTATGGTGAATGGGG‑CAACCAGCTTAATGAGGAGCAACCGCCTTACCTGCACCAGCCGTACCCACTGCTCGGCCACCCGCatgaa                                                                                                                      >  2:267692/1‑288 (MQ=11)                                       aagaaTGCCGGGGAGACCGGGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGt                                                                                                               <  1:517078/288‑1 (MQ=255)                                         gaaATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTCTAAGCGAAGACGGCAATCCGCCGGGGCGCGGGGAGGACGATGACGGGGG‑AAACCGGCTTAATGAGGAGGAGCCGCATCACCTCGACCATCCGGACCGTCGCCCGGGGCAGAAGCaggaacgggggcggg                                                                                                            >  2:416961/1‑278 (MQ=11)                                           aaTCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACGGGCCAGGCAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCgggg                                                                                                          <  1:554984/289‑1 (MQ=255)                                            agccGGGGAGACCGGCTGGGCAGTGAAAGCCGGGCGTGGCTGGCGCAGTGCGGGCTGACGGTGGGGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGc                                                                                                         <  1:355302/287‑1 (MQ=11)                                            aaccGGGCAGACCGCGTGAGCACGGAAAGCCGGGCGGGGCTTGCGCGGTGCGGGCTGGCGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGc                                                                                                         <  1:267692/287‑1 (MQ=11)                                             accggcgaggccGCGTGGGCAGGGGAAGGGGGGCGTGGGTTGGGGAGGGGGGGCTGACGGTGGGGCAACGGGCAAGAAAGGGAGAGCCGGAATACGCACCGGCGCGAAAAGCTGATTTTTATCA‑ATGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGc                                                                                                         <  2:464156/278‑1 (MQ=11)                                                             ggcgagggAAGGCCGGGGGTGGGTTGCGGCGTGCGGGCGTAAGGGGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAATCGCAGCTTTATCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCAc                                                                                        <  2:469214/282‑1 (MQ=11)                                                                             gCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCCGCCGTACCGTCTGCCAGGGCAGCAACAGGATGAGGACTCGGGCCGGTAGGAACACCGTCACCCGccactcgacacgt                                                                        >  2:75068/1‑277 (MQ=11)                                                                                     tGCGCAGGGCGGGCTGACGGTGGGGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAgggg                                                                <  1:463751/289‑1 (MQ=255)                                                                                      gcgcAGTGCGGGCTGACGGGGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACCCACCGGCGCGAAAAGCTCATCTTTATCC‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAgggg                                                                <  1:456808/288‑1 (MQ=255)                                                                                        gcAGTGCGGGCTGCCGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCg                                                              <  1:491396/288‑1 (MQ=255)                                                                                            tGCGGGCTGACGGTGGGGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATc                                                                           <  1:555660/271‑1 (MQ=255)                                                                                            tGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGATGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGtata                                                         >  1:264564/1‑289 (MQ=255)                                                                                                 gCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGTATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCatg                                                                                                                        >  1:74375/1‑221 (MQ=255)                                                                                                   tGACGGTGGAGCAACGGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCA‑CTGCGACCACCGGGGACTGCCGCGGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCAGCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCAccc                                                   <  1:183379/288‑1 (MQ=255)                                                                                                      cagTGGATCAACTGGCCAGACAGGTGGAGCCGGAAGACACACCGGCGCGAAAAGCTCAACTTTAGCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGa                                               <  1:168336/287‑1 (MQ=255)                                                                                                      cGGTGGTGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCgg                                                >  1:614705/1‑288 (MQ=255)                                                                                                      cGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGACGAGGGCGCGGGGCGGTACTAGAACCGTGACCGCAACTCAGAGCCGTGGCAGGGGCgggcgctgagacacc                                                >  2:555660/1‑275 (MQ=11)                                                                                                      cGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCAACACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTTTATCACCCCgg                                                >  1:642594/1‑288 (MQ=255)                                                                                                       gggtgAGCCACTGGCCAGACGTGGTGAGCCGGAAGACACCCCGGCGCGAAACGTGCAGCTGGTGCA‑CGGAGACCACCGGGGAGTGCCGCTGGCGCTTATCAGCGAAGACGGCAATGCGGCGGGGCGCGGGGAGTATGATGAATGGGG‑CAAACAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGa                                               <  2:689080/284‑1 (MQ=14)                                                                                                               aaCTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATCTCACCCCGGACCAGAGTgg                                      >  1:689080/1‑289 (MQ=255)                                                                                                                   gtccagcaggggggAGCCGGAACAACCACCGGCG‑GAAAAAGCCATGCTTGTCACCTGCGACAACCGGGGACGGCCGCTGGGGCTTAACAGCGAAGACGGGAATAAGGCGTGGGGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCAACACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGtt                                   <  2:642594/276‑1 (MQ=14)                                                                                                                             gTGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCA‑CGGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGtgga                         <  1:518674/288‑1 (MQ=255)                                                                                                                              tGGAGCCGGAATACACACCGGCGCGAAAAGCTCATCTTTATCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGt                            <  1:462392/284‑1 (MQ=255)                                                                                                                                        aTACACACCGGCGCGAAAAGCTCATCTTTATCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTAt             >  2:360476/1‑289 (MQ=255)                                                                                                                                                   cgcgAAAAGCTCATCTTTATCA‑CTGCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGTTGAGGAGTCGGGGCTGGACTATAACCGTTCCCGGCCCTACGTACCGTTGCAGTGGCGGTATCTCACCCCGGCCCCGATTGGGATGGCAGGTGGTGGGAATATGTATCAGGATCCGtt  >  2:231998/1‑289 (MQ=14)                                                                                                                                                                 |                                                                                                                                                                                                                                                                                     CCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCA‑CTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTT  >  NC_000913/1529020‑1529454

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