Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,510,395	0	C	A	48.4%	54.5 / 41.9	62	S131R (AGC→AGA)	ydcR	putative DNA‑binding transcriptional regulator and putative aminotransferase
Reads supporting (aligned to +/- strand):  ref base C (23/9);  new base A (0/30);  total (23/39)
Fisher's exact test for biased strand distribution p-value = 5.04e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACTCGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTC‑GGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCCAAC‑AACTAAACCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCATT  >  NC_000913/1510246‑1510541                                                                                                                                                          |                                                                                                                                                    aCTCGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTC‑GGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCCAAC‑AACTAAACCGCTCGCTGGCGCAGGT‑AAGc                                                                                                                                                     >  1:601725/1‑150 (MQ=255)     gagaTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTC‑GGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCCAAC‑AACTAAACCGCTCGCTGGCGCAGGT‑AAGCaaaa                                                                                                                                                 <  2:173198/150‑1 (MQ=255)     gagaTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTC‑GGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCCAAC‑AACTAAACCGCTCGCTGGCGCAGGT‑AAGCAAAA‑c                                                                                                                                               <  1:162094/151‑1 (MQ=255)     gagaTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTC‑GGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCCAAC‑AACTAAACCGCTCGCTGGCGCAGGT‑AAGCAAAA‑c                                                                                                                                               >  1:587853/1‑151 (MQ=255)     gagaTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTC‑GGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCCAAC‑AACTAAACCGCTCGCTGGCGCAGGT‑AAGCAAAA‑c                                                                                                                                               >  1:490763/1‑151 (MQ=255)      agaTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTC‑GGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCCAAC‑AACTAAACCGCTCGCTGGCGCAGGT‑AAGCAAAA‑cc                                                                                                                                              >  2:826694/1‑151 (MQ=255)       gaTGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTC‑GGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCCAAC‑AACTAAACCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCg                                                                                                                                             <  1:828372/151‑1 (MQ=255)           aaGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTC‑GGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCCAAC‑AACTAAACCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCca                                                                                                                                          <  2:577413/150‑1 (MQ=255)                  cGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTC‑GGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCCAAC‑AACTAAACCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATg                                                                                                                                  >  1:562646/1‑151 (MQ=255)                      atCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTC‑GGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCCAAC‑AACTAAACCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCg                                                                                                                              >  1:552369/1‑151 (MQ=255)                       tCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTC‑GGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCCAAC‑AACTAAACCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGt                                                                                                                             >  2:924451/1‑151 (MQ=255)                               tataTTTTTGATTTGTTGCAGCCCAGTCGCGATCCGTCGGTCGTCCCTTTTGCCTC‑GGCCTTC‑CCCGACCCGCGATTTTTCCCCCCCCAAC‑AACTAACCCGCTCGCTGGCGCAGGT‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGaa                                                                                                                      <  1:919442/150‑1 (MQ=255)                                     ttttattttttGCTGTCCAGTCGTAACCCGTCTTTCGTTCCTTTTGCCTG‑GTCTTTT‑CCCGCCCCGCTTTTTTCCCCCCCCCAAC‑ACCAATCCCGCTCGCGGCCGCATGA‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTa                                                                                                                <  2:910842/141‑1 (MQ=255)                                            tGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTC‑GGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCCAAC‑AACTAAACCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAgg                                                                                                        >  1:855708/1‑151 (MQ=255)                                             gTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTC‑GGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCCAAC‑AACTAAACCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAgg                                                                                                        >  1:306378/1‑150 (MQ=255)                                               tGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTC‑GGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCCAAC‑AACTAAACCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGaa                                                                                                      >  1:205762/1‑150 (MQ=255)                                               tGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTC‑GGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCCAAC‑AACTAAACCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGaa                                                                                                      >  1:406996/1‑150 (MQ=255)                                                gCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTC‑GGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCCAAC‑AACTAAACCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAAc                                                                                                    >  2:18577/1‑151 (MQ=255)                                                 caccccAGCCGCGCTCCGTCGGTCGTCCCTTTTCCCCC‑GGCCTTC‑CCCCCCCCGCGATTTCCCCCCCCCCAAA‑AACAAACCCCCTCGCGGGCGCAGGA‑AACAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAgg                                                                                                        <  2:145031/142‑1 (MQ=255)                                                  aGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTC‑GGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCCAAC‑AACTAAACCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGc                                                                                                  >  1:336808/1‑151 (MQ=255)                                                   ggCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTC‑GGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCCAAC‑AACTAAACCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGCGCGTGATTGAAAACTTACCGCCAGGAAACGc                                                                                                  >  1:325760/1‑150 (MQ=255)                                                   ggCCAGCCGCGATCCGTCGGTCTTTCCGTTTGCCTC‑GCCTTTC‑CCCGCCCCGCGATTTTCCCCCCCCCAAC‑AACTAACCCGCTCGCGGGCGAAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGaaa                                                                                                     <  1:200481/147‑1 (MQ=255)                                                     ccatccgcgTTCCGCCGTTCTTCCCGTTT‑CCCG‑GGCCTTTCCCCCCCCCGGGATTTTTCCCCCCCAAAA‑AAAAAACCCGCTCGCGGGCGCAGAA‑AAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAAc                                                                                                    <  1:748620/141‑1 (MQ=255)                                                                     ggtcg‑TCCGTTTTCCCCGGGCCTTC‑CCCGCCCCGCG‑CTTTTCCCCCCCCCACAAACAAACCCGCTCGCGGGCG‑AGGTAAAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGt                                                                                        <  2:981435/142‑1 (MQ=255)                                                                        cgTTCCGTTTGCCTC‑GGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCCAAC‑AACTATACCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCt                                                                            >  1:1012427/1‑151 (MQ=255)                                                                         gTTCCGTTTGCCTC‑GGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCCAAC‑AACTAAACCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCtc                                                                           >  2:160913/1‑151 (MQ=255)                                                                          tccgttttGCCCC‑GCCTTTC‑CCCCCCCCGCGTTTTTCCCCCCCCCAAC‑AACAAACCCGCCCGCGGGCGGAGAA‑AAAAAAAC‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTAtt                                                                               <  2:656089/141‑1 (MQ=255)                                                                          tcccGTTTCCCGC‑GCCCTTC‑CCCGCCCGCCGACTTTTCCCCCCCCACC‑AACTAACCCCCCCGCGGGCGCAGTA‑AGGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTg                                                                              <  1:100461/145‑1 (MQ=255)                                                                            tcccttttccctggCCTTTCCCCGCCCCTCGTTTTTCCCCCCACCAAA‑AACAAACCCGCCGGCGGGCGCAGAA‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCt                                                                            <  2:586839/137‑1 (MQ=255)                                                                              ccttttcccgggCCTTTCCCCGCCCCGCGATTTTCCCCCCCCCAAC‑AACAAACCCGCTGGCGGGGGGGGAA‑AAAAAAAT‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGc                                                                             <  1:341830/136‑1 (MQ=255)                                                                                   ccTC‑GCCCTTC‑CCCCCCCCGCGATTTTCCCCCCCCCAAC‑AAAAAACCCGCTCGCTGGCGCAGTA‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTAt                                                                   <  1:795768/149‑1 (MQ=255)                                                                                    ccccgcccTTC‑CCCGCCCCGCGTCTTTCCCCCCCCCAAC‑AACAAACCCGCTCGCGGGCGC‑GGT‑AAAAAAAACCCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCtcgtcg                                                                       <  1:778704/139‑1 (MQ=255)                                                                                          ccTTT‑CCCGACCCGCGACTTTTCCCCCTCCAAC‑AACTAAACCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTAc                                                            >  2:218265/1‑150 (MQ=255)                                                                                              t‑CCCGACCCGCGACTTTTCCCCCTCCAAC‑AACTAAACCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGc                                                       >  2:896518/1‑151 (MQ=255)                                                                                                   gACCCGCGACTTTTCCCCCTCCAAC‑AACTAAACCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATc                                                    >  1:715616/1‑150 (MQ=255)                                                                                                     cccGCGACTTTCCCCCCCCAAAC‑AACAAACCCGCTCGCTGGCGAAGGA‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCAc                                                  <  1:127226/150‑1 (MQ=255)                                                                                                      ccgcgttattcccccccccAAC‑AACAAAGCTGCCCGCGGGCGGAGAA‑AATCAAAC‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCa                                                <  2:663291/140‑1 (MQ=255)                                                                                                        gcgttttccccccccccAAA‑AACAAACCCGCCCGCGGGCGCAGGA‑AACAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCAtt                                              <  2:322585/142‑1 (MQ=255)                                                                                                            cTTTTCCCCCTCCAAC‑AACTAAACCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTtctc                                          >  1:40762/1‑151 (MQ=255)                                                                                                             tttcCCCCCCCCAAC‑AACAAACCCCCTCGCGGGCGCGGGA‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTtct                                           <  2:126919/149‑1 (MQ=255)                                                                                                                 ccccccccAAC‑AACTAACCCGCTCGCGGCCGAAGTT‑AAAAAAAC‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCt                                        <  2:637430/148‑1 (MQ=255)                                                                                                                 ccccccccAAC‑AAAAAACCCGCTGGCGGGGGCGGTT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCtgat                                     <  2:562597/151‑1 (MQ=255)                                                                                                                 ccccccccAAA‑AACAAACCCGCTCGCGGGCGCAGGA‑AAACAAAC‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCtg                                       <  2:546909/149‑1 (MQ=255)                                                                                                                  ccccccCAAC‑AACAAACCCGCCCGCGGGCGCAGGA‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCtgatg                                    <  2:379918/151‑1 (MQ=255)                                                                                                                  ccccccCAAA‑AACAAACCCGCCCGCTGGCGCAGGA‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCtgatg                                    <  1:774941/151‑1 (MQ=255)                                                                                                                   ccccccAAC‑AAAAAACCCGCCCGCGGGCCCAGGA‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCtgatga                                   <  1:74825/151‑1 (MQ=255)                                                                                                                    cccccAAC‑AACAAACCCGCTCGCGGGCGCAGGT‑AAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGaa                                  <  1:470514/148‑1 (MQ=255)                                                                                                                       ccAAC‑AACAAACCCCCGCGGTGGCGCAGGA‑AAGAAAAC‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAAtt                               <  2:38257/151‑1 (MQ=255)                                                                                                                             aacAAACCCGCTCGGGGGCGCAGGA‑AAAAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTc                            <  1:849307/149‑1 (MQ=255)                                                                                                                             aacAAACCCGCCCCCGGGCACAGAA‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGt                             <  2:978312/148‑1 (MQ=255)                                                                                                                                tAAACCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTAc                       >  2:685776/1‑151 (MQ=255)                                                                                                                                  aaCCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACt                      >  2:600962/1‑150 (MQ=255)                                                                                                                                   cccGCTCGCTGGCGAAGAA‑AAGAAAAC‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTg                     <  2:747713/149‑1 (MQ=255)                                                                                                                                    ccgcccgcTGGCGAGGGA‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGcc                   <  1:784998/151‑1 (MQ=255)                                                                                                                                    cccctcgcgGGCGCAGGA‑AAGAAAAC‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGc                    <  2:685604/147‑1 (MQ=255)                                                                                                                                        tcgctGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCgggg               >  2:184846/1‑151 (MQ=255)                                                                                                                                             ggCGCAGGT‑AAGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGc              <  1:378077/147‑1 (MQ=255)                                                                                                                                             ggCGCAGGT‑AACAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGtt           <  1:894398/150‑1 (MQ=255)                                                                                                                                             ggCGCAGGA‑AAAAAAAC‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGtt           <  1:637503/150‑1 (MQ=255)                                                                                                                                              ggcgcggaaaaGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTa          <  2:336782/142‑1 (MQ=255)                                                                                                                                                    gaaaaGAAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAgg      <  2:896764/146‑1 (MQ=255)                                                                                                                                                     t‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCAtt  >  1:465933/1‑151 (MQ=255)                                                                                                                                                      aaaGAAAAC‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAgg      <  2:537439/146‑1 (MQ=255)                                                                                                                                                          |                                                                                                                                                    ACTCGAGATGAAGCAGTCGATATCAACACTTATATTTTTGATATGTTGCAGGCCAGTCGCGATCCGTCGGTCGTTCCGTTTGCCTC‑GGCCTTT‑CCCGACCCGCGACTTTTCCCCCTCCAAC‑AACTAAACCGCTCGCTGGCGCAGGT‑AAGCAAAA‑CCGCCACAGCGATGAGCGTGATTGAAAACTTACCGCCAGGAAACGCAGAACTGCGTCAGGCTATTGCTCGTCGCTATGCCTTACAGGGCATCACCATTTCTCCTGATGAAATTGTCATTACTGCCGGGGCGTTAGAGGCATT  >  NC_000913/1510246‑1510541

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