Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,947,481	1	.	A	36.7%	94.9 / 117.1	60	intergenic (+18/‑16)	metZ/metW	tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base . (3/35);  new base A (2/20);  total (5/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.67e-01

CGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAAC  >  NC_000913/2947330‑2947571                                                                                                                                                         |                                                                                           gggggTATCAAACAGGTTCCATGAAAACTAGTGCGGAGTATTCTGCGCATCCGCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGt                                                                                             <  5:53726/150‑1 (MQ=14) gggggTAAACAACAGGGTCCATGAAAACGCGAGCGGAGTATTGTGCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGt                                                                                             <  3:570754/150‑1 (MQ=17) cgtgttAAAACAAAGGTTCCATGAAAACCCCAGAGGTGTATAGTGCGCCTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTGGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGt                                                                                             <  5:1563962/146‑1 (MQ=14)  ggggTAAAAAAAAGGGTGAATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑a                                                                                           <  3:8128/151‑1 (MQ=255)   ggttAAAATTAAGGTTGAATGAAATCTCGAGTGGAGTATAGTTCGATTCCTGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aa                                                                                          <  1:880808/148‑1 (MQ=11)   ggttAAAAAAAAGGTGCAATGAAAAATCCTGCGTAGTATTTTGCGCTTCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aa                                                                                          <  4:613612/148‑1 (MQ=11)    ggTAAACAAAAGGTGGCATGAAAACCCGAGCGGAGTATAGTCCGCGTCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aaa                                                                                         <  8:995072/151‑1 (MQ=17)    ggTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑aa                                                                                          <  3:856653/151‑1 (MQ=255)     gTAAAAAAATGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑aaa                                                                                         <  6:1228663/151‑1 (MQ=255)     gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑aaa                                                                                         <  2:634538/151‑1 (MQ=255)     gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑aaa                                                                                         <  2:956542/151‑1 (MQ=255)     gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑aaa                                                                                         <  1:1430461/151‑1 (MQ=255)     gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑aaa                                                                                         <  5:321800/151‑1 (MQ=255)       aaaaaaaaGGTTGCATGAAAACGCTAGCGGAGTATAGTGGGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aaaag                                                                                       <  1:1318018/150‑3 (MQ=17)       aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAg                                                                                        <  1:232605/150‑1 (MQ=255)       aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGc                                                                                       <  6:1261308/151‑1 (MQ=255)        aaaataaGGTTGCATGAAAACTCTATCGTAGTATTGTTCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aaaagca                                                                                     <  6:1036540/151‑5 (MQ=14)          aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGt                                                                                    <  4:1224884/151‑1 (MQ=255)          aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGt                                                                                    <  5:98857/151‑1 (MQ=255)           cggatcttGAATGAAAACTCTAGCGTAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGt                                                                                    <  1:1050308/144‑1 (MQ=14)           aaacGGGTGCATGAAAACGCCAGCGGAGTATATTGGCCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTa                                                                                   <  3:1409297/151‑1 (MQ=14)           aaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAACTAAAAGCAGt                                                                                    >  1:1226184/1‑150 (MQ=17)            aacggTTGCATGAAAAAGCGATCGGAGTATAGTGGGCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTAc                                                                                  <  6:307127/148‑1 (MQ=14)            aaaGGTTGCATGAAAACTCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTAc                                                                                  <  3:1223976/151‑1 (MQ=17)            aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTAc                                                                                  <  4:535577/151‑1 (MQ=255)            aaaGGTAGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTAc                                                                                  <  1:570886/151‑1 (MQ=255)             acggTTGAATGAATACTCTAGCGTTGTATATTGCGATTCCTGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACg                                                                                 <  4:1205717/149‑1 (MQ=11)             aaGGTTGCATGAAATATCTTGCGGAGTATATTGCGCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACg                                                                                 <  4:377885/151‑1 (MQ=11)             aaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACg                                                                                 <  8:239452/151‑1 (MQ=255)             aaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACg                                                                                 <  1:1240649/151‑1 (MQ=255)              aGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGTGCTCATAACCCGAAGGTCGTCGGTTCAAATCCTGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACgg                                                                                <  3:762905/151‑1 (MQ=255)              aGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACgg                                                                                <  1:817084/151‑1 (MQ=255)               ggTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGTGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGt                                                                               <  4:275028/151‑1 (MQ=255)               ggTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGt                                                                               <  5:91189/151‑1 (MQ=255)                 gtGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGa                                                                             <  5:1220562/150‑1 (MQ=255)                  tGCATGAAAACTCTAGCGGAGTATACTGCGCATCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGa                                                                             <  4:559046/150‑1 (MQ=11)                  tGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTAcatcgg                                                                             >  5:1044782/1‑145 (MQ=17)                  tGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCGCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGa                                                                             <  2:893255/150‑1 (MQ=255)                  tGCATGAAAAAGATATAGGTGTATATTGCGCATCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGa                                                                             <  7:132113/150‑1 (MQ=11)                  ggCATGAAAAAGCCTGCGCAGTATACTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGAc                                                                            <  3:853535/150‑1 (MQ=11)                   gaaTGAAATCTCGTGCGGAGTATAGTGCGCATACACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGAc                                                                            <  6:1547333/148‑1 (MQ=11)                   gCATGAATACTCGATCGGAGTATTCTGGGCTTCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGAcg                                                                           <  6:197785/151‑1 (MQ=11)                   gCATGAACAAGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGTTGGAGCAGCCTGGTAGCGCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCACCGCAACCAATTAAAATTTGATGAAGT‑ACAGCAGTACGGTGAc                                                                            <  4:618730/150‑1 (MQ=18)                   gCATGAAAACGCGATGGTAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAACTAAAAGCAGTACGGTGAcg                                                                           <  4:104517/151‑1 (MQ=14)                   gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGAcg                                                                           <  3:585357/151‑1 (MQ=255)                   gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGAcg                                                                           <  4:1569367/151‑1 (MQ=255)                   gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGAcg                                                                           <  4:750840/151‑1 (MQ=255)                   gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGGCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGAc                                                                            <  8:609724/150‑1 (MQ=255)                    cATGAAAACGCTAGGGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGAcgc                                                                          <  5:1573265/151‑1 (MQ=14)                     aTGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGAcgcg                                                                         <  1:1109363/151‑1 (MQ=255)                     aTGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGAcgcg                                                                         <  8:1662537/151‑1 (MQ=255)                      agaaaaatatatcgtaGTATATTGCGATTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCgg                                                                        <  6:70926/136‑1 (MQ=11)                       gAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCggg                                                                       <  1:1058833/151‑1 (MQ=255)                       gAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCgg                                                                        <  5:732766/150‑1 (MQ=255)                       gAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCgg                                                                        <  8:1279172/150‑1 (MQ=255)                       gAAAACCCGAGCGGAGTATAGTCCGCATCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCgg                                                                        <  1:1333118/150‑1 (MQ=11)                        aaaaCGAGTGGGCAGTATAGAGCGCTTACCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCgggg                                                                      <  7:986003/151‑1 (MQ=11)                        aaaaCAAGATCGGAGTATAGTCCGCATCCTGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCgggg                                                                      <  3:1465745/151‑1 (MQ=1)                          accgcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTg                                                                    <  4:687317/149‑1 (MQ=255)                           aCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTgg                                                                   <  7:763141/151‑1 (MQ=255)                            cgcgAGCGGAGTATAGTGCGCCTCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGAc                                                                            <  1:1479073/141‑1 (MQ=11)                            cgcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATGAAAATTTCTGGAGGT‑AAAGCAGTACGGTGAcggcgcggccg                                                                  >  5:1472026/1‑142 (MQ=17)                             gcgAGCGGAGTATAGTGCGCATCGACGGACGAGGGGTGGAGCAGCCTTGTAGCTCGTCGGGCTCATAACCCGAAGGGCGTCGGTGCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGGAGTACGGTGACGCGGGGTGGag                                                                 <  8:1508186/151‑1 (MQ=18)                              cgAGCGGAGTATAGTGATCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGGCGGTCTCATAACCCGAATGTCGGCGGTGCAAATCCGTCCCCCGCAACCAATAAAAATTTGATGAAGT‑AAAGCAGTACCGTGACGCGGGGTGGagc                                                                <  6:299514/151‑1 (MQ=17)                              cgAGCGGAGTATAGGGAGCATCCACTGACGCGGTGTGGAGCAGCCTGGTAGCTCTTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGagc                                                                <  1:138021/151‑1 (MQ=255)                              cgAGCGGAGTAGAGTGCGCATCCACGGACGCGTGCTGGAGCAGCCTGGTAGCTCGGCGGGCTCATAACCCGAAGGTCGTCGGTGCAAATCCGCCCCCCGCAACCAATTAAAAATTGATGAAGT‑AAAGCAGTACGCTGACGCGGGGTGGagc                                                                <  1:1118235/151‑1 (MQ=18)                              cgAGCGGAGTAGACGGCGCATCCACCGACGCGGGGTGGAGCAGCCTGGTAGCTCGGCGGTCTCATAAGCCTAATGTCGTCGGTTCAAATCCGTACCCCGCAACCAATTAAACTTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGagc                                                                <  1:420922/151‑1 (MQ=17)                                ggtggAGTATAGTGGGCTTGGCCGGACGAGGGTGGGTGGAGGCTGGTAGCTTGTATGGCTCATATGCCGAAGGGCGTAGGTTCAAATCCGGCCCGCGCAACCAATTAAAATTTGATGAAGT‑CACGCAGTACGGTGACGCGGGGTGGagcag                                                              <  2:1509810/148‑1 (MQ=9)                                 gCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGagcagc                                                             <  6:1640773/151‑1 (MQ=255)                                   gtaGTATAGTGCGCATCCACGGACTCTTGTTGCAGGAGCCTGGTCGCTCGTCGGTCTCATAACCCGAAGGTCGTCGCTTCTAATACGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCt                                                           <  4:1049851/149‑1 (MQ=17)                                    gAGTATAGTGCGCATCCCCGGACGCGGGGAGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGGTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTg                                                          <  5:1231089/151‑1 (MQ=255)                                                         aCGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTGGGATCAACTCCGGCCCGCGCAACAAAAGACAATAGAATGAAGG‑ACCGCAGGACGGTGCCGCGGGGAGGAGCAGCCTGGTAGCGCGTCGGGcgccggcc                                     >  7:362816/1‑144 (MQ=2)                                                                                        ggatCATATCCCTAAGGTGGTCTGTTGCAATCCTTCCCCCGGAACCAATTATAAATTGTTGTATT‑TACGGATTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCg      <  6:75134/148‑1 (MQ=11)                                                                                            aaTAACCCGAAGGTCTGCGGTTCAAATCCGGCCGCCGCAAGGAATTAAAATTTGCTGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGATCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAAc  <  7:1204703/150‑1 (MQ=14)                                                                                                        gtcgtcGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAgg                              >  1:165125/1‑111 (MQ=255)                                                                                                                                                         |                                                                                           CGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAAC  >  NC_000913/2947330‑2947571

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