Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	366,526	0	T	G	39.3%	89.5 / 21.0	58	T329P (ACC→CCC)	lacI	lactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (8/26);  new base G (22/0);  total (32/26)
Fisher's exact test for biased strand distribution p-value = 3.66e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGA  >  NC_000913/366379‑366673                                                                                                                                                     |                                                                                                                                                      aGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGggtggt                                                                                                                                                    <  1:776343/151‑1 (MQ=255)    cTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGGTTTTGGCGGCAGGGGGGgttt                                                                                                                                                 >  1:96740/1‑151 (MQ=255)     tGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGGATTGGGGGCCCGGGGGGgtttt                                                                                                                                                >  2:955054/1‑151 (MQ=255)       gggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCACTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTttttct                                                                                                                                              <  1:382674/151‑1 (MQ=255)        ggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTttttctt                                                                                                                                             <  2:1021114/151‑1 (MQ=255)            ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTCGTATTTGGCGCCAGGGGGGGttttctttt                                                                                                                                           >  2:587431/1‑149 (MQ=255)              tAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACATGAATTCCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGGTTGTGCGTTGGGGGCCAGGGAGGTgtttcttttttc                                                                                                                                        >  1:164263/1‑147 (MQ=255)               aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTGGGATTTGGCGCCAGGGGGGGTTTTCTTTTTAc                                                                                                                                        >  2:387416/1‑149 (MQ=255)               aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTTTTTGGGGCCCGGGGGGGttttttttttccc                                                                                                                                       >  1:875452/1‑146 (MQ=255)                aTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTTGGGGCCAGGGGGGGTTTTTTTTTTAcc                                                                                                                                       >  2:313090/1‑149 (MQ=255)                aTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAg                                                                                                                                     <  2:540225/151‑1 (MQ=255)                    gtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTTGGGACTTGGTGCCCCGGGGGGTTTTTTTTTTACCCAgg                                                                                                                                    >  1:238381/1‑148 (MQ=255)                       agCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTGTGTGTCGTTGGGGCCAGGGGGGGttttctttttccgacagag                                                                                                                                 >  1:439994/1‑138 (MQ=255)                            aCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCaa                                                                                                                        <  2:993857/151‑1 (MQ=255)                                     aaTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGGTTGCGTATTGGGCGCCCGGGTGGGTTTTCTTTTTCCCCGGGAGAC‑GGGCAACAG‑CCGAtt                                                                                                              >  1:141918/1‑151 (MQ=255)                                       ttgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGGTTTGGGGGCCAAGGGGGGTTTTTTTTTCACCCGGGAGAC‑GGGGAACCG‑cgtatt                                                                                                              >  1:651802/1‑144 (MQ=255)                                         gcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTTGGGGCGAGGGGGGGTTTTTTTTTTACCCGAGAGGC‑gggggcccgctgcttg                                                                                                              >  1:281133/1‑135 (MQ=255)                                            ttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCtt                                                                                                        <  1:594255/150‑1 (MQ=255)                                               cgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCAc                                                                                                     <  1:310226/150‑1 (MQ=255)                                                    aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTCGTGATGGGGGCCCGGGGGGGTTTTTTTTTTACCCGGGAGAC‑GGGGAACCG‑CGGATTGCCCTTCcccccc                                                                                                 >  1:405247/1‑149 (MQ=255)                                                                ccAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTTGGGGCCCGGGGGGGTTTTTTTTTTTCCCGTGGGGC‑GGGGAAAAG‑CTGGTTTGCCTTCCCCCCCTGGGCCTgaggg                                                                                     >  1:265460/1‑147 (MQ=255)                                                                        gAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTCGTAATGGGCGCCAGGGGGGGTTTTCTTTTTACCCGGGAGGC‑GGGCCACCGAGGGCTTGCCC‑TCTCCGCCTGGCCCTCAGAGAGTTtca                                                                              >  1:686345/1‑148 (MQ=255)                                                                          aaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAg                                                                         <  2:599211/151‑1 (MQ=255)                                                                            ccTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGc                                                                        <  1:517242/150‑1 (MQ=255)                                                                             cTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGTGTATTTGGCGCCAGGGGGGGTTTTTTTTTCACCAGGGAGAC‑GGGGAACAG‑CTGATTGCCCTTTACCGCCTTGCCCCGGGAGAGGTGGAGCAAACgg                                                                      >  2:189237/1‑151 (MQ=255)                                                                                 cGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCa                                                                  <  2:692611/151‑1 (MQ=255)                                                                                      cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGGATTTGGGGCCAGGGGGGTTTTTCTTTTCACCAGGGAGAC‑GGGCAACAG‑CTGATTGCCCCTCCCCCCCTCGGCCTGAGAGAGGTGCGGCAAGAGGTGCCCGc                                                               >  1:46945/1‑149 (MQ=255)                                                                                       aGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGCTTTGGTATTGGGTGCCAGGGGGGGTTTTCTTTTCGCCAGGGAAAC‑GGGCAACAG‑CTTGTTCCACTTCACCGCCCGGGCCTGAGAGAGGTGAAGCAAGCGGTCCAccttgg                                                            >  2:220076/1‑146 (MQ=255)                                                                                            cATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTCGTATTTGGCGCCAGGGGGGGTTTTTTTTTTACCAGTGAGAC‑GGGGAACAG‑CTGATTGCCCTTCACCCCCTGGCCCTGAGAGAGGTGCAGCAAGCGGGCCACCCTGGttt                                                         >  2:943590/1‑149 (MQ=255)                                                                                            cATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGATTGCGTTTTTGGTGGCCGGGGGGGTTTTCTTTTTATCCGTGGGAA‑GGGCAACAG‑AAGATTGCCCCTCACCGCCTCGGCCTGAGAGGGTTGCGGCAAGAGGGccc                                                                  >  1:407099/1‑139 (MQ=255)                                                                                                aaTGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCca                                                   <  2:680020/151‑1 (MQ=255)                                                                                                  tGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTTCGGTTTTGGCGCCCGGGGGGGTTTTTTTTTTACCCGTGAGAC‑GGGCAACAA‑CTGATTGCCCTTCACCGCCCGGCCCTGGGAGAGTTGCAGCAAGCGGTCCACGCTGGGTTTcccc                                                    >  2:740973/1‑148 (MQ=255)                                                                                                   gAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagca                                                <  1:800022/151‑1 (MQ=255)                                                                                                    aaTCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCCGGGGGGTTTTTCTTTTCACCAGGGAGGC‑GGGCCACAA‑GTGATTGCCCTTCACCGCCTGGCCCTGAGGGAGGTTTAGCAAGCAGTCCACCCTGGGTTGTCCcacc                                                 >  1:311805/1‑147 (MQ=255)                                                                                                     aTCGGCCAACGCGCGGGGAGAGGC‑GGATTTCGTGTTTGGGGGCCGGGGGGTTTTTTTTTTTCCCCGTGAGACGGGGCAACAC‑TGGTTTTCCCTTC‑CCCCCCCGCCCTGAGAGAATTGCGACCAATGGGCCACGCACGGTTGcccc                                                    >  2:903045/1‑145 (MQ=255)                                                                                                       cGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCg                                            <  2:869500/151‑1 (MQ=255)                                                                                                       cGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCg                                            <  2:640530/151‑1 (MQ=255)                                                                                                        ggCCAACGCGCGGGGAGAGGC‑GGTTGGCGTATTGGGCGCCAGGGTGTTTTTTCTTTTCTCCAGTGAGAC‑GCGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGgaga                                           <  2:516126/151‑1 (MQ=255)                                                                                                          ccAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaaa                                         <  1:803380/151‑1 (MQ=255)                                                                                                            aaCGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTGCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCGCCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAAt                                       <  1:713414/151‑1 (MQ=255)                                                                                                             aCGCGCGGGGAGAGGC‑GGTTTGGGGGTTGGGGGCCCCGGGGGTTTTTCTTTTCACCAATGAGAC‑GGGCAACAG‑CTGATTGGCCTTCACCCCCTGGCCCTTTGAGAGTTTCAGCAAGCGGGCCACGCCGGGTTGCCCCACCAGGCggaca                                        >  1:640834/1‑145 (MQ=255)                                                                                                                        gagGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtgg                           <  1:664996/151‑1 (MQ=255)                                                                                                                        gagGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtgg                           <  2:455535/151‑1 (MQ=255)                                                                                                                         agGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggt                          <  1:46583/151‑1 (MQ=255)                                                                                                                          ggC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtg                         >  1:877002/1‑151 (MQ=255)                                                                                                                            c‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATggtggt                       >  2:308070/1‑151 (MQ=255)                                                                                                                              ggTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGtt                      <  2:384566/151‑1 (MQ=255)                                                                                                                              ggTGTGCGCATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGtt                      >  1:808661/1‑151 (MQ=255)                                                                                                                               gTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                     <  1:871342/151‑1 (MQ=255)                                                                                                                                  tGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg                  >  2:837680/1‑151 (MQ=255)                                                                                                                                  tGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAc                   <  1:83011/150‑1 (MQ=255)                                                                                                                                    cGTATTGGGCGCCATGTTGTTTTTTCTTTTCACCAGTGAGAC‑GGCCAACAG‑CTGATTGCCCTTCACCGCCTGTCCCTGAGATAGTTGCAGCAAGCTGTCCACTCTCGTTTGCCCCAGCACTCGAAAATCCTATTTAAttctccttaactt                 <  1:323393/150‑13 (MQ=255)                                                                                                                                         tGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg             <  2:949153/149‑1 (MQ=255)                                                                                                                                         tGGGCGCCAGGGCGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg             >  1:357718/1‑149 (MQ=255)                                                                                                                                          gggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGat           <  1:249401/150‑1 (MQ=255)                                                                                                                                               ccAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa       >  2:588308/1‑149 (MQ=255)                                                                                                                                                 aGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGGTATaaaa     >  2:552171/1‑147 (MQ=255)                                                                                                                                                  gggTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGGTATaa       >  1:779901/1‑146 (MQ=255)                                                                                                                                                  gggTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGa  <  2:598286/151‑1 (MQ=255)                                                                                                                                                   ggtggtTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGGTATAAc      >  1:224092/1‑146 (MQ=255)                                                                                                                                                    gtggtTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatatta       >  2:789421/1‑142 (MQ=255)                                                                                                                                                     |                                                                                                                                                      AGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGC‑GGTTTGCGTATTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAC‑GGGCAACAG‑CTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGA  >  NC_000913/366379‑366673

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