Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,947,481	1	.	A	26.5%	150.5 / 92.0	68	intergenic (+18/‑16)	metZ/metW	tRNA‑Met/tRNA‑Met
Reads supporting (aligned to +/- strand):  ref base . (7/43);  new base A (0/18);  total (7/61)
Fisher's exact test for biased strand distribution p-value = 1.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.57e-01

CGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGG  >  NC_000913/2947330‑2947612                                                                                                                                                         |                                                                                                                                    ggtgttAAAAATAAGGTTGAATGAAAACGAGATTGGTGTATATTGCGATTCCTCGGACGCGGGGTGGAGCAGCCTGGTATCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGt                                                                                                                                      <  5:590296/146‑1 (MQ=14) ggggtTAAAAATAAGGATGCATGAAAACTCGAGCGGAGAATATTGCGCATCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGt                                                                                                                                      <  2:203689/150‑1 (MQ=17)  ggggTAAAAAAACGGCTGCATGAAAACCCCAGCGGAGTATACTCCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑a                                                                                                                                    <  3:1072169/151‑1 (MQ=17)  gggcTAAAAAAAACGGTGCATCAAAAACCGAGAGGAGTATAGTGCGCATCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑a                                                                                                                                    <  8:1240812/151‑1 (MQ=14)   ggtaaaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑a                                                                                                                                    <  2:915320/147‑1 (MQ=255)    gctAAAATACAGCTTGCATGAATAAGCGAGCGCAGTATAGTGCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aaa                                                                                                                                  <  8:110579/149‑1 (MQ=14)     gTAAAAAAACGCTGGCATGAAAAAGCGAGAGATGAATACTCCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aaaa                                                                                                                                 <  8:1038746/151‑2 (MQ=11)     gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑aaa                                                                                                                                  <  4:1088878/151‑1 (MQ=255)     gTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑aaa                                                                                                                                  <  4:252587/151‑1 (MQ=255)     gTAAAAAAAAGGGTGCATGAATACCCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aaaa                                                                                                                                 <  6:1003628/151‑2 (MQ=17)      tAAAAAAAAGGTTGCATGAAAACGCTAGCGGAGTATAGTCCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aaaa                                                                                                                                 <  8:25693/150‑2 (MQ=255)      aaaaacaaaGGCTGAATGAATACCAGATAGTAGTATATTGCGCATACCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aaaag                                                                                                                                <  5:382297/150‑3 (MQ=1)       aaaaaaaaGGTTGCATGAAAACGCGAGCGTAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aaaagc                                                                                                                               <  6:1374503/151‑4 (MQ=255)       aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAg                                                                                                                                 <  5:28962/150‑1 (MQ=255)       aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAg                                                                                                                                 <  5:999633/150‑1 (MQ=255)       aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAg                                                                                                                                 <  8:997139/150‑1 (MQ=255)       aaaaaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGc                                                                                                                                <  3:219763/151‑1 (MQ=255)         aaaaaaGGGTGAATGAAAAACCCAGGGTAGTATAGTCCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑aaaagcag                                                                                                                             <  1:1567372/151‑6 (MQ=11)          aacaaGGGTGAATGAAATAGCCATCGGAGTAGAGTCCGACTCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGt                                                                                                                             <  8:391640/148‑1 (MQ=11)          aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGt                                                                                                                             <  4:297831/151‑1 (MQ=255)          aaaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGt                                                                                                                             <  2:2217/151‑1 (MQ=255)          aaaaaGGTTGAAGGAAAAATAGAGCGGAGTATAGTCCGCATCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGt                                                                                                                             <  4:193631/151‑1 (MQ=11)           aaaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTa                                                                                                                            <  4:565105/151‑1 (MQ=255)            aacgggtGCATGAAATATCGATCGGTGTAGTGAGCGATTCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGGCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTAc                                                                                                                           <  4:183795/145‑1 (MQ=11)            aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGTGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTAc                                                                                                                           <  5:1338954/151‑1 (MQ=255)            aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTAc                                                                                                                           <  1:1696085/151‑1 (MQ=255)            aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTAc                                                                                                                           <  8:1676019/151‑1 (MQ=255)            aaaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTAc                                                                                                                           <  7:314795/151‑1 (MQ=255)            aaaGGTTGCATGAAAAAGAGAGCGGTGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTAc                                                                                                                           <  1:1028983/151‑1 (MQ=14)             aaGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACg                                                                                                                          <  2:1590606/151‑1 (MQ=255)              cgggggaatcaaaaCGCCATCGCTGTATTCTCCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACgg                                                                                                                         <  2:1287531/141‑1 (MQ=11)              cggatcaaTGAAAAAGAGATGGGAGTATTTTGCGCATCCTCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACgg                                                                                                                         <  4:396081/144‑1 (MQ=1)               ggTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGt                                                                                                                        <  6:1639750/151‑1 (MQ=255)                ggtGCACGAAATCGCGATCGTAGTATAGTGGGCTTCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTg                                                                                                                       <  2:1517835/149‑1 (MQ=11)                 gtGCACGAAAAAGCGAGCGGAGTATACTGCGCATCCCCGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGa                                                                                                                      <  4:695193/150‑1 (MQ=14)                  tGCATGAAAACTCTAGTGTTGTATTTTGCGCATCCTGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGa                                                                                                                      <  6:575058/150‑1 (MQ=11)                  tGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGa                                                                                                                      <  8:623164/150‑1 (MQ=255)                  tGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGa                                                                                                                      <  2:1179951/150‑1 (MQ=255)                  tGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTAGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TAAAATTTTATTAAGt                                                                                                                                      >  5:817232/1‑134 (MQ=255)                  ggCAGGAAAACGCGTGGGGAGTATTGTTCGCCTCCAGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGAc                                                                                                                     <  2:1584724/150‑1 (MQ=11)                   gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAATTTACATCAACA‑AGCGCAGTACTGTGAcg                                                                                                                    >  3:103107/1‑151 (MQ=18)                   gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGAcg                                                                                                                    <  7:1020801/151‑1 (MQ=255)                   gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGAc                                                                                                                     <  8:1081260/150‑1 (MQ=255)                   gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGAc                                                                                                                     <  8:1379276/150‑1 (MQ=255)                   gCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGAcg                                                                                                                    <  8:1291116/151‑1 (MQ=255)                   gCATGAAAAAGCTAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGAcg                                                                                                                    <  2:1086112/151‑1 (MQ=14)                    cATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGAcgc                                                                                                                   <  2:1740116/151‑1 (MQ=255)                    cATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGAcgc                                                                                                                   <  5:440318/151‑1 (MQ=255)                     atgaaaaatcttgCGGAGTATAGTCCGACTCCCGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGAcgcg                                                                                                                  <  2:1743734/139‑1 (MQ=11)                     aTGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGAcgcg                                                                                                                  <  8:1501173/151‑1 (MQ=255)                       gAAAACCCGTGGGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCgg                                                                                                                 <  6:256442/150‑1 (MQ=11)                       caatactagtgtggAGTATAGTGCGCCTCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCggg                                                                                                                <  1:436864/139‑1 (MQ=11)                        aaaaCTCGAGCGGAGTATAGTGGGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCgggg                                                                                                               <  4:751870/151‑1 (MQ=14)                         aatggcgATCGGAGTATTGTTCGCATGCAGGGTCGCGGTGTTTAGGAGGCTGGTAGCTGGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCGGGCCCGCGCATGCAATTTAAATTTGATGAAGT‑ATAGCTGTAGGGTGAGGGGGGGt                                                                                                              <  5:1667558/147‑1 (MQ=9)                         aaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGt                                                                                                              <  4:1397786/151‑1 (MQ=255)                         aaaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGt                                                                                                              <  4:1405813/151‑1 (MQ=255)                          aaCGGGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTGGTCGGTTCAAATCCGGCCCCCTCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTg                                                                                                             <  3:700546/151‑1 (MQ=255)                          aaCGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGGCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTg                                                                                                             <  6:1034596/151‑1 (MQ=255)                           aCGCGAGCGGAGTATAGTGCCCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTgg                                                                                                            <  2:1626874/151‑1 (MQ=255)                            cgcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGa                                                                                                           <  1:1606405/151‑1 (MQ=255)                             tcctgCGGAGTATAGTGCGCATCCAGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGTAAAAGCAGTACGGTGAcg                                                                                                                    <  3:1704502/137‑1 (MQ=11)                             gcgAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATAAAAATTTATTAAATA‑AAAGCAGTACGGTGAcg                                                                                                                    >  7:1700159/1‑141 (MQ=17)                                aGCGTAGTATTCTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TCAAATTTGATGAAGT‑a                                                                                                                                    <  1:4906/121‑1 (MQ=11)                                aGCGGAGTATAGTTCGCATCCACTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCTGTCTCATAACCCGAAGGTCGTCGGTTCAAATACGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGagcag                                                                                                       <  6:202412/151‑1 (MQ=255)                                aGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTGGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGagcag                                                                                                       <  3:968020/151‑1 (MQ=255)                                   ggAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAA‑TAAAATTTTATtaaat‑aaaagcag                                                                                                                             >  6:176889/1‑120 (MQ=14)                                   ggAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGGTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCt                                                                                                    <  2:199790/151‑1 (MQ=255)                                    gAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGGCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTg                                                                                                   <  7:385971/151‑1 (MQ=255)                                     aGTATAGTGATAATCCACTGAAGAGTGGTGGAGCAGCCTGGTCGCTCTTCGGTCTAATAAGCCTAAGGTGGTCGGTTCAAATACGTACCCCGCAAACAATTAAACTTTGATGAAGT‑AAAGCAGTACTGTGACGCGGGGTGGAGCAGCCTgg                                                                                                  <  5:619657/151‑1 (MQ=9)                                     aGTAGAGTGCTCATCCACTGAATCTGTTTGGTTTAGCGTGTTAGCTTGTGGCTCTCATAACCCGAATCTCGGCGGTTGATATCCTGCCCCCGGAACCAATTAACCTTTGATGAAGT‑CAACCACTACCCTGACGCGGGGTGGAGCAGCCTg                                                                                                   <  1:695891/150‑1 (MQ=9)                                      ttttaggtcgcATCCTCCGACGAGGTGTGGATCAGGGTGGTAGCTCGTCGGTATCATAACCCGAATGTCGTCGGTTGAAATCATGCCCGGGCAACCAAATAAAATTTGATGAAGT‑AAAGCACTACGGTGACGCGGGGTGGAGCAGCCTgg                                                                                                  <  4:314521/142‑1 (MQ=9)                                        attgTGCTCATCCACGGAAGCTGGGTGGAGCAGCCTGGTAGCTCGTCGGGATCATAAGCCGAATTGCGTCGGGTCAAATCAGTCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCACTACGGTGACGCGGGGTGGAGCAGCCTGGTa                                                                                                <  4:1156151/147‑1 (MQ=17)                                             gcgcATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCtcgt                                                                                          <  3:100719/151‑1 (MQ=255)                                             ccggctCCTCGGCCGCGGGTTGGATCAGCGTTGTATCTCGGCTGTCTAATAACCCTAAGGGCGGCGCTTCAAAACCGGCCCCCGCAACCAATTAAATATTGCAGAATT‑AACCGCCTACGGGGACGCGGGGTGGAGCAGCCTGGTAGCtc                                                                                            <  3:1537193/144‑1 (MQ=2)                                                    cACGGACGAGGGGTGGAGCAGGCTGGTATCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCTAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTc                                                                                   <  3:975763/151‑1 (MQ=255)                                                                                    gtttgtatCATATGCCTAAGGTCTTCGCTTGAAATCCGTCCGCCGGAAGGAATTAAATTTTGTTGAAGT‑ATAGGACTCCGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGcc                                                   <  6:286971/144‑1 (MQ=2)                                                                                                              ggTTCAAATCCGGCCCCCGCAACCAATTAAATTTTGATCAAGT‑AACGCAGTACGGTGACGCGGGGTGTAGCAGCCTGGTAGCTCGTCGGGCTCATGACCCGAAGGTCGTCGGTGGAAATCCGGCCCCCGCAACCAAtgatagtagaccctc                         >  8:573373/1‑137 (MQ=11)                                                                                                                   aaaTACTGCCGCCGGAATCAATTAAAATTTGATGAATT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTaaaa                    <  7:686573/151‑1 (MQ=17)                                                                                                                     aTCCGGCCCCCGCAACCAATTAAAATTTTTTTAATT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCgg                                                     <  5:703211/116‑1 (MQ=255)                                                                                                                     aTCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCgg                                                     >  2:705000/1‑116 (MQ=255)                                                                                                                                  aaCGAATTATATTTTGTTGAAGT‑TTCGCATTACGGGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGAcgc     <  1:283879/151‑1 (MQ=14)                                                                                                                                     cAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGcc                                                   >  7:268855/1‑102 (MQ=40)                                                                                                                                     cAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGCTCGTCGCTTCAAATCCGGCCCCCGCAACCAATCAAATAAGAAGAACAAAAAGCAGTAATGTGACGAggg  >  7:390548/1‑151 (MQ=17)                                                                                                                                       aTTAAATTTTTATTACGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGccccc                                                <  5:1615905/103‑1 (MQ=11)                                                                                                                                       aTTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGccccc                                                >  2:1619163/1‑103 (MQ=39)                                                                                                                                                         |                                                                                                                                    CGGGGTAAAAAAAAGGTTGCATGAAAACGCGAGCGGAGTATAGTGCGCATCCACGGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATTAAAATTTGATGAAGT‑AAAGCAGTACGGTGACGCGGGGTGGAGCAGCCTGGTAGCTCGTCGGGCTCATAACCCGAAGGTCGTCGGTTCAAATCCGGCCCCCGCAACCAATCAAATTTGATGAAGTAAAAGCAGTACGGTGACGCGGG  >  NC_000913/2947330‑2947612

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: ‑