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

CCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCA‑GCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAG  >  NC_000913/366381‑366674                                                                                                                                                   |                                                                                                                                                       ccTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTT‑TTTGGGGCCCgggggg                                                                                                                                                      >  2:319439/1‑148 (MQ=255)   tGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTATTTGGGGGCCGGGGGGGtttt                                                                                                                                                  >  1:631616/1‑151 (MQ=255)       gTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTATTTGGGGGCCGGGGGGGtttttttt                                                                                                                                              >  1:964540/1‑151 (MQ=255)       gTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTttttcttt                                                                                                                                             <  1:107815/151‑1 (MQ=255)         gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTT‑TTTGGGGGCACGGGGGGttttttttt                                                                                                                                            >  1:876100/1‑150 (MQ=255)         gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTttttctttt                                                                                                                                            <  1:18748/150‑1 (MQ=255)             aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTAcc                                                                                                                                        >  2:463296/1‑150 (MQ=255)                gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTg                                                                                                                                    <  2:150507/151‑1 (MQ=255)                 agtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGCGCCCGGGGGGGTTTTCTTTTTCCCAGTg                                                                                                                                    >  2:950122/1‑150 (MQ=255)                 agtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTga                                                                                                                                   <  2:964021/151‑1 (MQ=255)                  gtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTgag                                                                                                                                  <  1:63154/151‑1 (MQ=255)                    gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGAc                                                                                                                                <  2:470541/151‑1 (MQ=255)                     agCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAGTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACg                                                                                                                               <  2:183519/151‑1 (MQ=255)                     agCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTTCCAGTg                                                                                                                                    >  2:688565/1‑146 (MQ=255)                          aCTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTTGGA‑TTTTGCGGCGGGGGGGGTTTTTTTTTTAGCAggg                                                                                                                                    >  2:862622/1‑139 (MQ=255)                             cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACa                                                                                                                        <  1:729257/150‑1 (MQ=255)                              acaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTACCAGTGAGAAGGGCAACa                                                                                                                        >  2:645099/1‑149 (MQ=255)                              acaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTGGGGGCCAGGGGGGGTTTTTTTTTTACCCGTGAGACGGGCaa                                                                                                                          >  2:693506/1‑147 (MQ=255)                               caTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTTGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCt                                                                                                                     <  2:919709/151‑1 (MQ=255)                               caTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGTA‑TTTGGGGGCAGGGGGGGTTTTTTTTTTACCCGTGAGACggg                                                                                                                             >  1:344483/1‑143 (MQ=255)                                         gttgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTGGGGGCCCGGGGGGGTTTTTTTTTTTCCCGCGGGACCGGCCACAGcccgttggc                                                                                                              >  2:80514/1‑140 (MQ=255)                                           tgcgCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTTACAGTGAGACGGGGAAAAACTGA‑TTTcccct                                                                                                          >  2:877802/1‑147 (MQ=255)                                              gcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCg                                                                                                     <  1:422066/151‑1 (MQ=255)                                                 cACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTACCCGTGAGACGGGGAACAGCCGA‑TTGCCCTTCCCCGcc                                                                                                   >  2:651692/1‑150 (MQ=255)                                                         gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTTGTA‑ATTGGGGCCCGGGGGGGTTTTCTTTTTACCACTGAGACGGGGCACAGCTGA‑TTGCCCTTCACCGCCTCGccccg                                                                                           >  2:123785/1‑148 (MQ=255)                                                         gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTTGGGA‑GTGGGCGCCCAGGGGGGTTTTCTTTTTTCCCGCGGGACGGGGAACAACCGA‑TTTCCCCTCTCCCCCCGGcccc                                                                                            >  1:608289/1‑148 (MQ=255)                                                          cTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGGGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCCGGccc                                                                                             >  2:78782/1‑147 (MQ=255)                                                              ccAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTTGTA‑TTTGGGGCCCGGGGGGGTTTTCTTTTTAACAGTGAGGCGGGGAAAAACTGATTTGCCCTTTA‑CGCCCGGCCCTggggg                                                                                       >  1:241261/1‑145 (MQ=255)                                                               cAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTAACAGTGAGACGGGCAACAGGCGA‑TTTCCCTTCACCCCCCGGGCCTGagagag                                                                                     >  2:117900/1‑150 (MQ=255)                                                                aGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCGTA‑TTTGGCGCCAGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCCCCCCCTGGGCCTGAGAGGGtt                                                                                   >  1:194912/1‑151 (MQ=255)                                                                   cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGTA‑TTTGGGGCCAGGGGGGGTTTTCTTTTTACCAATGAGACGGGGAACAGCTGG‑TTGCCCTTTCCCGCCTGGGCCCGggagaggt                                                                                   >  1:835653/1‑146 (MQ=255)                                                                       aaaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGTGTA‑TTTGGGGCCCGGGGGGGTTTTTTTTTTTCCAAGGAGACGGGCAACAGCTGA‑TTGCCCCTTACCCCCTGGCCCCGGGAGAGGTgca‑g                                                                              >  1:658432/1‑148 (MQ=255)                                                                        aaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCA‑GCAAg                                                                          <  1:446639/151‑1 (MQ=255)                                                                         aCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTGGGGGCCCGGGGGGGTTTTTTTTTTACCCGTGAGAAGGGGAACAGCTGATTTGCCCTT‑ACCCCCTGGCCCTGGGAGAGGTGCA‑GCaa                                                                           >  2:88522/1‑149 (MQ=255)                                                                                 tGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTGGGCGCCCGGGGGGGTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCCCCTGGCCCTGGGGGAGTTGCA‑GCAAACGGGCCAc                                                                  >  2:164533/1‑150 (MQ=255)                                                                                  gCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCCGA‑TTGGGCGCCAGGGGGGGTTTTCTTTTTACCCGTGAGAAGGGCAACAGCTGA‑TTGCCCCTTACCCCCCGGCCCCGTGAGAGTTGGC‑GCCAGCGGGcc                                                                    >  2:708890/1‑147 (MQ=255)                                                                                         gCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGTTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCTGGCCCTGAAAGAGTTGCA‑GCAAGCGGTCCACGCTGGttt                                                          <  2:59457/150‑1 (MQ=255)                                                                                          cATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCA‑GCAAGCGGTCCACGCTGGTTTGc                                                        <  2:831669/151‑1 (MQ=255)                                                                                              aaTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCA‑GCAAGCGGTCCACGCTGGTTTGCCCca                                                    <  2:244136/151‑1 (MQ=255)                                                                                                     cGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCA‑GCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCg                                             <  1:677246/151‑1 (MQ=255)                                                                                                        ccAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGCGCCCGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTTACCCCCCGGCCCCGAGAGAGGTTCA‑GCAAGCGGGCCCCGCTGGTTTGGCCCCGCCGGCg                                             >  2:299354/1‑148 (MQ=255)                                                                                                          aaCGCGCGGGGAGAGGCGGTTTTCGGA‑TTTGGGGCCAGGGGGGGTTTTTTTTTTACCCGGGAGGAGGGGAACAGCTGA‑TTTCCCCTCCCCCCCCGGCCCCGAGAGAGTTGCAGGCAAG‑GGGCCACGCTGGTTTGCCCCCGCAGGCGaaaa                                         >  1:929738/1‑150 (MQ=255)                                                                                                          aaCGCGCGGGGAGAGGCGGTTTGCGGA‑TTTGGGGCCAGGGGGGGTTTTTTTTTTTCCACGGAGACGGGCAACAGCCGA‑TTGCCCTTCACCCCCCGGCCCCGAGAGAGGTGCA‑GCCAGCGGGCCACGCTGGGTTTCCCCCGCCGGGGaaaa                                         >  1:401617/1‑150 (MQ=255)                                                                                                                 ggggAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCA‑GCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGttt                                 <  2:533373/151‑1 (MQ=255)                                                                                                                    gagaGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCA‑GCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAt                              <  1:117898/151‑1 (MQ=255)                                                                                                                          cGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCA‑GCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATggtggt                        <  2:4398/151‑1 (MQ=255)                                                                                                                           ggTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCGGTGAGGCGGGCAACAGCTGA‑TTGCCCTTCGCCGCCTGGCCCTGAGAGAGTTGCA‑GCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGtt                       <  1:125186/151‑1 (MQ=255)                                                                                                                           ggTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCA‑GCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGtt                       >  2:732479/1‑151 (MQ=255)                                                                                                                            gTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCA‑GCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtg                             <  1:124662/144‑1 (MQ=255)                                                                                                                              ttGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCA‑GCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTaa                     <  2:551008/150‑1 (MQ=255)                                                                                                                                 cGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCA‑GCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc                 >  2:614369/1‑151 (MQ=255)                                                                                                                                 cGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCA‑GCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc                 >  1:290743/1‑151 (MQ=255)                                                                                                                                 cGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCA‑GCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcgg                  >  2:47873/1‑150 (MQ=255)                                                                                                                                   tA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCA‑GCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg               >  1:433990/1‑151 (MQ=255)                                                                                                                                   tA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCA‑GCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGCGGTTAAcggc                 >  1:702770/1‑149 (MQ=255)                                                                                                                                    a‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCA‑GCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg              <  2:631622/151‑1 (MQ=255)                                                                                                                                      ttGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCA‑GCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg               >  2:798490/1‑149 (MQ=255)                                                                                                                                        gggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCA‑GCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGat            <  2:99534/150‑1 (MQ=255)                                                                                                                                         ggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCA‑GCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat          <  1:363705/151‑1 (MQ=255)                                                                                                                                         ggCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCA‑GCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGatat          <  1:950080/151‑1 (MQ=255)                                                                                                                                           cgcCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCA‑GCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa        <  2:24636/151‑1 (MQ=255)                                                                                                                                                 ggtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCA‑GCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAg  <  1:673695/151‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                       CCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCAGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGA‑TTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCA‑GCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAG  >  NC_000913/366381‑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: ‑