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

GGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAG  >  NC_000913/366384‑366674                                                                                                                                                 |                                                                                                                                                     ggggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTttttc                                                                                                                                              <  1:754723/151‑1 (MQ=255)       ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCa                                                                                                                                        <  2:86592/151‑1 (MQ=255)       ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCa                                                                                                                                        <  2:37766/151‑1 (MQ=255)             gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGGA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCCGTg                                                                                                                                  >  1:683675/1‑151 (MQ=255)               gtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTgag                                                                                                                                <  1:24430/151‑1 (MQ=255)                tgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGGGCCAGGGGGGGTTTTTTTTTCACCAGTg                                                                                                                                  >  2:231552/1‑148 (MQ=255)                 gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGATTTGGGGG‑TGGGGGGCCGGGGGGGTTTTTTTGTTAACCAggag                                                                                                                                 >  1:553624/1‑146 (MQ=255)                 gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAc                                                                                                                              <  2:683691/151‑1 (MQ=255)                  agCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGTGTA‑TTGGGGGCCCGGGGGGGTTTTCTTTTTACCCGt                                                                                                                                   >  2:148286/1‑145 (MQ=255)                   gCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACgg                                                                                                                            <  2:514436/151‑1 (MQ=255)                   gCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACgg                                                                                                                            <  1:510832/151‑1 (MQ=255)                        cTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTG‑TTGGGGGCCAGGGGGGGTTTTCTTTTTTCCAATGGGACGGGGaa                                                                                                                        >  1:156621/1‑150 (MQ=255)                          cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGGGCCAGGGGGGGTTTTTTTTTTACCAGTGAGAAGGGCaa                                                                                                                        >  2:486811/1‑148 (MQ=255)                          cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAAc                                                                                                                       >  2:801254/1‑149 (MQ=255)                          cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGGA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGGGGGAAggggga                                                                                                                        >  1:904519/1‑145 (MQ=255)                          cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACa                                                                                                                      <  1:585026/150‑1 (MQ=255)                                 aTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGAT‑TTGGGGGGCTGGGGGGGTTTTTTTTTTTCCAGCGGGAGGGGCAACAGCTgtgt                                                                                                               >  2:925637/1‑147 (MQ=255)                                   tgcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGCC‑TTTGGCGCCCGGGGGGTGTTTCTTTTTTCCCGGGGGAGGGGCGACAGCTGact                                                                                                               >  1:941425/1‑146 (MQ=255)                                    gcgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGTGGA‑TTTGGGGCCAGGGGGGGTTTTTTTTTCTCCAGTGAGACGGGCAACAACTGATTGccc                                                                                                           >  2:877120/1‑151 (MQ=255)                                     cgttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑GTTGCGTG‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGGGAGACGGGCAACAGCTGAGTGccc                                                                                                           >  1:461371/1‑150 (MQ=255)                                      gttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGTGGA‑TTTGGCGCCAGGGGGGTTTTTCTTTTTACCAGTGAGACGGGCAACAGCTGATTTCCCtt                                                                                                         >  2:249792/1‑151 (MQ=255)                                           gcTCACTGCCCGCTTTCCAGTCGGGAAATCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGGGCCAGGGGGGGTTTTTTTTTTACCCGTGAGACGGGGAACAGCTGATTGCCCCTCAccc                                                                                                    >  2:918115/1‑150 (MQ=255)                                           gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGGA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAAAAGCTTATTTCCCCTCAccc                                                                                                    >  2:382653/1‑150 (MQ=255)                                            cTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTTGTATTTGGGGGCCGGGGGGGTTTTT‑TTTTTTCCCGTTGGACGGGGAACAACTGATTGCCCTTCTcc                                                                                                     >  2:810548/1‑149 (MQ=255)                                              cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCt                                                                                                 <  2:102129/151‑1 (MQ=255)                                               aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTg                                                                                                <  1:73059/151‑1 (MQ=255)                                               aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGGA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGTGAGAAGGGGGACAACTGGTTGCCCTTTACcccc                                                                                                  >  1:97038/1‑149 (MQ=255)                                                    ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGGGCCAGGGGGGGTTTTCTTTTCACCAGTGAGACGGGCAACAGCCGATTGCCCCTCACCCCCCGGccc                                                                                            >  2:110119/1‑150 (MQ=255)                                                         ttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgaga                                                                                       <  1:332102/150‑1 (MQ=255)                                                            cAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGGGCCAGGGGGGGTTTTTTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCCGGCCCTgag                                                                                        >  2:947170/1‑146 (MQ=255)                                                               tCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTTGGA‑TTTGGGGCCCGGGGGGGTTTTCTTTTTACCAGTGATAAGGGCGACAGCTGGTTTTCCTTCACCCCCCTGCCCCGggagaggt                                                                                  >  2:712222/1‑147 (MQ=255)                                                                cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgca                                                                               <  2:178272/151‑1 (MQ=255)                                                                 gggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGGGCCAGGGGGGTTTTTTTTTTTACCCGTGAGACGGGCAACAGCTGATTTCCCTTTCCCGCCTGGCCCTGAGAGAGtt                                                                                  >  2:251209/1‑147 (MQ=255)                                                                  ggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgcagc                                                                             <  2:61064/151‑1 (MQ=255)                                                                   gAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgcagca                                                                            <  1:9432/151‑1 (MQ=255)                                                                                ccAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGTTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCt                                                               <  1:341917/151‑1 (MQ=255)                                                                                ccAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCt                                                               <  2:865986/151‑1 (MQ=255)                                                                                   gCTGCATTAATGAATCGGTCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTgg                                                             <  1:519215/150‑1 (MQ=255)                                                                                      gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTg                                                         <  1:15638/151‑1 (MQ=255)                                                                                        aTTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTCGTA‑TTTGGCGCCAGGGGGGGTTTTTTTTTTACCCATGAGAAGGGCAACAGCTGATTGCCCTTCACCCCCCGGGCCTGAGAGAGTTTCAGCAAACGGGCCACGCTGGGTTggc                                                       >  1:898965/1‑149 (MQ=255)                                                                                         ttAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTTGTA‑TTTGGGGCCCCGGGGGGTTTTTTTTTTTCCAGTGAGACGGGCAACAGCCGATTGCCCTTCACCCCCCTGCCCTGAGGGAGGTGTCGCAAGAGGGCCACGCTGGTTTGtcc                                                      >  1:805190/1‑151 (MQ=255)                                                                                           aaTGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCca                                                    <  1:382642/151‑1 (MQ=255)                                                                                            aTGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTTGGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGGAAAAGCTGATTGCCCCTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAACGGTCCACGCTGGGTTTccccc                                                    >  2:855326/1‑149 (MQ=255)                                                                                              gAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagc                                                  <  2:382406/150‑1 (MQ=255)                                                                                                    gCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaa                                           <  1:489554/151‑1 (MQ=255)                                                                                                      cAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGGGCCAGGGGGGTTTTTCTTTTTACCAGGGAGACGGGGAACAGCTGATTTCCCCTCACCGCCCGGCCCTGTGAGAGTTGCCGGAAACGGGCCACCCTGGTTTTCCCCCGCAgg                                               >  2:543785/1‑145 (MQ=255)                                                                                                      cAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGaaaa                                         <  2:730574/151‑1 (MQ=255)                                                                                                            gcgGGGAGAGGCGG‑TTTGCGTA‑TTTGGCGCCCGGGGGGGTTTTCTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCCCCTGGCCCTGAGAGAGTTGCAGCAAGCGGGCCACGCTGGTTTGCCCCCGCAGGCGGAAATCCTg                                    >  1:264757/1‑150 (MQ=255)                                                                                                            gcgGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGGCCGGGGGGGTTTTCTTTTCACCAGTGAGACGGGCCACAGCTGATTGCCCTTTACCCCCTCGCCCTGAGGGAGTTGCAGCAAGCGGGCCACGCTGGGTTGCCCCCGCAGGCGAAAATCCTGt                                   >  1:699958/1‑151 (MQ=255)                                                                                                              ggggAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGttt                                 <  2:75650/151‑1 (MQ=255)                                                                                                                        gg‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGtt                       >  1:538303/1‑151 (MQ=255)                                                                                                                         g‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                      <  1:249789/151‑1 (MQ=255)                                                                                                                           tttGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                      >  1:788075/1‑150 (MQ=255)                                                                                                                              gCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg                   <  1:148288/150‑1 (MQ=255)                                                                                                                              gCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcg                   <  1:527002/150‑1 (MQ=255)                                                                                                                                gTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcg                <  1:371259/151‑1 (MQ=255)                                                                                                                                  a‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAACTCCTGTTTGATGGTGGTTAAcggcgg               <  1:493602/150‑1 (MQ=255)                                                                                                                                  a‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg              <  2:97036/151‑1 (MQ=255)                                                                                                                                     tGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGa             <  1:135191/150‑1 (MQ=255)                                                                                                                                           ccAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa        >  1:665345/1‑149 (MQ=255)                                                                                                                                             aGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa      >  2:598976/1‑149 (MQ=255)                                                                                                                                             aGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACAt     >  1:881534/1‑150 (MQ=255)                                                                                                                                             aGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATg    <  1:877102/151‑1 (MQ=255)                                                                                                                                             aGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATg    <  1:810532/151‑1 (MQ=255)                                                                                                                                                gtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAg  >  2:373667/1‑150 (MQ=255)                                                                                                                                                 |                                                                                                                                                     GGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGG‑TTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAG  >  NC_000913/366384‑366674

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