Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,604,237	0	C	T	43.5%	85.1 / 31.2	70	S64L (TCG→TTG)	lsrD	autoinducer 2 import system permease protein
Reads supporting (aligned to +/- strand):  ref base C (13/26);  new base T (30/0);  total (44/26)
Fisher's exact test for biased strand distribution p-value = 1.22e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.78e-01

TATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTT  >  NC_000913/1604094‑1604385                                                                                                                                                |                                                                                                                                                     tATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTgg                                                                                                                                                <  1:1340698/150‑1 (MQ=255)   ttGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGGTTTTTCTTTTTGt                                                                                                                                               >  2:1380340/1‑149 (MQ=255)    tGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTc                                                                                                                                             <  2:1035075/150‑1 (MQ=255)     gagaTTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTttgttttgtt                                                                                                                                              >  1:623598/1‑148 (MQ=255)     gagaTTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGtttg                                                                                                                                            >  1:629014/1‑148 (MQ=255)        aTTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGGCGATATTTTGttttggttggcc                                                                                                                                         >  2:594998/1‑141 (MQ=255)               cATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTCGTTTTGTTCGACCATCGGcc                                                                                                                                  >  1:1264246/1‑150 (MQ=255)                aTTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGAGATTGTCAGTGGCGGGGTCGATATTTAGTTTAGTTTGTCCATCAGcc                                                                                                                                  >  1:1404884/1‑149 (MQ=255)                 tttGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCCATATTTTGTTTGGGTccgcccgcagac                                                                                                                                  >  2:394897/1‑137 (MQ=255)                    ggTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTTGACCATCGGcccc                                                                                                                                >  1:1276630/1‑145 (MQ=255)                     gTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGGTATTTTGTTTGGGTCGTCCATCGGcc                                                                                                                                  >  2:552920/1‑144 (MQ=255)                     gTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGTTTTGTTCTACCATCGGGctc                                                                                                                                >  1:1317878/1‑146 (MQ=255)                       gCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTTGACCATCCGCctc                                                                                                                                >  2:143906/1‑144 (MQ=255)                        cAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCATCGGCCCCTGCGcc                                                                                                                          >  1:488089/1‑149 (MQ=255)                            tAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTCTTTGGGTTGACCATCGGCCTCTGCGCcatt                                                                                                                       >  2:21007/1‑148 (MQ=255)                            tAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCcattgc                                                                                                                     <  1:295923/150‑1 (MQ=255)                             aaCCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCcattgca                                                                                                                    <  2:766405/150‑1 (MQ=255)                              aCCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTTGGTATTTTGTTTTGTTTGACCCTCCGCCCCTGCGGccttgca                                                                                                                    >  2:1149768/1‑149 (MQ=255)                              aCCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCcattgcatt                                                                                                                  <  2:269154/151‑1 (MQ=255)                               cccGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATTGATATTTTGTTTTGTTTGGCCATCAGCCTCTTCGCcctttcatt                                                                                                                  >  2:180781/1‑150 (MQ=255)                                ccGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTgg                                                                                                                <  2:1121816/151‑1 (MQ=255)                                 cgcgAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTggg                                                                                                               <  1:1140259/151‑1 (MQ=255)                                  gcgAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGc                                                                                                              <  2:248375/151‑1 (MQ=255)                                   cgAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCg                                                                                                             <  2:438445/151‑1 (MQ=255)                                      aTGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGAAGATATTTTGTTTTGTTCGACCATCCGCCCCTCCGCCCTTTTATTGGGCGGAc                                                                                                          >  2:1369882/1‑151 (MQ=255)                                       tGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTTTTTTGTTTGACCATCCGCCTCTGCGCCATTGCATTTGGCGGAc                                                                                                          >  1:879653/1‑150 (MQ=255)                                        gTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTg                                                                                                        <  2:61212/151‑1 (MQ=255)                                           aGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTTGATATTTTGTTTTGTTTTACCATCGGCCCCTGGGCCCTTTCATTGGGCggatttgtt                                                                                                     >  1:875257/1‑143 (MQ=255)                                           aGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATTGATATTTTGTTTGGTTTGACCATCGGCCTCTGCGCCCTTTCATTGGGCGGACTGttt                                                                                                     >  2:885245/1‑151 (MQ=255)                                           aGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGtt                                                                                                      <  1:288985/150‑1 (MQ=255)                                              tctcAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGTTTGACCCTCCGCCTCTGCGCCCTTGCCTTGGGGGGACTGttt                                                                                                     >  1:1209012/1‑148 (MQ=255)                                                     atGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTGGTTCGACCATCGGCCTCTTCGCCCTTTGATTGGGCGGGCTGTTTTAATgtgg                                                                                             >  1:520107/1‑149 (MQ=255)                                                            tGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATGGGTATTTTGTTTGGTTCGACCATCTGCCTCTTCGCCATTGCATTTGGGGGACTGTTTTAAAATGGggtg                                                                                         >  1:347894/1‑146 (MQ=255)                                                            tGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTTGACCATCCGCCTCTTCGCCATTTCATTGGGCGGACTGTTTTAAAATGGggtggc                                                                                       >  1:301796/1‑146 (MQ=255)                                                               tCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGc                                                                                  <  2:1100387/150‑1 (MQ=255)                                                                cAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCg                                                                                <  2:398214/151‑1 (MQ=255)                                                                  gCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCt                                                                              <  1:1324347/151‑1 (MQ=255)                                                                   cACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCt                                                                              <  2:135612/150‑1 (MQ=255)                                                                      cAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCg                                                                          <  2:416625/151‑1 (MQ=255)                                                                        gTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGAATGGTATTTTGTTTTGTTCGACCATCGGCCTCTTCGCCCTTTGCTTGGGGGGTCTGTTTTAAAGGGGGGTGCCGGGGCCGCTgg                                                                            >  2:677547/1‑147 (MQ=255)                                                                               tATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGGTATTTTGTTTTGTTCGACCATCGGCCCCTGCGCCATTGCATTTGGGGTACTGTTTCCAAGGTGGGTGCCCATGCCGCCGGCGGTACTcc                                                                   >  1:578851/1‑149 (MQ=255)                                                                                         ggcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATAGGTATTTTGTTTTGTTCGACCATCCGCCTCTCCGCCCTTGCATTGGGGGGACTGTTTTTAAGTGGGGGGCCGCGGCCGCTGGGCGTACTCCTTACCtt                                                            >  2:95524/1‑146 (MQ=255)                                                                                          gcattgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCt                                                      <  1:65942/151‑1 (MQ=255)                                                                                              tgTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTTGACCATCGGCCCCTGCGCCATTGCATTGGGCGTACTGTTTTAAAATGGGGGGCCGATGCCGCTGGCGGTACTCCTGACCTTACTGCTc                                                     >  1:482947/1‑148 (MQ=255)                                                                                              agTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATTGATATTTTGTTTTGTTGGACCATCGGGCCCTGCGCCAATGGCTTGGGGGGACTGTTTTGAAGTGGGGGGGCGGTGCCGCTGGGGGTACTTCTGGCCTTTCTGCTCgg                                                   >  2:554653/2‑150 (MQ=255)                                                                                               gTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGgcg                                                 <  1:199956/151‑1 (MQ=255)                                                                                                   cccTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTTTTTTGTTTCACCATCGGCCCCTGCGCCCTTGCATTGGGCGGACTGTTTTAAAGTGGTGTGGCCCTGCCCCTGGGGGTACTCCTTACCTTACTGCTCGGCGCAtt                                             >  2:1327593/1‑151 (MQ=255)                                                                                                    ccTACCGCTAACGATGGTGATTGTCAGTGGCGGGGTCGATATTTTGTTTGGTTCGACCATCGGCCTCTGCGCCCTTTCATTGGGCGTACTGTTTTAAAGTGGGGTGCCGATGCCGCCGGCGATACTCCTGACCTTACTGCTCGGCGCAtt                                             >  1:300797/1‑150 (MQ=255)                                                                                                      tACCGCTAACGATGGTGATTGTCAGTGGCGGGATTGATATTTTGTTTGGTTTGTCCATCGGCCTCCTCGCCATTGTATTGGGGGGTCTGTTTCAAAAAGGGGTGGCGACGCCCCTGGGGAAACTCCTGGCCTTACTACTCGGCGCaat                                             >  1:1335686/1‑146 (MQ=255)                                                                                                      tACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTTTTTGGTTTGACCATCGGCCTCTGCGCCCTTTGCTTGTGCGTTCTGTTTTAAAGTGGGGTGCCGGTGGCGCTGGCGGTACCCCTGGCCCTACTGCCCGGGGCATtgt                                           >  2:80243/1‑150 (MQ=255)                                                                                                      tACCGCTAACGATGGTGATTGTCAGTGGCGGGAACGATATTTTGGTTGGTTTGCCCACAAGCCTCTGCGCCCATGCCTTGGGGGTACTTTTTTCAAGTGGGGTGGCGAAGCCGCTTGCGGTAATCCTGCCCTTACTtctctgggcct                                              >  1:769471/1‑140 (MQ=255)                                                                                                       aCCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATtgtg                                          <  2:270840/150‑1 (MQ=255)                                                                                                         cGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGGACTGTTTTAAAATGGTGTGCCGATGCCGCTGGCGGTACTCCTTACCTTACTGCTCGGCGCATTGTGCgg                                       >  1:320463/1‑151 (MQ=255)                                                                                                          gCTAACTATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGTCCTTACTGCTCGGCGCATTGTGCgg                                       <  2:1049049/150‑1 (MQ=255)                                                                                                          gCTAACGATGGTGATTGTCAGTGGCGGGATCGAAATTTTGTTTTGTTTGACCATCGGCCTCTTCGCCCTTGCATTGGGCGGACTGGTTTAAAGTGGTGTGGCGACGGCCCTGGCGATACTCCCGCCCTTACTGCTCCGCGCATTGTGCggg                                      >  1:363353/1‑151 (MQ=255)                                                                                                             aaCGATGGTGATTGTCAGTGGCGGGATCGATATTTTGTTTTGTTCGACCAACGGCCTCCGCGCCATTGCATTTGGGGTTCTGTTTTAAAGTGGGGTGCCGGTGCCGCGGGCGGTACTCCTGACCTTTCTGCTCGGCGCATTGTGCgg                                       >  2:657968/1‑147 (MQ=255)                                                                                                              aCGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGa                                  <  1:885266/151‑1 (MQ=255)                                                                                                                gATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATc                                <  2:623585/151‑1 (MQ=255)                                                                                                                               tGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGattaatt                 >  2:1216789/1‑151 (MQ=255)                                                                                                                               tGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGattaatt                 >  2:739399/1‑151 (MQ=255)                                                                                                                               tGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGattaatt                 >  1:226343/1‑151 (MQ=255)                                                                                                                                  cGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATc              <  2:309491/151‑1 (MQ=255)                                                                                                                                   gggATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCt             >  2:843252/1‑151 (MQ=255)                                                                                                                                     gATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCtat           <  2:764521/151‑1 (MQ=255)                                                                                                                                        cGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATAcc        <  2:1276602/151‑1 (MQ=255)                                                                                                                                          atatTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCa       >  2:674947/1‑150 (MQ=255)                                                                                                                                           tatTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCaa      >  2:315765/1‑150 (MQ=255)                                                                                                                                           tatTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCaa      >  1:476416/1‑150 (MQ=255)                                                                                                                                             tttCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGt   >  2:537023/1‑151 (MQ=255)                                                                                                                                               tCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGtt  >  1:1037357/1‑150 (MQ=255)                                                                                                                                               tCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGtt  <  2:312937/150‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                     TATTGAGATTGTCGCATTTGGTGCAATTAACCCGCGAATGTTAGATCTCAATATGTTGCTGTTCAGCACCAGTGACTTTATCTGCATTGGCATTGTCGCCCTACCGCTAACGATGGTGATTGTCAGTGGCGGGATCGATATTTCGTTTGGTTCGACCATCGGCCTCTGCGCCATTGCATTGGGCGTACTGTTTCAAAGTGGTGTGCCGATGCCGCTGGCGATACTCCTGACCTTACTGCTCGGCGCATTGTGCGGGCTGATCAACGCCGGATTAATTATCTATACCAAAGTT  >  NC_000913/1604094‑1604385

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 13 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑