Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	366,526	0	T	G	40.7%	152.0 / 20.8	81	T329P (ACC→CCC)	lacI	lactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (9/39);  new base G (33/0);  total (42/39)
Fisher's exact test for biased strand distribution p-value = 8.35e-15
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGA  >  NC_000913/366377‑366673                                                                                                                                                          |                                                                                                                                                     aaaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGggtg                                                                                                                                                     <  2:460813/151‑1 (MQ=255)         gggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTttttc‑t                                                                                                                                            <  2:816543/151‑1 (MQ=255)             gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTttttc‑tttt                                                                                                                                         <  2:1257420/150‑1 (MQ=255)               cTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGGTTTCGTA‑TTTGGG‑G‑CCCGGGGGGGgtttc‑tttttcc                                                                                                                                      >  1:419742/1‑148 (MQ=255)                tAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTCGTA‑TTTGGG‑G‑CCCGGGGGGGTTTTT‑TTTTTAcc                                                                                                                                     >  1:551572/1‑151 (MQ=255)                 aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGCTTTCGCG‑TTGGGC‑GCCCAGGGGGGTTTTT‑‑TTTTCAcc                                                                                                                                     >  1:1329232/1‑150 (MQ=255)                  aTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTTCGGAA‑TTGGGG‑G‑CCGGGGGGGTttttt‑ttttccc                                                                                                                                      >  1:783419/1‑147 (MQ=255)                   tgagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTT‑TTTGGCGG‑CCGGGGGGGTttttt‑ttttccca                                                                                                                                     >  1:611953/1‑146 (MQ=255)                    gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGG‑G‑CCAGGGGGGGTTTTT‑TTTTTACCAg                                                                                                                                   >  2:1031835/1‑149 (MQ=255)                     agtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTTGTA‑TTTGGC‑G‑CCAGGGGGGGTTTTT‑TTTTTACCAGt                                                                                                                                  >  2:1223797/1‑149 (MQ=255)                     agtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTga                                                                                                                                <  1:832892/151‑1 (MQ=255)                      gtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTgag                                                                                                                               <  1:13555/151‑1 (MQ=255)                       tgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGGGGG‑TTTGGG‑G‑CCCGGGGGGGTTTTT‑TTTTTAcccg                                                                                                                                   >  1:159773/1‑144 (MQ=255)                       tgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGGT‑TTTGGC‑G‑CCCGGGGGGGttttc‑ttttttccaatg                                                                                                                                 >  2:826699/1‑140 (MQ=255)                        gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTCGGA‑TTTGGC‑G‑CCCGGGGGGTTTTTT‑TTTTCACCAGTg                                                                                                                                 >  1:310162/1‑147 (MQ=255)                         agCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACg                                                                                                                            <  1:124057/151‑1 (MQ=255)                         agCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACg                                                                                                                            <  1:134627/151‑1 (MQ=255)                          gCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACgg                                                                                                                           <  1:93181/151‑1 (MQ=255)                          gCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACg                                                                                                                            <  1:1176065/150‑1 (MQ=255)                             aaCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTCGTC‑GTTGGG‑G‑CCCGGGGGGGTTTTT‑TTTTTACCCGCgaga                                                                                                                              >  1:734937/1‑145 (MQ=255)                             aaCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGc                                                                                                                         <  2:1217652/150‑1 (MQ=255)                                tCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGGTTTTGGA‑TTGGGA‑G‑GCCGGGGGGGTTTTT‑TTTTTACCCGGGAGAAGGGGaa                                                                                                                       >  1:1024632/1‑149 (MQ=255)                                 cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGG‑G‑CCACGGGGGGTTTTT‑TTTTTACCAGTGGGACggggca                                                                                                                       >  1:42013/1‑145 (MQ=255)                                 cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAAc                                                                                                                      <  2:1075760/149‑1 (MQ=255)                                 cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAg                                                                                                                    <  2:253773/151‑1 (MQ=255)                                  acaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGc                                                                                                                   <  2:1088168/151‑1 (MQ=255)                                  acaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGc                                                                                                                   <  2:1331910/151‑1 (MQ=255)                                     ttAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGGATTTGGGG‑G‑CCGGGGGGGTTTTT‑‑TTTTCACCAGTGGGACGGGGaa                                                                                                                       >  1:1074230/1‑144 (MQ=255)                                       aaTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGAtt                                                                                                              <  1:1047810/151‑1 (MQ=255)                                        aTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCCGGGGGGGTTTTT‑TTTTTACCAGTGAGACGGGCAAAAGGTGAtt                                                                                                              >  2:16701/1‑150 (MQ=255)                                          tgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGG‑G‑GCAGGGGGGGTTTTC‑TTTTCACCCCTGGGGAGGGGAACAGCTGAATggc                                                                                                           >  1:668594/1‑149 (MQ=255)                                          tgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGGTTGCGCA‑TTGGGG‑G‑GCCGGGGGGTTTTTC‑TTTTTCACAGTGCGAGGGGCGACAggcgatagc                                                                                                            >  1:628653/1‑142 (MQ=255)                                              ttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGGTTGCGGT‑TTTGGA‑G‑GCACAGGGGGGTTTG‑TTTTCACCAAAGAGACAGGCAACAGCtttgttggcct                                                                                                        >  1:770964/1‑140 (MQ=255)                                                      aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGGA‑TTTGGC‑GCCCGGGGGGGTTTTTT‑TTTTCCCCCG‑GAGACGGGGAAACACTGATTGCCCCTcccc                                                                                                    >  1:847051/1‑146 (MQ=255)                                                           ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCt                                                                                          <  2:598043/151‑1 (MQ=255)                                                              cTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGCATCTGCCAACGCGCGGGGAGAGGC‑GGTGTAACGA‑TTGGGA‑G‑TCAGCGGGGTTTTTC‑TTTTCCCCAGGGAGGCGGGCAACAGCTGAATTGCCTTCTCCCCATCGCCCttgg                                                                                       >  1:173894/1‑148 (MQ=255)                                                                ttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgag                                                                                       <  2:644622/149‑1 (MQ=255)                                                                  ccAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGagagag                                                                                   <  1:1171761/151‑1 (MQ=255)                                                                   cAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGGGTA‑TTTGGG‑G‑CCCGGGGGGGTTTTT‑TTTTTACCAGTGAGACGGGCCACAGCTGATTGCCCTTCACCCCCTGGCCCTGagagag                                                                                   >  1:720146/1‑150 (MQ=255)                                                                     gTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGGGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGtt                                                                                 <  2:1273412/150‑1 (MQ=255)                                                                       cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgca                                                                              <  2:346397/151‑1 (MQ=255)                                                                       cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGCTTTTGTG‑AGTTGC‑G‑CCAGGGTTGGTTTTT‑TTTTCACCAATTAAATGAGCGACAGCCTAGTGCCCATCTCCGCCTCGGCCtcgg                                                                                       >  2:944842/1‑139 (MQ=255)                                                                               cTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTCGTA‑TTTGGC‑G‑CCAGGGGGGGTTTTT‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCCGGCCCCGAGAGAGTTGCAGCAAGCgg                                                                      >  2:1032207/1‑151 (MQ=255)                                                                                tGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTTGGC‑G‑CCCGGGGGGTTTTTT‑TTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGGTGCAGCAAGCgg                                                                      >  2:851561/1‑150 (MQ=255)                                                                                    gTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTCGTA‑TTTGGG‑G‑CCCGGGGGGGTTTTT‑TTTTCACCAGTGAGACGGGCAACAGCTGATTTCCCTTCCCCCCCTGGCCCTGAGAGAGGTGCAGCAAGCGGGCCAc                                                                 >  2:1204748/1‑151 (MQ=255)                                                                                      gCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGTGTATTTGGGG‑G‑CCGGGGGGGTTTTT‑‑TTTTTACCCGGGAGAAGGGGAACAGCTGATTGCCCTTTACCCCCTGGCCCTGAGGGAGGTTCAGCCAGCGGGCCACGc                                                               >  2:410156/1‑151 (MQ=255)                                                                                          gCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGGGGA‑TTTGGG‑G‑CCCGGGGGGGTTTTCTTTTTCCCCAG‑GAGAAGGGGAACAGCCGATTGCCCCTCACCCCCCCGCCCCGCGAGAGTTGCCGCaaaggggcccccctgg                                                            >  1:447210/1‑136 (MQ=255)                                                                                              cATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTCGTA‑TTTGGG‑G‑CCAGGGGGGGTTTTT‑TTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCCCCTGGCCCTGAGAGAGTTGCAGCAAGCGGGCCACGCCGGttt                                                         >  2:759128/1‑149 (MQ=255)                                                                                                ttAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGTGTA‑TTGGGC‑G‑CCAGGGGGGTTTTTT‑TTTTCACCCGGGAAACCGGCCAAAGCTGAATGCCCTTCACCGCCCGGCCCTGGGAGAGTTGCCGGAAGCGGTCCACGCTGGGTTGccc                                                     >  1:732142/1‑151 (MQ=255)                                                                                                            ccAACGCGCGGGGAGAGGC‑GGTTTGTGGA‑TTGGGG‑G‑CCCGGGGGGGTTTTT‑TTTTTACCAGTGAGACGGGGCACACATGAATGCCCCTCACCGCCCGGCCCCGGGAGAGGTGGAGCCAACGGGCCACCCTGGTTTGGCCCCGCAGGCg                                            >  2:1049000/1‑148 (MQ=255)                                                                                                            ccAACGCGCGGGGAGAGGC‑GGCTTGCGGC‑TTGGGC‑G‑CCCGGGGGGGTTTTT‑TTTTTCCCACAGAGGCGGGCAACAGCCGCTTGCCCCTCCCCCCCCGGGCCCGAGGGAGTTGCAGCAAGCGGGACCCGCCGGGTTGCGCCCGCCGGc                                             >  1:841634/1‑147 (MQ=255)                                                                                                              aaCGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaaaa                                        <  2:612038/150‑1 (MQ=255)                                                                                                              aaCGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAAt                                       <  2:238667/151‑1 (MQ=255)                                                                                                                 gcgcgGGGAGAGGC‑GGTTTGCGTA‑TGGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGGGAGGCGGGCAACAGCTGTTTGCCATTCGCCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAAGcc                                     <  1:1261958/150‑1 (MQ=255)                                                                                                                 gcgcgGGGAGAGGC‑GGGTTTCGTA‑TTGGGC‑G‑CCCGGGGGGGTTTTC‑TTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCCCCTGGCCCTGAGAGAGTTGCCGCAAGCGGTCCACGCTGGTTTGCCCCCGCAGGCGCAAATcc                                     >  1:73085/1‑150 (MQ=255)                                                                                                                  cgcgGGGAGAGGC‑GGTTTTCGTA‑TTGGGG‑G‑CCAGGGGGGGTTTTC‑TTTTCACCAGTGAGACGGGGAACAGCTGATTGCCCTTCACCCCCCGGGCCTGGGAGAGGTGCAGCAAGCGGGCCACGCCGGTTTGCCCCCGCAGGCGAAAAtac                                     >  1:1022585/1‑147 (MQ=255)                                                                                                                  cgcgGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTg                                   <  2:696905/151‑1 (MQ=255)                                                                                                                   gcgGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGt                                  <  2:461639/151‑1 (MQ=255)                                                                                                                         agagGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtg                            <  2:645602/151‑1 (MQ=255)                                                                                                                         agagGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtg                            <  1:898724/151‑1 (MQ=255)                                                                                                                          gagGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtgg                           <  2:1024767/151‑1 (MQ=255)                                                                                                                            ggC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtg                         <  2:757881/151‑1 (MQ=255)                                                                                                                            ggC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggt                          <  1:528584/150‑1 (MQ=255)                                                                                                                             gC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtgg                        >  2:413926/1‑151 (MQ=255)                                                                                                                                  tttGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCATTGAGACGGGAAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGATGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggggg                        >  1:1295512/1‑147 (MQ=255)                                                                                                                                   ttGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTaa                    <  2:991791/150‑1 (MQ=255)                                                                                                                                   ttGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTaa                    <  1:885198/150‑1 (MQ=255)                                                                                                                                     gCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg                  <  2:131011/150‑1 (MQ=255)                                                                                                                                     gCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg                  <  2:735041/150‑1 (MQ=255)                                                                                                                                      cGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc                >  1:283798/1‑151 (MQ=255)                                                                                                                                       gTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc                <  1:1010281/150‑1 (MQ=255)                                                                                                                                        tA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg              >  1:974876/1‑151 (MQ=255)                                                                                                                                        tA‑TTGGGA‑G‑CGAGGGGGGTTTTTC‑TTTTCACCCGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGATGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGAGGTTAAcggcgg              >  1:259189/1‑151 (MQ=255)                                                                                                                                             gggC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGat           <  1:901608/150‑1 (MQ=255)                                                                                                                                              ggC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg             >  2:1118977/1‑147 (MQ=255)                                                                                                                                              ggC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat         <  1:646886/151‑1 (MQ=255)                                                                                                                                               gc‑g‑cCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGAtata        <  2:406808/151‑1 (MQ=255)                                                                                                                                                c‑g‑cCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat         >  2:1301854/1‑149 (MQ=255)                                                                                                                                                  g‑cCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAAc      <  2:486309/151‑1 (MQ=255)                                                                                                                                                    ccAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa     <  1:1193520/151‑1 (MQ=255)                                                                                                                                                     cAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa     >  2:21229/1‑150 (MQ=255)                                                                                                                                                      aGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGGTATAACAt    >  1:355588/1‑150 (MQ=255)                                                                                                                                                         gtggtTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACAtaa  >  1:479073/1‑147 (MQ=255)                                                                                                                                                         gtggtTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACAt    >  2:192981/1‑147 (MQ=255)                                                                                                                                                          |                                                                                                                                                     AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTA‑TTGGGC‑G‑CCAGGGTGGTTTTTC‑TTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGA  >  NC_000913/366377‑366673

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