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

AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCA‑CGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACA  >  NC_000913/366377‑366670                                                                                                                                                       |                                                                                                                                                  aaaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGggt                                                                                                                                                   <  2:183980/150‑1 (MQ=255)    gCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTTGGGGCCCggggggg                                                                                                                                                >  2:825456/1‑149 (MQ=255)      cTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGCGCCCGGGGGGgttt                                                                                                                                             >  1:100648/1‑151 (MQ=255)           gTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTttttcttt                                                                                                                                        <  2:418630/151‑1 (MQ=255)           gTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTttttcttt                                                                                                                                        <  2:308160/151‑1 (MQ=255)           gTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTttttcttt                                                                                                                                        <  1:942391/151‑1 (MQ=255)             gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGGGCCCGGGGGGGttttttttt                                                                                                                                       >  1:648347/1‑150 (MQ=255)             gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTttttctttt                                                                                                                                       <  1:375829/150‑1 (MQ=255)              ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGGGCCCGGGGGGGttttttttt                                                                                                                                       >  1:407866/1‑149 (MQ=255)              ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGGGCCCGGGGGGGttttttttt                                                                                                                                       >  1:476716/1‑149 (MQ=255)                 aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCAc                                                                                                                                    <  2:21055/149‑1 (MQ=255)                 aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCa                                                                                                                                  <  2:188754/151‑1 (MQ=255)                  aTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTAcc                                                                                                                                   >  1:82592/1‑149 (MQ=255)                        gagCTAACTCACATTAATTGCGTTGCGCTCACTGGCCGCGTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAGTGAATCTGCCAACGCGAGGGGAGAGGCGGCTTGGGAT‑TAAGGGGCCCGGGGGGTATTTTTTTTAAGCAGGgaga                                                                                                                            >  1:907533/1‑150 (MQ=255)                         agCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACg                                                                                                                          <  1:220415/151‑1 (MQ=255)                          gCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACgg                                                                                                                         <  2:793977/151‑1 (MQ=255)                                tCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTAGTA‑TTTGGCGCCAGGGGGGGTTTTTTTTTTACCAGggggaaggggaa                                                                                                                     >  1:251197/1‑149 (MQ=255)                                tCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCGTA‑TTTGGCGCCCGGGGGGTTTTTCTTTTTCCCAGTGAGACggg                                                                                                                        >  1:554365/1‑146 (MQ=255)                                  acaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGGGGTTTTTTTTTTTACCAGTGAGACGGGGAACa                                                                                                                   >  1:396330/1‑149 (MQ=255)                                      tAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGa                                                                                                              <  1:107719/150‑1 (MQ=255)                                       aaTTGCGTTGCGCTCACTGCCCGCTTTCCAGCCGGGAAACCTGTCGTGCCAGCTGCATTCATGCATCGGCCAACGCGCGGGGAGCGGGGGTTTGCGGA‑TTTTGGGCCGGGGGGGGTTTTCTTTTTACCCGTGAGATGGGCAACTGc                                                                                                                 >  2:280913/1‑146 (MQ=255)                                        aTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGT‑TTGGGGGCCCGGGGGGGTTTTTTTTTTACCCGTGGGACGGGGAACAGCTGAtt                                                                                                            >  2:796737/1‑150 (MQ=255)                                            cgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTTCGGA‑TTGGGCGCCCGGGGGGGTTTTTTTTTTACCACTGCGATGGGCAACAGCTGAATGtcc                                                                                                        >  2:188076/1‑150 (MQ=255)                                             gttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCtt                                                                                                      <  1:261692/151‑1 (MQ=255)                                                    tCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGGGCCCGGGTGGGTTTTTTTTTTTCCCGTGGGAAGGGCAACAACTGATTTCCCTTCACCGcc                                                                                               >  1:672962/1‑151 (MQ=255)                                                          cccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGTT‑TTGGGGGCCAGGGGGGGTTTTTTTTTTCCCAGTGGGACGGGGAAAACCTGATTTCCCTTCCCCCCCCGGccc                                                                                         >  1:337616/1‑151 (MQ=255)                                                           ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGGCCGGGGGGGTTTTT‑TTTTTACCAGTGAGACGGGGAACAGCTGATTTCCCTTCCCCCCCCGGccc                                                                                         >  1:272611/1‑150 (MQ=255)                                                                ttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgaga                                                                                    <  1:341494/150‑1 (MQ=255)                                                                  ccAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCCGTGAGACGGGCAACAGCTGATTGCCCTTCACCCCCCGGCCCTgagag                                                                                   >  1:46512/1‑149 (MQ=255)                                                                   cAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGGGGCAGGGGGGGTTTTCTTTTTACCAGGGAGACGGGGAACAGCTGATTGCCCTTCACCGCCTGGCCCTGagagag                                                                                 >  1:258030/1‑150 (MQ=255)                                                                   cAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGt                                                                                <  2:86269/151‑1 (MQ=255)                                                                     gTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGG‑TTTGGGGCCAGGGGGGGTTTTTTTTTTTCCAAGGGGACGGGGCAAAACTGATTGCCCTTCACCGCCCGGGCCTGGGAGAGtt                                                                               >  1:438492/1‑150 (MQ=255)                                                                       cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGtt                                                                               >  1:840079/1‑148 (MQ=255)                                                                       cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgc                                                                             <  1:25306/150‑1 (MQ=255)                                                                        gggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTGGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTTCCCCTCACCCCCTGGCCCCGggagaggt                                                                               >  1:178554/1‑145 (MQ=255)                                                                         ggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGCGCCCGGGGGGGTTTTCTTTTCACCAGTGGGACGGGCAACAGCTGATTGCCCTTCACCGCCCGGCCCTGAGAGAGGTgca                                                                            >  1:507529/1‑149 (MQ=255)                                                                                  tCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑ATGGGGGCCAGGGTGGGTTTTTTTTTTACCACTGGGACGGGGAACAGCTGATTGCCCCTCCCCCCCCGGCCCTGAGAGAGTTGCAGCAAGCGGGcc                                                                 >  2:756813/1‑151 (MQ=255)                                                                                     tGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTGGGCGCCCGGGGGGGTTTTCTTTTCACCAGTGAGAAGGGGAACAGCTGATTGCCCTTCACCGCCTGGCCCCGAGAGAGGTTCAGCAAGCGGGcc                                                                 >  1:390687/1‑148 (MQ=255)                                                                                          gCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCA‑CGCTGGt                                                        <  1:669750/151‑1 (MQ=255)                                                                                           cTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGGGCCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGGAACAGCTGATTTCCCTTCACCGCCTGGCCCTGAGAGAGGTGCAGCAAGCGGGCCA‑CGCCgggt                                                       >  1:397871/1‑149 (MQ=255)                                                                                            tGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGCGCCCGGGGGGGTTTTTTTTTCACCAGTGAGACGGGGAACAGCTGATTTCCCCTCACCGCCTGGCCCTGAGAGAGTTTCAGCAAGCGGTTCA‑CGCTGGgtt                                                      >  1:187967/1‑151 (MQ=255)                                                                                              cATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGCGGCAGGGGGGTTTTTCTTTTCACCACTGAGAAGGGCAACAACTGATTGCCCTTCTCCGCCTGGCCCTGAGGGAGGTGCAGCAAGCGGGCCA‑CGCTgggt                                                       >  1:59703/1‑146 (MQ=255)                                                                                                 tAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGGGAGAAGGGGAACAGCTGATTGGCCTTCACCCCCCGGCCCCGAGAGGGTTGCAGCAAGCGGGCCC‑CCCCGGGTTGcccc                                                 >  2:832782/1‑151 (MQ=255)                                                                                                 tAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGCA‑TTTGGCGCCAGGGGGGTTTTTCTTTTTACCAGTGAGACGGGCAACCACCTACTGCCCCTCACCCGCTGGCCCTCAGAGGGGGGCAGCAATCGGTCCG‑CGCTGGGTTGcccc                                                 >  1:672931/1‑151 (MQ=255)                                                                                                     gAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGCGCCCGGGTGGGTTTTCTTTTCACCAGTCCGACAGGGCACAGCCGATTGGCCTTTCCCCCCCTGCCCTGCGAGAGTTTCAGCAAACGGGGCA‑CGCCGGGTTGtcc                                                  >  2:155582/1‑146 (MQ=255)                                                                                                      aaTCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCCCCTGGGCCTGAGAGAGGTGCAGCAAACGGTCCACCCCGGGTTTGccccgc                                               >  1:896729/1‑147 (MQ=255)                                                                                                      aaTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCCGGGGGGTTTTTCTTTTCCCCAGGGGGAAGGGGCACCGCTGATTGCCCCTCACCCCCCGGCCCTGAGGGGGGTGCCACAAGCGGGCCA‑CGCTGGGTTGcccc                                                 >  2:238636/1‑146 (MQ=255)                                                                                                       aTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCA‑CGCTGGTTTGCCCCAGCAgg                                           <  2:409601/151‑1 (MQ=255)                                                                                                         cGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCA‑CGCTGGTTTGCCCCAGCAGGc                                          <  2:626569/150‑1 (MQ=255)                                                                                                           gCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTGGGCGCCCGGGGGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCCCCCGGCCCTGAGAGAGTTGCAGCAAGCGGGCCA‑CGCCGGGTTTCCCCCGCAgg                                           >  2:271337/1‑147 (MQ=255)                                                                                                            ccAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCA‑CGCTGGTTTGCCCCAGCAGGCGaaa                                      <  1:402267/151‑1 (MQ=255)                                                                                                               aCGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGGCCA‑CGCTGGGTTTCCCCCGCAGGCGgaaa                                     >  1:879003/1‑149 (MQ=255)                                                                                                                   gcgGGGAGAGGCGGTTTGCGTA‑TTGGGCGCAAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCA‑CGCTGGTTTGCCCCAGCAGGCGAAAATCCTGt                               <  2:852626/151‑1 (MQ=255)                                                                                                                                 tttGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCA‑CGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                  >  1:257836/1‑150 (MQ=255)                                                                                                                                  ttGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCA‑CGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTaa                 <  1:972915/150‑1 (MQ=255)                                                                                                                                  ttGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCA‑CGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGAGGTTAAc                >  1:165616/1‑151 (MQ=255)                                                                                                                                   tGCGTA‑TTGGGCGCCGGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCA‑CGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAc                <  2:687846/150‑1 (MQ=255)                                                                                                                                   tGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCA‑CGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg               >  1:277599/1‑151 (MQ=255)                                                                                                                                    gCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCA‑CGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcgg              <  1:222785/151‑1 (MQ=255)                                                                                                                                     cGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCA‑CGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc             <  2:232963/151‑1 (MQ=255)                                                                                                                                        a‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCA‑CGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg           >  2:299294/1‑150 (MQ=255)                                                                                                                                              gcgcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCA‑CGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGAtata     <  2:867560/151‑1 (MQ=255)                                                                                                                                                   aGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCA‑CGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa  >  1:700076/1‑149 (MQ=255)                                                                                                                                                       |                                                                                                                                                  AAAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCA‑CGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACA  >  NC_000913/366377‑366670

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