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

AAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGC  >  NC_000913/366378‑366675                                                                                                                                                       |                                                                                                                                                      aaGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTC‑GGAAACCTGTCGTGCCAGCTGCCTTAATGAATCGGCCAACGCGGGGGGAGAGGCGGGTTTCGGA‑TTGGGGGCC‑Ggggggg                                                                                                                                                     >  2:512492/1‑150 (MQ=255)  aGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGggtggt                                                                                                                                                    <  2:111689/151‑1 (MQ=255)         ggTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGAGTA‑TTGGGCGCC‑AGGGTGGTttttctt                                                                                                                                             <  2:662971/151‑1 (MQ=255)            gCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTttttcttttc                                                                                                                                          <  1:650508/151‑1 (MQ=255)             ccTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGGGGTttttctttt                                                                                                                                           >  1:508811/1‑149 (MQ=255)                aaTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCAcc                                                                                                                                       <  1:286628/150‑1 (MQ=255)                 aTGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAg                                                                                                                                     <  2:288369/151‑1 (MQ=255)                   gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGGGCC‑CGGGGGGGTTTTTTTTTTACCAGTg                                                                                                                                   >  1:615550/1‑151 (MQ=255)                   gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTC‑TTTTGGGGC‑CGGGGGGTTTTTCTTTTCTccccg                                                                                                                                    >  2:1129495/1‑147 (MQ=255)                   gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTg                                                                                                                                   <  1:145948/151‑1 (MQ=255)                   gagtgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTg                                                                                                                                   <  2:320847/151‑1 (MQ=255)                      tgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGGGCC‑CGGGGGGGTTTTTTTTTTACCAGTg                                                                                                                                   >  1:1108142/1‑148 (MQ=255)                      tgagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGGGCC‑CGGGGGGGTTTTTTTTTTACCAGt                                                                                                                                    >  1:356431/1‑147 (MQ=255)                       gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGAc                                                                                                                               <  2:571449/151‑1 (MQ=255)                       gagCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGAc                                                                                                                               <  1:716345/151‑1 (MQ=255)                        agCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACg                                                                                                                              <  2:64000/151‑1 (MQ=255)                        agCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACg                                                                                                                              <  1:965077/151‑1 (MQ=255)                           tAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTA‑TTTGGCGGC‑AGGGGGGGTTTTTTTTTTACCAATGAGACggg                                                                                                                            >  1:826618/1‑150 (MQ=255)                                cacaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGGA‑TTGGGGGCC‑CGGGGGGGTTTTCTTTTTCTCCGTGAGAAGGGCaa                                                                                                                         >  2:917026/1‑148 (MQ=255)                                 acaTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCCGCGGGGGGTTTTTTTTTTC‑CCAGTgcggcg                                                                                                                              >  2:586523/1‑137 (MQ=255)                                   aTTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGCGGC‑AGGGGGGGTTTTCTTTTTACCAGTGAGACGGGCAACa                                                                                                                       >  1:1085739/1‑147 (MQ=255)                                      aaTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTATTTGGGGGCC‑GGGGGGGTTTTT‑TTTTTACCAGGGAGACGGGGAACAGCTGAtt                                                                                                                >  1:336218/1‑151 (MQ=255)                                       aTTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTg                                                                                                               <  1:679808/151‑1 (MQ=255)                                               gcgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTTCGTA‑TTGGGGGGC‑CGGGGGGGTTTTTTTTTTACCAGGGAGACGGGGAACAGCTGATTGCCCTTCAc                                                                                                       >  2:480253/1‑151 (MQ=255)                                                cgcTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGGGCC‑AGGGGGGTTTTTCTTTTCACCCGTGAGACGGGCAACAGCTGATTGCCCTTCAcc                                                                                                      >  1:1183128/1‑151 (MQ=255)                                                   tCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGGTTGCGTA‑TTGGGGGCC‑AGGGGGGGTTTTTTTTTTTCCAGGGGGGCGGGGAAAAGCCGATTGCCCTTCACcccc                                                                                                   >  2:94133/1‑151 (MQ=255)                                                     aCTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTGCGCT‑TTAGGCGCC‑AGGGGGGTTTTTTTTTTTAACCCGGTAACGGGCAACAACCGTTTGGCCttttccccc                                                                                                   >  2:986016/1‑142 (MQ=255)                                                      cTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTg                                                                                                 <  2:137248/150‑1 (MQ=255)                                                         cccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTGGGGGCC‑CGGGGGGGTTTTTTTTTTCCCAGTGAGACCGGCAACAGCTGATTGCCCTTCCCCGCCTGGgcc                                                                                             >  1:396474/1‑151 (MQ=255)                                                          ccGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGT‑TTTGGGGCC‑CGGGGGGGTTTTTTTTTTACCCGTGAGACGGGCGACAGCCGATTGCCCCTCCCCACCCCGccc                                                                                             >  2:1158559/1‑150 (MQ=255)                                                            gCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTga                                                                                          <  2:829539/151‑1 (MQ=255)                                                               ttCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgag                                                                                         <  2:681617/149‑1 (MQ=255)                                                                tCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGGGTA‑TTGGGGGGC‑AGGGGGGGTTTTTTTTTTACCAGTGAGAAGGGGAAAAGCTGATTGCCCCTCACCGCCCGGCCCTgaggg                                                                                       >  1:919886/1‑148 (MQ=255)                                                                tCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTgagag                                                                                       <  2:207901/150‑1 (MQ=255)                                                                  cAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGTT‑TTTGGGGCC‑CGGGGGGGTTTTTTTTTTTCCCGGGGGACGGGCAACAGCCGGTTTCCCTTTCCCGCCCGGcccctggggag                                                                                     >  2:627375/1‑142 (MQ=255)                                                                   aGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTTGGGGCC‑AGGGGGGGTTTTTTTTTTACCAGGGAGACGGGGCACAGCTGATTGCCCTTCACCCCCCGGCCCTGagagaggt                                                                                   >  2:789811/1‑149 (MQ=255)                                                                      cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGGGCC‑AGGGGGGGTTTTTTTTTTACCCGTGAGACGGGCAACAGCTGATTGCCCTTCACCCCCTGGCCCTgagaga                                                                                      >  2:297957/1‑145 (MQ=255)                                                                      cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgca                                                                                <  1:1131594/151‑1 (MQ=255)                                                                      cGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGA‑TTTGGGGCC‑AGGGGGGGTTTTTTTTTCACCACTGAGACGGGCAACAACCGATTTCCCCTCACCGCCCGGCCCTGggagaggt                                                                                   >  2:119893/1‑146 (MQ=255)                                                                       gggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGCGCC‑CGGGGGGGTTTTTTTTTTACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCCGGCCCCGAGAGGGtt                                                                                   >  1:288478/1‑147 (MQ=255)                                                                        ggAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGGA‑TTTGGGGCC‑AGGGGGGGTTTTTTTTTTACCAGACAGCCGGGCCACACGTGTTAGCCCTTCCCCATCCGGCCCCGGGAGAGTTgcag                                                                               >  2:806317/1‑150 (MQ=255)                                                                         gAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTgcagca                                                                             <  1:287089/151‑1 (MQ=255)                                                                          aaaCCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCaa                                                                            <  1:921894/151‑1 (MQ=255)                                                                                gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGt                                                                       <  2:796307/150‑1 (MQ=255)                                                                                gTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGGT‑TTTGGGGCC‑CTGGGGGTTTTTCTTTTCACCCGAGAAACGGGCAACAACTGGCTGCCCCTCTCCCCCCGGCCCTGAGAGAGGTTCAGCCAGCggg                                                                       >  1:727107/1‑149 (MQ=255)                                                                                     gCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCACCAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCAc                                                                   <  2:83090/149‑1 (MQ=255)                                                                                     gCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGCTTGCGGA‑TTTGGGGCC‑CGGGGGGGTTTTCTTTTTAACAGTGAGACGGGGCACAGCCGATTGCCCTTCCCCCCCTGGCCCCGGGGGAGTTGCAGCAAACGGGGCACg                                                                  >  2:590977/1‑150 (MQ=255)                                                                                      ccAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCt                                                                <  2:845714/151‑1 (MQ=255)                                                                                       cAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGATTGGGGA‑TTTGGCGCC‑CGGGGGGGTTTTTTTTTTACCAGTGAGACGGGGAACAGATGGTTGGCCCTCACCCCCCGGGCCTGCGAGGGGTGCAGCAAACGGGCCAccccg                                                               >  2:747851/1‑147 (MQ=255)                                                                                         gCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTTGGCGCC‑AGGGGGGTTTTTTTTTTTTCCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCCCCCGGCCCTGAGAGAGGTGCAGCAAACGGGGCACGc                                                                 >  2:47857/1‑147 (MQ=255)                                                                                           tGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTGGTA‑TTTGGGGCC‑CGGGGGGGTTTTTTTTTTACCAGGGAGACGGGGAACAGCTGGTTGCCCCTCCCCCCCCGGCCCTGGGAGAGGTGCAGCAAGCGGTCCCCCCCGGttt                                                           >  2:652417/1‑151 (MQ=255)                                                                                                 aaTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCca                                                     <  1:12254/151‑1 (MQ=255)                                                                                                  aTGAATCGGCCAACGCGCGGGGAGAGGCGGTTTTCGTA‑TTTGGCGGC‑AGGGTGGTTTTTCTTTTTCCCAATGAGACGGGCAACAGCTGATTGCCCTTCACCCCCTGGCCCTGAGAGAGTTGCAGCAAGAGGTCCACGCTGGGTTGcccc                                                      >  1:382846/1‑149 (MQ=255)                                                                                                     aaTCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCcagcag                                                 <  2:783297/151‑1 (MQ=255)                                                                                                        cGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGc                                               <  2:94661/150‑1 (MQ=255)                                                                                                        cGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGc                                               <  2:349451/150‑1 (MQ=255)                                                                                                        cGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGc                                               <  2:865880/150‑1 (MQ=255)                                                                                                         ggCCAACGCGCGGGGAGAGGCGGTTTGGGTA‑TTGGGGGCC‑AGGGGGGGTTTTCTTTTCACCAGGGAGACGGGCAACAGCTGATTGCCCTTCACCCCCTGGCCCTGAGAGAGTTGCCGCAAGCGGTCCACGCCGGTTTTCCCCCGCCgg                                                >  2:173485/1‑148 (MQ=255)                                                                                                              aCGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑CGGGGGGGTTTTCTTTTTACCCGGGAGACGGGCAACAGCTGGTTGCCCCTCACCGCCTGGCCCTGGGGGGGGTGCAGCAAGCGGGCCACGCCGGGTTGCCCCCGCAGGGGaaaa                                          >  2:413456/1‑149 (MQ=255)                                                                                                                 cgcgGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAAt                                         <  2:259340/147‑1 (MQ=255)                                                                                                                   cgGGGAGAGGCGGTTTTCGTA‑TTTGGCGCC‑CGGGGGGGTTTTTTTTTCACCAGTGAGACGGGGCACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGGGAGTTGCAGCAAGCGGGCCCCGCCGGGTTGCCCCAGCAGGCGGAAATCCTGtt                                   >  1:518792/1‑151 (MQ=255)                                                                                                                         gagGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACTGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtgg                             <  1:139984/151‑1 (MQ=255)                                                                                                                          agGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggt                            <  1:693268/151‑1 (MQ=255)                                                                                                                          agGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtgg                             >  2:169352/1‑150 (MQ=255)                                                                                                                           ggCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGAtggtg                           >  1:22049/1‑151 (MQ=255)                                                                                                                              ggTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATggtggt                         <  1:717183/150‑1 (MQ=255)                                                                                                                              ggTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGtt                        <  1:383198/151‑1 (MQ=255)                                                                                                                               gTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTa                       <  2:479167/151‑1 (MQ=255)                                                                                                                                    cGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcgg                   >  2:40614/1‑150 (MQ=255)                                                                                                                                     gTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggc                  <  2:528928/150‑1 (MQ=255)                                                                                                                                      tA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAAcggcgg                >  2:928924/1‑151 (MQ=255)                                                                                                                                         ttGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGTggg               <  2:805639/150‑1 (MQ=255)                                                                                                                                           gggCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCggg               >  2:640254/1‑148 (MQ=255)                                                                                                                                               gcC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAAc        <  1:820640/151‑1 (MQ=255)                                                                                                                                               gcC‑AGGGGGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGAAATaa         >  1:1027580/1‑150 (MQ=255)                                                                                                                                                   aGGGTGGTTTTTCTTTTCCCCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACAt      <  1:61236/150‑1 (MQ=255)                                                                                                                                                   aGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATaa         >  1:719397/1‑147 (MQ=255)                                                                                                                                                   aGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACa       >  1:29729/1‑149 (MQ=255)                                                                                                                                                    gggTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGACCTAACa       >  2:512693/1‑148 (MQ=255)                                                                                                                                                     ggtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAg   <  2:908225/151‑1 (MQ=255)                                                                                                                                                      gtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGc  <  1:283909/151‑1 (MQ=255)                                                                                                                                                      gtggtTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGc  <  2:295829/151‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                      AAGCCTGGGGTGCCTAATGAGTGAGCTAACTCACATTAATTGCGTTGCGCTCACTGCCCGCTTTCCAGTCGGGAAACCTGTCGTGCCAGCTGCATTAATGAATCGGCCAACGCGCGGGGAGAGGCGGTTTGCGTA‑TTGGGCGCC‑AGGGTGGTTTTTCTTTTCACCAGTGAGACGGGCAACAGCTGATTGCCCTTCACCGCCTGGCCCTGAGAGAGTTGCAGCAAGCGGTCCACGCTGGTTTGCCCCAGCAGGCGAAAATCCTGTTTGATGGTGGTTAACGGCGGGATATAACATGAGC  >  NC_000913/366378‑366675

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