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	97.4%	100.2 / NA	39	pseudogene (1302/2037 nt)	rhsE	pseudogene, Rhs family
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (16/22);  total (17/22)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.67e-01

CCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCA‑CTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAAT  >  NC_000913/1529020‑1529481                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                        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)                                                                                                                                                                     atcA‑CTGCGACCCCCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGAt           <  1:181753/287‑1 (MQ=255)                                                                                                                                                                       cA‑CTGCGACCACCGGGGACTGCCGCTGGGGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAg         <  1:75068/289‑1 (MQ=255)                                                                                                                                                                            gCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTCCGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGGGGGGGGTGGAATATGGATCAGGCTCCGTTGAATCCCAGCCAAGGGATAGAccc     >  2:303008/1‑289 (MQ=17)                                                                                                                                                                            gCGACCACCGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGACTGGGTTGAGAGGTGGATGGAATCCGTATCAGTCTCCGTTGAATCCCATACAAGTGATAGAccc     >  1:462371/1‑289 (MQ=255)                                                                                                                                                                               accaccGGGGACTGCCGCTGGCGCTTATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGAACCGAGTGGGTTGAGAGGTGGATGGAATATGTAACAGTATCAGGTGAATCCCATACAAGTGATAGACCCAAt  >  1:357993/1‑289 (MQ=255)                                                                                                                                                                                         |                                                                                                                                                                                                                                                                                        CCCGGCTGAACTGGTGCGGCTGCTGGACAGGCTGGAGGAAGAAATCCGGGCAGACCGCGTGAGCAGTGAAAGCCGGGCGTGGCTTGCGCAGTGCGGGCTGACGGTGGAGCAACTGGCCAGACAGGTGGAGCCGGAATACACACCGGCGCGAAAAGTTCATTTTTACCA‑CTGCGACCACCGGGGCCTGCCGCTGGCGCTCATCAGCGAAGACGGCAATACGGCGTGGCGCGGGGAGTATGATGAATGGGG‑CAACCAGCTTAATGAGGAGAACCCGCATCACCTGCACCAGCCGTACCGTCTGCCAGGGCAGCAGCATGATGAGGAGTCGGGGCTGTACTATAACCGTCACCGGCACTACGATCCGTTGCAGGGGCGGTATATCACCCCGGACCCGATTGGGTTGAGAGGTGGATGGAATATGTATCAGTATCCGTTGAATCCCATACAAGTGATAGACCCAAT  >  NC_000913/1529020‑1529481

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