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%	88.7 / 14.4	60	T329P (ACC→CCC)	lacI	lactose‑inducible lac operon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (9/26);  new base G (24/0);  total (34/26)
Fisher's exact test for biased strand distribution p-value = 2.15e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGC‑G‑TA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCT  >  NC_000913/366378‑366676                                                                                                                                                         |                                                                                                                                                        aaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGC‑G‑TA‑TTGGGCGCCAGggtgg                                                                                                                                                       <  1:597690/151‑1 (MQ=255) aaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGC‑G‑TA‑TTGGGCGCCAGggtgg                                                                                                                                                       <  2:1132486/151‑1 (MQ=255)  aGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑ATTGC‑G‑CA‑TAGGGTGCCACggggg                                                                                                                                                       >  1:340317/1‑150 (MQ=255)    ccTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGT‑TTTGC‑G‑TA‑TTGGGCGCCAGGGTGGtt                                                                                                                                                     <  2:892400/150‑1 (MQ=255)    ccTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGG‑GT‑TTGGGGGCCGGGGGGGtt                                                                                                                                                     >  1:172283/1‑151 (MQ=255)         ggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGC‑G‑TA‑TTGGGCGTCAGGGTGGTttttctt                                                                                                                                               <  1:302284/151‑1 (MQ=255)           tGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGC‑G‑TA‑TTGGGGGCCCGGGGGGGttttttttt                                                                                                                                             >  1:628500/1‑151 (MQ=255)            gCCTAATGAGTGAGCTAACTCACATTAATTGCGTCGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGC‑G‑TA‑TTGGGCGCCAGGGTGGTttttctttt                                                                                                                                             <  1:115205/150‑1 (MQ=255)               tAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTTTC‑T‑TT‑TTGGG‑GCCGGGGGGGGtttttttttctc                                                                                                                                          >  2:195824/1‑148 (MQ=255)                 aTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGT‑T‑TA‑TTTGGCGCCCGGGGGGGGTTTTTTTTCatccg                                                                                                                                       >  2:573333/1‑147 (MQ=255)                   gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGC‑G‑TA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGt                                                                                                                                      <  1:70164/150‑1 (MQ=255)                    agtgagCTAACTCTCATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGC‑G‑TA‑TTGGGGGCCCGGGGGGGTTTTTTTTTCAcc                                                                                                                                         >  2:975611/1‑146 (MQ=255)                    agtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGC‑G‑TA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTga                                                                                                                                    <  2:212411/151‑1 (MQ=255)                          cTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑CTTGT‑G‑TATTTGGGGGCCGGGGGGGtttttttttttccagggggaaggg                                                                                                                               >  1:1298962/1‑137 (MQ=255)                            aaCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGC‑G‑TA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCa                                                                                                                           <  1:1232016/151‑1 (MQ=255)                                cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGC‑G‑TA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAAc                                                                                                                         <  1:192761/149‑1 (MQ=255)                                cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGC‑G‑TA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACa                                                                                                                        <  1:1215720/150‑1 (MQ=255)                                  caTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGT‑G‑TA‑TTGGGGGCCAGGGGGGTTTTTTTTTTCACCCGCGGG‑ACggggcaca                                                                                                                        >  2:631065/1‑143 (MQ=255)                                   aTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGC‑G‑TA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTg                                                                                                                    <  2:256472/151‑1 (MQ=255)                                   aTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGC‑G‑TA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTg                                                                                                                    <  1:749343/151‑1 (MQ=255)                                   aTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGC‑G‑GATTTGGGGGCCGGGGGGGTTTTT‑TTTTCTCCAGTGAG‑AAGGGCAAAAGCTg                                                                                                                    >  2:629165/1‑151 (MQ=255)                                    ttAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTT‑GTTA‑TTGGGGCCCGGGGGGGTTTTTTTTTTAACC‑GGGGG‑ACGGGCaa                                                                                                                          >  2:1263065/1‑144 (MQ=255)                                     tAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑GTTGG‑G‑TT‑TTTGGCGCCCGGGGGGGTTTTTTTTTTTCCAGTGAG‑ACgggggaca                                                                                                                        >  1:1110883/1‑140 (MQ=255)                                           cgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGC‑G‑TA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCt                                                                                                            <  2:922240/151‑1 (MQ=255)                                                     aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGC‑G‑TA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTTCCAGGGAG‑ACGGGCAACAGCTGATTGTCCTTCACcccc                                                                                                    >  2:203491/1‑149 (MQ=255)                                                             cTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTC‑G‑TA‑TTTGGGGGCCGGGGGGTTTTTTTTTTTACCCAAGAG‑ACGGGGGACAACTGATTTCCCCTCCCCCCCTGGccc                                                                                              >  2:251150/1‑147 (MQ=255)                                                                tCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTC‑G‑GA‑TTGGGGGCCCGGGGGGGTTTTTTTTTCACCCGGGAG‑ACCGGGAACAGCCGATTTCCCCTCACCCCCCGGCCCTgag                                                                                          >  1:61148/1‑148 (MQ=255)                                                                tCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGC‑G‑TA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTTCCAGTGAG‑AGGGGGAACAGCTGATTGCCCCTCACCCCCTGGgcccg                                                                                            >  2:993189/1‑144 (MQ=255)                                                                      cGGGAAACTTGTCGTACCGGCTGCATTAATGAATCGGCCGACACGCGGGGAGAGGCGG‑TTTGC‑G‑TA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGGGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgc                                                                                  <  2:356863/150‑1 (MQ=255)                                                                      cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTC‑G‑TT‑TTTGGGGGCAGGGGGGGTTTTCTTTTTTCCCATGAG‑ACGGGCAACAGCTGATTTCCCTTCACCGCCCGGCCCTgag                                                                                          >  2:360194/1‑142 (MQ=255)                                                                        ggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTG‑G‑TA‑TTTGGGGCCCGGGGGGGGTTTTTTTTTCCCCGTGAG‑GAGGGGAACAGCTGATTGGCCCTCACCCCCCCGCCCTgggaggg                                                                                      >  1:185422/1‑142 (MQ=255)                                                                          aaaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTC‑G‑TG‑TTTGGGGCCCGGGGGGGTTTTTTTTTTTCCCGTGAGAACGGGC‑ACAACCGGTTGCCCCTCCCCCCCTGGCCCTggg                                                                                          >  1:668873/1‑136 (MQ=255)                                                                           aaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGG‑G‑TA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCCGTGAG‑ACGGGCAACAGCTGGTTTGCCCTCACCCCCCGGCCCCGAGAGAGGTgcag                                                                                >  1:350381/1‑147 (MQ=255)                                                                           aaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGC‑G‑TA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAg                                                                            <  1:1021114/151‑1 (MQ=255)                                                                           aaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGC‑G‑TA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAg                                                                            <  1:352838/151‑1 (MQ=255)                                                                           aaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGC‑G‑GG‑TTTGGGGCCCGGGGGGGTTTTTTTTTTTCCAGTGAG‑AAGGGGGAAAGACGATTGGCCCTCACCCCCCTGGCCTGGGAGAGtttggggaa                                                                             >  1:509102/1‑143 (MQ=255)                                                                                       cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTC‑G‑TC‑TTTGGGGCCCGGGGGGGTTTTTTTTTTCCCAGTGGG‑AAGGGGAACAGCTGGTTGCCCCTCACCCCCCCGGCCCGAGAGAGtttgaggcaacgggcc                                                                      >  1:1033439/1‑131 (MQ=255)                                                                                           tGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTC‑G‑GA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGTGAG‑ACGGGGCACAGCTTATTTCCCTTCCCCCCCTGGCCCCGGGGGAGGTGCAGCAAACGGGcccccccgggtt                                                            >  2:654333/1‑141 (MQ=255)                                                                                             cATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGC‑G‑TA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGc                                                          <  2:1134036/151‑1 (MQ=255)                                                                                             cATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGC‑G‑GA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTACCCGGGAG‑ACGGGCCACAGCTGCTTGCCCCTCACCCCCCCGCCCCGGGAGAGGTGCAGCAAGCCGGCCAccccgggtt                                                            >  1:1003704/1‑141 (MQ=255)                                                                                                    gAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGC‑G‑TA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagc                                                    <  2:623356/150‑1 (MQ=255)                                                                                                          gCCAACGCGCGGGGAGAGGCGG‑TTTTC‑G‑GA‑TTTGGCGCCAGGGGGGGTTTTCTTTTCCCCCGTGGG‑ACGGGCGACAGCTGATTGCCCCTCACCCCCTCGCCCTGGGAGAGGTGCAGCAAGCGGTCCACGCCGGGTTGCCCCCGCAgg                                                 >  1:597139/1‑147 (MQ=255)                                                                                                             aaCGCGCGGGGAGAGGCGG‑TTTGC‑G‑TC‑TTTGGCGCCCGGGGGGGTTTTCTTTTTTCCCGCGAG‑AAGGGCCACAGCTTGTTGCCCCTCTCCCCCTTGGCCCGCGAGGGGTGCTGCAAGCGGGCCCCGCCTGGTTGCCCCCGCCGGCGaaa                                            >  1:88707/1‑149 (MQ=255)                                                                                                             aaCGCGCGGGGAGAGGCGG‑TTTGC‑G‑TA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAAt                                          <  2:1110944/151‑1 (MQ=255)                                                                                                                   cgGGGAGAGGCGG‑TTTGC‑G‑TA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGtt                                    <  2:172288/151‑1 (MQ=255)                                                                                                                           ggCGG‑TTTGC‑G‑TA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggt                             >  1:697309/1‑150 (MQ=255)                                                                                                                             cGG‑TTTGC‑G‑TA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCCGCAGGCGAAAATCCTGTTTGAtggtgg                           >  1:1058972/1‑150 (MQ=255)                                                                                                                             cGG‑TTTGC‑G‑TA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATggtggt                          >  1:590858/1‑151 (MQ=255)                                                                                                                               g‑TTTGC‑G‑TA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCCCCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                        <  1:251142/151‑1 (MQ=255)                                                                                                                                     c‑G‑TA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc                   >  1:1112627/1‑151 (MQ=255)                                                                                                                                         tA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg                 >  1:513969/1‑151 (MQ=255)                                                                                                                                          a‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg                 <  2:377998/150‑1 (MQ=255)                                                                                                                                             tGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGa               <  1:1165613/150‑1 (MQ=255)                                                                                                                                              gggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGata             <  1:852737/151‑1 (MQ=255)                                                                                                                                              gggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGata             <  2:509123/151‑1 (MQ=255)                                                                                                                                              gggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGat              <  1:1240034/150‑1 (MQ=255)                                                                                                                                               ggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat            >  1:1149998/1‑151 (MQ=255)                                                                                                                                                 cgcCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGAtata           >  2:361278/1‑150 (MQ=255)                                                                                                                                                   ccAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa        <  1:631042/151‑1 (MQ=255)                                                                                                                                                     aGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa          >  2:1235929/1‑147 (MQ=255)                                                                                                                                                     aGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACAt       >  2:321719/1‑150 (MQ=255)                                                                                                                                                     aGGGCGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGCTTGATGGTGGTTAACGGCGGGGAATaa          >  1:1075320/1‑147 (MQ=255)                                                                                                                                                      gggTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGa     <  2:671145/151‑1 (MQ=255)                                                                                                                                                       ggtggtTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAg    <  2:974916/151‑1 (MQ=255)                                                                                                                                                         tggtTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCt  >  2:167557/1‑151 (MQ=255)                                                                                                                                                         |                                                                                                                                                        AAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGC‑G‑TA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAG‑ACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGCT  >  NC_000913/366378‑366676

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