Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,555,109	0	T	A	43.1%	62.4 / 30.6	58	N178K (AAT→AAA)	uxuR	fructuronate‑inducible hexuronate regulon transcriptional repressor; autorepressor
Reads supporting (aligned to +/- strand):  ref base T (11/22);  new base A (25/0);  total (36/22)
Fisher's exact test for biased strand distribution p-value = 3.95e-08
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATTGCAACTGGAAGAGCGTGAACTGGCTTCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAATCCAATGT‑GGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCGCAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGCCCGAGCGGCGAAGCTGGCAATGTGGCAGCATCTGGAAA  >  NC_000913/4554959‑4555251                                                                                                                                                        |                                                                                                                                                aTTGCAACTGGAAGAGCGTGAACTGGCTTCTAGTGCGCCGGGCAGCAGCGAAAGCGTTGACATGCAGTTCCCTCTCGCTATTGCCTAAGAAACGCATAACAGCATGCTGGTGGCGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACCAt                                                                                                                                                 <  1:143451/151‑1 (MQ=255)    gCAACTGGAAGAGCGTGAACTGGCTTCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAATcc                                                                                                                                               <  2:363196/150‑1 (MQ=255)    gCAACTGGAAGAGCGTGAACTGGCTTCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAATCCa                                                                                                                                              <  1:1197974/151‑1 (MQ=255)     cAACTGGAAGAGCGTGAACTGGCTTCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAATCCaa                                                                                                                                             <  1:779960/151‑1 (MQ=255)       aCTGGAAGAGCGTGAACTGGCTTCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAATCCaa                                                                                                                                             >  2:192535/1‑149 (MQ=255)                  gTGAACTGGCTTCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAATCCAATGT‑GGATCCAGtt                                                                                                                               <  1:1548/151‑1 (MQ=255)                   tGAACTGGCTTCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGGAAAAAATCCAAAGT‑GGATcccggt                                                                                                                               >  1:1046875/1‑146 (MQ=255)                    gAACTGGCTTCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAATCCAATGT‑GGATCCAGTTGc                                                                                                                             <  2:1108313/151‑1 (MQ=255)                    gAACTGGCTTCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAAAAATCCAATGG‑GGGTTCAGt                                                                                                                                >  2:81952/1‑148 (MQ=255)                     aaCTGGCTTCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAAACCAATGG‑GGATCCCGTTGc                                                                                                                             >  2:1150825/1‑150 (MQ=255)                       cTGGCTTCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGCGGGGAAAAAAATCCAA‑GT‑GGATcc                                                                                                                                   >  1:600179/1‑142 (MQ=255)                        tGGCTTCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAAAAAACCAATGT‑GGGTCCAGGTGgac                                                                                                                           >  1:74099/1‑149 (MQ=255)                          gCTTCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGGAAAAAAAACAAAGT‑GGATCCAATTGCACAGcc                                                                                                                       >  1:438630/1‑151 (MQ=255)                           cTTCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAAACCCATGT‑GGGGCCCGTTGCCCAGcc                                                                                                                       >  2:1188627/1‑150 (MQ=255)                            ttCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAAACCAATGT‑GGATCCCGGTGCACAGcc                                                                                                                       >  1:158206/1‑149 (MQ=255)                             tCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAAACCAATGT‑GGGTCCAGTTTCACAGcc                                                                                                                       >  1:859636/1‑148 (MQ=255)                               cAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAATCCAATGT‑GGATCCAGTTGCACAGCCATCTg                                                                                                                  <  1:747772/151‑1 (MQ=255)                               cAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAATCCAATGT‑GGATCCAGTTGCACAGCCATCTg                                                                                                                  <  1:1020021/151‑1 (MQ=255)                                        gggCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAATCCAATGT‑GGATCCAGTTGCACAGCCTTCTGGATGACAGc                                                                                                         <  1:867938/151‑1 (MQ=255)                                        gggCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAAAAAACCAATGT‑GGGGCCCGGTGGACAGCCATCTGGGTGACAg                                                                                                          >  2:157512/1‑150 (MQ=255)                                         ggCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAAAAATCCCAAGT‑GGAGCCCGTTGCACAGCCATCTTGGTGACAGcc                                                                                                        >  1:417986/1‑151 (MQ=255)                                            agcagcGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGCAAAAAATTCAATTT‑GGGTCCCGTTGCACAGCCATCCTGGTGACAGcc                                                                                                        >  2:435201/1‑148 (MQ=255)                                             gcagcGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAATCCAATGT‑GGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTa                                                                                                    <  1:798999/151‑1 (MQ=255)                                               agcGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAAAAAACCCATGT‑GGGTCCCGTTGCACAGCCCTCTGGGTGACAGCCTGTa                                                                                                    >  2:590899/1‑149 (MQ=255)                                                gcGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAATCCAATGT‑GGGTCCAGTTGCACAGCCATCTTGATGACAGCCTGTAt                                                                                                   >  2:350228/1‑149 (MQ=255)                                                 cGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGAGAAACAAACCAATAT‑GGTTACTCTTGTCCATCCATCTGGATTACAAcccggatt                                                                                                  >  2:822181/1‑143 (MQ=255)                                                        ggTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAATCCAATGT‑GGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCGCAAAGAGt                                                                                         <  1:980985/151‑1 (MQ=255)                                                        ggTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAAAAAACCAATGT‑GGGTCCAGTTGGACCGCCATCTGGGTGACAGCCCGTATCGCAAaga                                                                                           >  1:72083/1‑149 (MQ=255)                                                              aTGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAATCCAATGT‑GGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCGCAAAGAGTGGTTgg                                                                                   <  2:615249/151‑1 (MQ=255)                                                                gCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAATCCAATGT‑GGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCGCAAAGAGTGGTTGGGc                                                                                 <  1:429580/151‑1 (MQ=255)                                                                       cATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAATCCAATGT‑GGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCGCAAAGAGTGGTTGGGCGATcaca                                                                          <  2:407892/151‑1 (MQ=255)                                                                       cATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAAAAAACCCATGT‑GGGTCCCGTTTCACAGCCATCTGGATGACAGCCCGTATCGCAAAGAGGGGGTGGGCGGTc                                                                             >  2:55282/1‑148 (MQ=255)                                                                         tctcGCTATTGCCGAAGCAACGCATAACAGCATGCTGGCGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAAAAAAACAATGT‑GGGTCCAGTTGCACAGCCATCTGGGTGACAGCCCGTATCGCAAAAAGGGGGTGGGCGGTCACaa                                                                         >  1:859224/1‑150 (MQ=255)                                                                          ctcGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAAAAAACCAATGT‑GGGTCCAGTTGCACAGCCACCTGGATGACAGCCTGTTTCGCAAAGAGGGGGTGGGCGATCACaaa                                                                        >  2:224547/1‑150 (MQ=255)                                                                            cGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAAAAAACCAAAGT‑GGGTCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCGCAAAGAGGGGGTGGGCGGTCACAAACa                                                                      >  2:92771/1‑150 (MQ=255)                                                                               tATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAAAAAACCCAAGT‑AGGTCCAGTTGGCCACCCCACTGGGATAAAGCCCGGCTCGCAAAAAGAGGTTggggggccccaaaaacaa                                                                   >  1:450377/1‑135 (MQ=255)                                                                                  tGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAAACCAATGTAGGAT‑CAGATGGACAGGCAGCAGGATGACAGACCGTAAACCAAAGAGTGGGTGGGAGATCACAAACCGATcc                                                                 >  1:451266/1‑149 (MQ=255)                                                                                    ccGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAAACCCATGT‑GGATCCAGTTTCACAGCCATCTGGGTGACAGCCTGTATCGCAAAGAGTGGGTGGGCGGTCACAAACAGATCCTcgc                                                             >  2:389424/1‑151 (MQ=255)                                                                                      gAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAAAAAACCAATGT‑GGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCGCAAAGAGTGGTTGGGCGGTCACAAAAAGATCCTcgcc                                                            >  2:278705/1‑150 (MQ=255)                                                                                       aaGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAAAAATCCAATGT‑GGGTTCAGGTGCACAGCCCTCTGGGTGACAGCCTGTATTGCAAAGAGGGGTTGGGGGGTCACAAACAAATCCTcgccgc                                                          >  1:701599/1‑151 (MQ=255)                                                                                           aaCGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAATCCAATGT‑GGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCGCAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGtt                                                       <  1:586799/150‑1 (MQ=255)                                                                                            aCGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAATCCAATGT‑GGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCGCAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTaa                                                     <  2:590475/151‑1 (MQ=255)                                                                                              gCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAAAAAACAAAAGC‑GGGTGCCATTGTAAAACCATCTGGGAGAGAAGCTGTTTTGCAAAGGGAGGTTTGGCGATTTCAAAAAAAGTCTCGCCGGGtt                                                       >  2:843079/1‑147 (MQ=255)                                                                                                aTAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAAACCCATGT‑GGGTCCAGTTGCACAGCCATCTGGGTGACAGCCTGTATCGCCAAGAGGGGGTGGGCGATCACAAACAGATCCTCCCCGCGTTAATCa                                                  >  1:628/1‑150 (MQ=255)                                                                                                       cATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAATCCAATGT‑GGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCGCAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGAt                                          <  1:157510/151‑1 (MQ=255)                                                                                                        aTGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAATCCAATGT‑GGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCGCAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATg                                         <  1:524327/151‑1 (MQ=255)                                                                                                           cTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGGAAACAAACCAATGA‑TGGTCCAAGTGGACCAACATCTGGAAGAAAAAACATATCTGACAGAGTGAGTCGGCGGTCACAAAACGATACCCGCCCCGGTAATAaaaaaa                                             >  1:1180126/1‑144 (MQ=255)                                                                                                               tggAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAAAAAACCAATGT‑GGGTCCAGTTGGCCAGCCCTCTTGAGGGCAACCCGTATCGCCAAAAGCGGGTTGGGGATCACAAACAGATCCCCGCCGCGGTTATCAAAAAAAAAGccc                                      >  1:199786/1‑147 (MQ=255)                                                                                                                 gAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAATCCAATGT‑GGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCGCAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGCCCGAgcgg                                <  2:746668/151‑1 (MQ=255)                                                                                                                   gCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGGAAAAAAACCAATGT‑GGATGCAGTTTTCCAGCCATCTGGAAGTTAGCCTGTAACGCAAAGGGAGGGTTGGGGCACACAAACAGAACCCCCCCCCGTTTATCaaaaaaaaaggcccag                                   >  2:791837/1‑136 (MQ=255)                                                                                                                   gCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAAAAATCCAATGT‑GGGTCCAGTTTCACAGCCACCTGGGTGACAGACTGTTTCGCCAAGAGGGGGTGGGGGGTCACAAACAAATCCCCCCCGCGTTAATCAAAAAAAAAGCCCGAgcggcg                              >  2:387562/1‑151 (MQ=255)                                                                                                                      gTTCCGTCAGTCCTGGCAGTGGC‑GGGGAAAAAAACCAATGT‑GGGTCCAGGTGCCCAGCCCTCTGGATGAGAGCCTGTATCGCCAAAAGGGGTTGGGGGGTCACAAAAAGATCCTCCCCGCGTTAATCAAAAAAAATGTCCGAgcggcg                              >  2:586275/1‑148 (MQ=255)                                                                                                                        tCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAATCCAATGT‑GGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCGCAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGCCCGAGCGGCGAAGc                          <  2:1071333/150‑1 (MQ=255)                                                                                                                         ccGTCAGTCCTGGCAGTGGC‑GGGAAAACAATCCAATGT‑GGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCGCAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGCCCGAGCGGCGAAGCTg                        <  2:395148/151‑1 (MQ=255)                                                                                                                            tCAGTCCTGGCAGTGGC‑GGGAAAAAAAACCCATGT‑GGGTCCCGGTGCACAGCCCTCTGGGTGGCAGCCTGTATCGCAAAAAGGGGGTGGGGGGTCACAAACAGATCCTCGCCGCGTTTATCAAAAAAAAAGCCCCAGCGGCGAAGATgg                       >  1:224754/1‑149 (MQ=255)                                                                                                                                 ccTGGCAGTGGC‑GGGAAAACAATCCAATGT‑GGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCGCAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGCCCGAGCGGCGAAGCTGGCAAtgtg                <  1:46729/151‑1 (MQ=255)                                                                                                                                                gAAAACAATCCAATGT‑GGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCGCAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGCCCGAGCGGCGAAGCTGGCAATGTGGCAGCATCTGGaaa  <  2:266297/151‑1 (MQ=255)                                                                                                                                                gAAAACAATCCAATGT‑GGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCGCAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGCCCGAGCGGCGAAGCTGGCAATGTGGCAGCATCTGGaa   >  1:1155942/1‑150 (MQ=255)                                                                                                                                                gAAAACAATCCAATGT‑GGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCGCAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGCCCGAGCGGCGAAGCTGGCAATGTGGCAGCATCTGGaa   >  1:647623/1‑150 (MQ=255)                                                                                                                                                        |                                                                                                                                                ATTGCAACTGGAAGAGCGTGAACTGGCTTCCAGTGCGCCGGGCAGCAGCGAAAGCGGTGACATGCAGTTCCATCTCGCTATTGCCGAAGCAACGCATAACAGCATGCTGGTGGAGCTGTTCCGTCAGTCCTGGCAGTGGC‑GGGAAAACAATCCAATGT‑GGATCCAGTTGCACAGCCATCTGGATGACAGCCTGTATCGCAAAGAGTGGTTGGGCGATCACAAACAGATCCTCGCCGCGTTAATCAAAAAAGATGCCCGAGCGGCGAAGCTGGCAATGTGGCAGCATCTGGAAA  >  NC_000913/4554959‑4555251

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