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

CTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGA  >  NC_000913/366382‑366673                                                                                                                                                  |                                                                                                                                                    cTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGGGCCCGGGGGGgttt                                                                                                                                               >  2:341556/1‑151 (MQ=255)  tGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTC‑TTTGGCGGCCGGGGGGttttt                                                                                                                                              >  2:645364/1‑151 (MQ=255)  tGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTTCGTA‑TTTGGGGCCAGGGGGGgtttt                                                                                                                                              >  1:609781/1‑151 (MQ=255)       tGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTGCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGTGTGGTttttct                                                                                                                                            <  1:673583/148‑1 (MQ=255)        gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTttttcttttc                                                                                                                                        <  1:38561/151‑1 (MQ=255)         ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTGGGGGCCCGGGGGGGttttctttt                                                                                                                                         >  2:440613/1‑149 (MQ=255)            aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGTA‑TTTGGGGGCCGGGGGGGttttctttttccc                                                                                                                                     >  1:923841/1‑146 (MQ=255)                  tgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGTGGT‑TTTGGCGCCCGGGGGGGttttcttttttcccgaggg                                                                                                                               >  1:324086/1‑140 (MQ=255)                   gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAc                                                                                                                             <  1:618324/151‑1 (MQ=255)                    agCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACg                                                                                                                            <  2:295358/151‑1 (MQ=255)                     gCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGACCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACgg                                                                                                                           <  1:267787/151‑1 (MQ=255)                        aaCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACggg                                                                                                                          <  1:175707/149‑1 (MQ=255)                         aCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCaa                                                                                                                       <  1:533199/151‑1 (MQ=255)                            cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACa                                                                                                                     <  2:201683/150‑1 (MQ=255)                             acaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGGCCGGGGGGGTTTTT‑TTTTCACCAGGGAGGAGGGcca                                                                                                                       >  1:466541/1‑145 (MQ=255)                               aTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTC‑TTTGGCGCCAGGGGGGGTTTTTTTTTCACCAGTGAGAAGGGCaa                                                                                                                       >  1:920619/1‑145 (MQ=255)                                 tAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGa                                                                                                                <  2:805668/150‑1 (MQ=255)                                 tAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑ATGGGCGCCCGGGGGGGTTTTCTTTTCATCAGTGCAACGGGCAACCGcc                                                                                                                  >  1:103333/1‑147 (MQ=255)                                  aaTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGAtt                                                                                                              <  1:117367/151‑1 (MQ=255)                                         ttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCtt                                                                                                        >  2:176510/1‑150 (MQ=255)                                            cgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGCGCCCGGGGGGTTTTTTTTTTTACCCGTGAGACGGGGAACAGCTGATTTCCCTTcccc                                                                                                    >  2:126439/1‑151 (MQ=255)                                             gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCg                                                                                                   <  2:985207/151‑1 (MQ=255)                                               tCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGT‑TTGGGGGCCCGGGGGGGTTTTCTTTTCTCCCGGGGGACGGGGAACAGCCGAGTTCCCTTCACCGcc                                                                                                 >  1:38354/1‑151 (MQ=255)                                                cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGTA‑TTGGGGGCCAGGGGGGGTTTTTTTTTCACCAATGAGACGGGGAACAGCTGATTGCCCTTCAACGcc                                                                                                 >  2:363771/1‑150 (MQ=255)                                                cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGcc                                                                                                 <  1:1187/150‑1 (MQ=255)                                                    gCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTgg                                                                                              <  1:473632/149‑1 (MQ=255)                                                     cccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAAAAGCTGATTGCCCTTCACCGCCTGGccc                                                                                           <  2:350145/151‑1 (MQ=255)                                                      ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCt                                                                                          <  1:583346/151‑1 (MQ=255)                                                                        aCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGc                                                                        <  2:347860/151‑1 (MQ=255)                                                                           tGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTTCGGA‑TTGGGCGCCCGGGGGGGTTTTTTTTTTACCAGCGAGACGGGCAACAGCTGATTGGCCTTCACCGCCTGGCCCTGGGGGAGGTGGAGGAAGCgg                                                                      >  1:515552/1‑150 (MQ=255)                                                                              cGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCa                                                                  <  2:580605/151‑1 (MQ=255)                                                                                tGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGCGGCCGGGGGGGTTTTTTTTTTACCAGGGGGAAGGGCCACAGCTGATTGCCCTTCACCCCCCGGCCCTGGGAGAGGTGCAGCAAGCGGGCCAc                                                                 >  2:387037/1‑150 (MQ=255)                                                                                   cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGGGAGAAGGGGAAAAGCTGATTGCCCTTCACCGCCCGGCCCCGGGAGGGGTGCAGCAAGCGGGcccccccg                                                             >  2:78486/1‑145 (MQ=255)                                                                                   cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCCACAGCTGATTGCCCTTCACCCCCTGGCCCTGAGAGAGGTGCAGCAAGCGGGCCACGc                                                               >  2:993130/1‑149 (MQ=255)                                                                                    aGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGA‑TTTGGCGCCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAAAAGCTGAATTCCCTTCACCCCCTGGCCCTGAGGGAGTTGCAGCAAACGGTCCACGCTgg                                                            >  2:210619/1‑151 (MQ=255)                                                                                       tGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGCGCCCGGGGGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCCGGCCCTGGGAGAGTTGCAGCAAACGGTCCACCCTGGttt                                                         >  2:443041/1‑151 (MQ=255)                                                                                         cATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGA‑TTGGGGGCCCGGGGGGGTTTTTTTTTCACCAGTGAGACGGGCAACAGCTGGTTTCCCTTCACCCCCCGGCCCTGAGAGAGTTGCAGCAAACGGGCCACGCTGGgtt                                                         >  2:471442/1‑149 (MQ=255)                                                                                             aaTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTGGGCGCCCGGGTGGGTTTTTTTTTCACCAGTGAGACGGGCAACAGCTGGTTGCCCTTCTCCGCCTTGCCCTGGGAGGGGTGTAGCAAAGGGGCCACGCCGGTTTGcccc                                                    >  1:819941/1‑150 (MQ=255)                                                                                             aaTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGGGCCAGGGGGGGTTTTCTTTTCACCCGTGAGACGGGCCACAGCTGATTGCCCCTCACCCCCCGGCCCCGAGAGAGTTGCAGCAAGCGGGCCCCGCTGGGTTGccccc                                                   >  2:622125/1‑150 (MQ=255)                                                                                                 aaTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagca                                                <  2:408156/150‑1 (MQ=255)                                                                                                    cGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCCACAGCTGATTGCCCTTTACCCCCTGGCCCTGAGAGAGGTTCAGCAAGCGGGCCACGCTGGTTTGCCCCCGCAgg                                              >  2:869866/1‑149 (MQ=255)                                                                                                    cGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGc                                             <  1:812329/150‑1 (MQ=255)                                                                                                       ccAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGCGCCAGGGGGGTTTTTTTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCCCCCGGCCCTGAGAGAGTTGCAGCAAGCGGGCCACGCTGGTTTGCCCCAGCAGGCGaa                                          >  2:284894/1‑150 (MQ=255)                                                                                                          aCGCGCGGGGAGAGGCGGTTTGCGGA‑ATTGGGGCCCGGGGGGGTTTTCTTTTTACCCGTGAGACGGGGAACAGGTGATTGCCCTTCACCGCCTGGCCCCGGGAGAGGTGCAGCAAGCGGGCCACGCCGGTTTGGCCCCGCAGGGGGAAATc                                      >  1:39097/1‑151 (MQ=255)                                                                                                            gcgcgGGGAGAGGCGGTTTTCGTA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGGTGCAGCAAGCGGTCCACGCCGGGTTGCCCCAGCAGGCGGAAAAcc                                     >  2:776881/1‑150 (MQ=255)                                                                                                              gcgGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGt                                  <  1:819488/151‑1 (MQ=255)                                                                                                               cgGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGtt                                 <  1:98023/151‑1 (MQ=255)                                                                                                                ggggAGAGGCGGTTTGCGTA‑TTGGGCGCCCGGGGGGTTTTTCTTTTCACCACCGAGACGGGCAACAACTGCTTGCCCTTCCCCGCCTGGCCCCGAGAGAGTTGCAGCAAGCGGTCCCCGCTGGTTTGCCCCAGCCCGCGAAAATCCTCttt                                >  1:585466/1‑151 (MQ=255)                                                                                                                 gggAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTg                               <  1:345454/151‑1 (MQ=255)                                                                                                                      agGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtgg                           >  2:379874/1‑150 (MQ=255)                                                                                                                            tttGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                     >  1:813236/1‑150 (MQ=255)                                                                                                                               gCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcgg                 <  1:614810/151‑1 (MQ=255)                                                                                                                                  tA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg              >  2:282775/1‑151 (MQ=255)                                                                                                                                   a‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg              <  2:140415/150‑1 (MQ=255)                                                                                                                                        ggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat         <  2:875760/151‑1 (MQ=255)                                                                                                                                            ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa       >  2:705006/1‑149 (MQ=255)                                                                                                                                            ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAAc      >  2:254154/1‑150 (MQ=255)                                                                                                                                            ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa     <  1:650205/151‑1 (MQ=255)                                                                                                                                             cAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa     >  2:860517/1‑150 (MQ=255)                                                                                                                                             cAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACAt    <  2:594165/151‑1 (MQ=255)                                                                                                                                              aGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa     >  2:337116/1‑149 (MQ=255)                                                                                                                                               gggTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGa  <  2:565425/151‑1 (MQ=255)                                                                                                                                               gggTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGa  <  2:46592/151‑1 (MQ=255)                                                                                                                                                  |                                                                                                                                                    CTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGA  >  NC_000913/366382‑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: ‑