Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	535,918	0	A	G	43.5%	82.4 / 13.3	63	E70G (GAG→GGG)	hyi	hydroxypyruvate isomerase
Reads supporting (aligned to +/- strand):  ref base A (6/29);  new base G (27/0);  total (33/30)
Fisher's exact test for biased strand distribution p-value = 4.39e-12
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGC‑GC‑GGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGC‑GCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCAC  >  NC_000913/535769‑536066                                                                                                                                                      |                                                                                                                                                         ttttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTAGGAGGGGGGGGa                                                                                                                                                          >  2:913297/1‑150 (MQ=255)   ttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAgc                                                                                                                                                        <  1:98405/150‑1 (MQ=255)   ttGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAgc                                                                                                                                                        <  1:526090/150‑1 (MQ=255)               cgcAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGC‑GC‑GGTATTGCCTg                                                                                                                                         <  2:269849/151‑1 (MQ=255)               cgcAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGC‑GC‑GGTATTGCCTg                                                                                                                                         <  2:13493/151‑1 (MQ=255)                gcAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGGCTGGGCGGCCGGGGGGC‑GC‑ggcagt                                                                                                                                              >  2:756866/1‑141 (MQ=255)                     gtgGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGC‑GC‑GGTATTGCCTGTATTcc                                                                                                                                   <  1:588476/151‑1 (MQ=255)                     gtgGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACCCGGCGGCGGGGGGGC‑GC‑GGGAATGCCTGTATTcc                                                                                                                                   >  2:880658/1‑151 (MQ=255)                       ggTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGC‑GC‑GGTATTGCCTGTATTCCTg                                                                                                                                 <  1:71659/151‑1 (MQ=255)                         ttttCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGC‑GC‑GGGATTGCCTGTATTTCttg                                                                                                                                >  2:916470/1‑148 (MQ=255)                          tttCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGC‑GC‑GGTATTGCCTGTATTCCTGGc                                                                                                                               <  1:1146965/150‑1 (MQ=255)                            tCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGC‑GC‑GGTATTGCCTGTATTCCTGGCCg                                                                                                                             <  1:490053/150‑1 (MQ=255)                               cggcgTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGC‑GC‑GGTATTGCCTGTATTCCTGGCCGTgaa                                                                                                                         <  1:917555/151‑1 (MQ=255)                                   gTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGATGGGCGGGGGGGGGGC‑GC‑GGGGTTTTCTGTATTTCttg                                                                                                                                >  1:23021/1‑138 (MQ=255)                                        aTTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGTGC‑GC‑GGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGttt                                                                                                               <  1:823267/151‑1 (MQ=255)                                         tttATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGC‑GC‑GGGATTTGCTGTATTTCTTGGCGTGAAG‑AAGAGgtt                                                                                                               >  2:853166/1‑150 (MQ=255)                                               tttCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGAATGGGCGGGGGGGGGGG‑GG‑GGGCTTGCCCGTATTCCTGGCCGTGAAG‑AAGAGCTTTGGGGTg                                                                                                        >  2:763797/1‑151 (MQ=255)                                                ttCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGC‑GC‑GGGATTGCCTGTATTTCTGGGCGTGAAG‑AAGAGGTTTGGGGTgg                                                                                                       >  2:940179/1‑151 (MQ=255)                                                ttCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGC‑GC‑GGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATg                                                                                                        <  1:1001043/150‑1 (MQ=255)                                                  ccTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGC‑GC‑GGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGc                                                                                                      <  1:1099891/150‑1 (MQ=255)                                                  ccTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGG‑GC‑GGGGCTGCCTGTATTTCTGGCCGGGGAG‑GAGAGGTTCgggg                                                                                                          >  2:587034/1‑145 (MQ=255)                                                   cTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGC‑GC‑GGGATTGCCTGTATTTCTGGCCGTGGAG‑AAGAGGTTTGGGGTGGCg                                                                                                     >  2:413603/1‑150 (MQ=255)                                                       tGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGGCTGGGGGGCGGGGGAGC‑GC‑GGGATTGACTGTATTTCTGGCCGTGAAG‑AAGAGTTTCggggtggggt                                                                                                    >  1:147881/1‑147 (MQ=255)                                                        gACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGC‑GC‑GGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTagcag                                                                                               <  1:910543/151‑1 (MQ=255)                                                        gACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGC‑GC‑GGGATTGGCTGTTTTCCTGGCCGTGAAG‑AAGAGGTTTGGGGTgg                                                                                                       >  2:670861/1‑143 (MQ=255)                                                          cTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGCGGCGGGGGGGC‑GC‑GGTATTGCCCGGATTTCCGGCCGGGGAG‑AAAAGTTTCGGGGTGGCGGag                                                                                                  >  2:1194067/1‑146 (MQ=255)                                                           tACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGGGGCGGGGGGGC‑GC‑GGGATTTCCTGTATTCCTTGGCGGGAAG‑AAGAGTTTTGGGGTGGGGTagca                                                                                                >  2:375142/1‑147 (MQ=255)                                                            aCGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGC‑GC‑GGGATTGCCTGTATTTCTGGGCGGGAAG‑AAGGGGTTCGGGGTGGCGTagcagca                                                                                             >  2:27850/1‑149 (MQ=255)                                                             cGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGTGGAGGGCGGCGGGGGAGC‑GC‑GGGAATTACTGTATTCGTGCGCGTGGAGAAAGAGTTCCGGGAggggtagacaca                                                                                              >  1:1128356/1‑139 (MQ=255)                                                               aCATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGC‑GC‑GGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGAt                                                                                        <  2:845928/151‑1 (MQ=255)                                                                 aTTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGGGGGGGGGG‑GC‑GGGATTGCCTGTATTTCTGGCCGTGAAG‑AAGAGGTTTGGGGTGGGGTAGCAGCAGCGGtt                                                                                       >  2:705573/1‑150 (MQ=255)                                                                        aaTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGC‑GC‑GGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGc                                                                               <  2:800924/151‑1 (MQ=255)                                                                         aTTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGGGC‑GC‑GGGATTGCCTCTTTTCCCTGCCGTGAAG‑AAGAGGTTTGGGAGGGCGGAGCAGCAGCGAGTAGTTATGcc                                                                              >  2:173629/1‑151 (MQ=255)                                                                         aTTAAAACATGTGCTGGCGAGGAATAAACTCGAACATACGCTGCACAACTTACCGGCGGGTGACTGGGCGGCGGGGGAGCGGG‑GGTATTCCCTGATTCCCTTGCCG‑GATG‑AAGAGGTTCGGGGGGGCGCtgcagcccagatcagttatg                                                                                >  1:706946/1‑135 (MQ=255)                                                                           tAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGGGGGGGGGG‑GG‑GGGATTACTTGTTTTCCCGGCCGTGAAG‑AAGAGGTTCGGGGGGGGGGAGCAGCGGGGATTTGGTCTGccc                                                                             >  1:614383/1‑150 (MQ=255)                                                                             agaCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGC‑GC‑GGTATGGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTAGTTATGCCCg                                                                            <  1:118794/147‑1 (MQ=255)                                                                               aCATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGC‑GC‑GGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTgc‑g                                                                       <  2:998897/151‑1 (MQ=255)                                                                                   gtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGC‑GC‑GGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGC‑GCTgg                                                                   <  2:27720/151‑1 (MQ=255)                                                                                   gtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGC‑GC‑GGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGC‑GCTgg                                                                   <  1:1124332/151‑1 (MQ=255)                                                                                   gtgCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGC‑GC‑GGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGC‑GCTgg                                                                   <  2:600713/151‑1 (MQ=255)                                                                                      cTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGGGGGGGGGC‑GG‑GGGATTGTCTGTATTCCTGGGCGGGGAG‑AAGAGGTTGGGGGGGGGGTGGGAGGAGGGATTTGGTATGCCCGTGC‑GCTggggg                                                                >  1:917111/1‑149 (MQ=255)                                                                                       tGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGGGGCTGGGGGGCGGGGGGGG‑GC‑GGGATTTCCTGTTTTCCTGGGCGGGAAG‑AAGAGGTTTGGGGTGGGGGAGCAGCAGCGGTTCGTTTTGCCCGGGC‑GCTggg                                                                  >  2:225924/1‑148 (MQ=255)                                                                                         gCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGAATGGGCGGGGGGGGGGA‑GCGGGAATTCCCTTTTTTCC‑CGCCGTGGAG‑AAGAGTTTTGGGGGGGGGTAGGAGCAGCGATTCGTTATTCCCCGGG‑GCTgggggat                                                              >  1:1193177/1‑146 (MQ=255)                                                                                           gAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGAATGGGCGGCGGGGGGGC‑GC‑GGGATTGCCTGTATTCCTGGGCGGGGAG‑AAGAGGTTCGGGATGGCGGAGCAGGAGCGATTCGTTATGCCCGTGG‑GCTGGGgaataa                                                            >  1:559100/1‑150 (MQ=255)                                                                                             gTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGTGCGGCGGGGGAGC‑GC‑GGTATCGCCTGTATTCCTGGCCGTGAAG‑AAGATTTTCGGGATGGCTTAGCAGCAGCGATTCGTTATGCCCGTGC‑GCTGGGTAATaaaaa                                                         <  2:600008/151‑1 (MQ=255)                                                                                              taataaACTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGC‑GC‑GGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGGTTTGGGGTGGCGGAGCAGCAGCGATTCGTTATGCCCGTGC‑GCTGGGTAATaaaaa                                                         >  2:422960/1‑150 (MQ=255)                                                                                                ataaACTCGAACATACGCTGCACAATTTACCGGCGGGTGGACGCGCGGCGGGGGGGG‑GC‑CGGTTTTCCTGCATTTCTCGCCGGGTAG‑AAGAGGTTCGGGGTGGCGCGGCATCAGCGATTCGCTATGCCCACGA‑GCCGGGTAATaaaaaa                                                        >  2:299382/1‑149 (MQ=255)                                                                                                ataaACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGC‑GC‑GGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTCTGCCCGTGC‑GCTGGGTAATaaaaa                                                         <  1:1175381/148‑1 (MQ=255)                                                                                                ataaACTCGAACATACGCTGCACAATTTACCGGCGGGTGAATGGGCGGCGGGGGGGG‑GC‑GGGATTGCCCGTATTCCCGGCCGTGAAG‑AAGAGGTTCTGGATGGCGTAGCAGCCGCGATTTGTTATGCCCGTGC‑GCTTGGGAATaaaaaaa                                                       >  2:216137/1‑150 (MQ=255)                                                                                                   aaCTCGAACATACGCTGCACAATTTACCGGCGGGTGACCGGGCGGGGGGGGAGC‑GC‑GGGATTTTCTTTATTCCTGGCCGTGAAG‑AGGGGGTTTGGGGGGGCGTAGCAGCAGCGATTTGTTATGCCCCTGC‑GCTGGGGAATaaaaaaa                                                       >  1:478563/1‑147 (MQ=255)                                                                                                     cTCGAACATACGCTGCACAATTTACCGGCGGGTGGCTTGGCGGCGGGGGAGC‑GC‑GGGATTGCTTTTTTTCCCGGCCGTGGAG‑AAGAGGTTTGGGGAGGGGGAGGAGCAGCGATTTGGTTTGGCCGTGG‑GGTGGGGGAGCAAAAAAATaa                                                   >  1:473146/1‑149 (MQ=255)                                                                                                     cTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGC‑GC‑GGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGAATCGTTATGCCCGTGC‑GCTGGGTAATAAAAAAATTAAc                                                  <  1:521278/150‑1 (MQ=255)                                                                                                     cTCGAACATACGCTGCACAATTTACCGGCGGGGGACTGGGCGGCGGGGGGGC‑GC‑GGGATTGCCTGTATTTCTGGCCGTGAAG‑AAGAGGTTTGGGGTGGCGTAGCAGCAGCGATTTGTTATGCCCGTGG‑GCTGGGGAATAAAAAAAATaa                                                   >  2:881273/1‑149 (MQ=255)                                                                                                           cATACGCTGCACAATTTACCGGCGGGGGACTAGGCGGGGGGGGGGG‑GC‑GGGATTTTCTTTGTTTCTCGCCCGGGTG‑AAGCGGTTTTGGATGGCGGAGCCCGAGGGATTTGTTTTTACCTTGCGGCGGGGTATAAAAAAAATTaa                                                   >  2:270848/1‑144 (MQ=255)                                                                                                            aTACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGC‑GC‑GGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGC‑GCTGGGTAATAAAAAAATTAACTGTCTGGt                                          <  2:801235/151‑1 (MQ=255)                                                                                                            aTACGCTGCACAATTTACCGGCGGGGGGCTTGGGGGGGGGGGGGG‑GC‑GGGAATGCTTCTATTCCCGGGGCTGGAG‑GAGAAGTTTGGGGTGGGGGAGCAGGAGCGATTTGGTTTTGCCGGGC‑GCTGGGGGATAAAAAAAATaaa                                                  >  1:1013812/1‑142 (MQ=255)                                                                                                            aTACGATGGATAATTTACCGGCGGGTGACTGGGGGGGGTGTGATC‑GC‑GGTTTTGCCTGTATTCCTGGCGGGAAAG‑AAGAGTTTCGGGATGGCGTAGCAGCATCAATTCTTTATGCCAATGC‑CCGTGGTAATAAAAAAATTTCCTGTCTgtt                                          <  1:816833/151‑3 (MQ=255)                                                                                                                   gCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGC‑GC‑GGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGC‑GCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTaa                                    <  1:664323/150‑1 (MQ=255)                                                                                                                   gCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGC‑GC‑GGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTAGGCCCGTGC‑GCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTaaa                                   <  1:291812/151‑1 (MQ=255)                                                                                                                        aTTTACCGGCGGGTGATTGGGAGGCGGGGGGGC‑GC‑GGCATAATCTGTATTCCTGGCCGGGAAG‑AAGGGGTTCGGGTTGGCGTAGCAGCAGCGATCTGTTTTGACCGTGC‑GCTGGGTAATAAAAAAATTAACTGTCTGGTTGGGAAAACGcc                              >  1:1092307/1‑151 (MQ=255)                                                                                                                               ggcggGTGACTGGGCGGCGGGGGAGC‑GC‑GGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGC‑GCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGt                       >  2:989029/1‑151 (MQ=255)                                                                                                                                    gTGACTGGGCGGCGGGGGAGC‑GC‑GGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGC‑GCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTcag                  <  1:875208/151‑1 (MQ=255)                                                                                                                                            gcggcgGGGGAGC‑GC‑GGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGC‑GCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAAc          <  1:565393/151‑1 (MQ=255)                                                                                                                                                    ggAGC‑GC‑GGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGC‑GCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCAc  <  2:440590/151‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                         TTTTGAGAAAGCTGCGCAGTGTGGTTTTCGCGGCGTTGAATTTATGTTTCCTTATGACTACGACATTGAAGAATTAAAACATGTGCTGGCGAGTAATAAACTCGAACATACGCTGCACAATTTACCGGCGGGTGACTGGGCGGCGGGGGAGC‑GC‑GGTATTGCCTGTATTCCTGGCCGTGAAG‑AAGAGTTTCGGGATGGCGTAGCAGCAGCGATTCGTTATGCCCGTGC‑GCTGGGTAATAAAAAAATTAACTGTCTGGTCGGTAAAACGCCGGCTGGTTTCAGCAGTGAACAGATTCAC  >  NC_000913/535769‑536066

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