Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	517,031	A→C	6.6%	intergenic (‑40/+28)	PP_0429 ← / ← PP_0430	short‑chain dehydrogenase/reductase family oxidoreductase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	517,031	0	A	C	6.6%	176.9 / 2.1	57	intergenic (‑40/+28)	PP_0429/PP_0430	short‑chain dehydrogenase/reductase family oxidoreductase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (28/25);  new base C (2/2);  total (30/27)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGCAGGCGCTCGATCAATGCGCTGGCGCTGGGCAGGTCATGGATGGCACCGACGCCCTGGCCC  >  NC_002947/516886‑517179                                                                                                                                                  |                                                                                                                                                      gTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACtt                                                                                                                                                    >  2:86462/1‑149 (MQ=255) gTTGGCGTGCCACCAGGATAAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACtt                                                                                                                                                    >  1:224194/1‑149 (MQ=255)    ggCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACACCTTGca                                                                                                                                                 >  3:98786/1‑148 (MQ=255)           caccaGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGc                                                                                                                                            >  6:129571/1‑147 (MQ=255)           caccaGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCaa                                                                                                                                          <  3:109053/149‑1 (MQ=255)           caccaGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCaa                                                                                                                                          <  2:83781/149‑1 (MQ=255)            accaGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAg                                                                                                                                         >  5:101624/1‑149 (MQ=255)                 gATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACACCTTGCACGGGCAAGttgta                                                                                                                                    <  2:75772/149‑2 (MQ=255)                  aTCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTg                                                                                                                                   <  8:191488/149‑1 (MQ=255)                       aTTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCttt                                                                                                                              <  5:227696/149‑1 (MQ=255)                        ttGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCttt‑t                                                                                                                            <  4:180543/149‑1 (MQ=255)                           cgcCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGc                                                                                                                          <  5:196040/148‑1 (MQ=255)                            gcCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAg                                                                                                                                         <  4:132375/133‑1 (MQ=255)                            gcCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAg                                                                                                                                         >  3:132375/1‑133 (MQ=255)                            gcCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTa                                                                                                                        <  4:273392/149‑1 (MQ=255)                                gcgcCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAg                                                                                                                    >  5:112323/1‑149 (MQ=255)                                        ccAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAggcgcgggc                                                                                                            >  7:179089/1‑149 (MQ=255)                                        ccAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAggcgcgggc                                                                                                            >  7:24486/1‑149 (MQ=255)                                                   gCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCt                                                                                                 >  2:256825/1‑149 (MQ=255)                                                    ccAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTg                                                                                                >  1:53517/1‑149 (MQ=255)                                                       aggccaggccCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCa                                                                                             <  6:101624/149‑1 (MQ=255)                                                       aggccaggccCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCa                                                                                             <  2:239230/149‑1 (MQ=255)                                                        ggccaggccCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAg                                                                                            >  5:216682/1‑149 (MQ=255)                                                           caggcccaggccGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGc                                                                                         <  1:86462/149‑1 (MQ=255)                                                                ccaggccGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAg                                                                                    >  3:231263/1‑149 (MQ=255)                                                                 caggccGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAgg                                                                                   <  6:112323/149‑1 (MQ=255)                                                                    gccGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCt                                                                                >  7:263433/1‑149 (MQ=255)                                                                      cGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTg                                                                              >  1:284038/1‑149 (MQ=255)                                                                       gCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACACCTTGCACGGGCAAg                                                                                                                                         >  2:213390/1‑90 (MQ=255)                                                                       gCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGc                                                                             >  3:195617/1‑149 (MQ=255)                                                                       gCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGc                                                                             >  1:303212/1‑149 (MQ=255)                                                                       gCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGCGCGCCTAGTGTACACCTTGCACGGGCAAg                                                                                                                                         <  1:213390/90‑1 (MQ=255)                                                                            aaGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTa                                                                        <  5:129571/149‑1 (MQ=255)                                                                                 ccAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGc                                                                   <  2:284038/149‑1 (MQ=255)                                                                                     tGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATg                                                               >  3:55515/1‑149 (MQ=255)                                                                                       aTCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGc                                                              <  8:179089/148‑1 (MQ=255)                                                                                       aTCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGc                                                              <  8:24486/148‑1 (MQ=255)                                                                                        tCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGCAg                                                            >  6:115409/1‑149 (MQ=255)                                                                                         cATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGCAgg                                                           >  3:97389/1‑149 (MQ=255)                                                                                               gTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGCAGGCGCTc                                                      >  5:33932/1‑148 (MQ=255)                                                                                                tAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGCAGGCGCTCGa                                                    <  2:303212/149‑1 (MQ=255)                                                                                                          tGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGCAGGCGCTCGATCAATgcgct                                          <  2:53517/149‑1 (MQ=255)                                                                                                                    cTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGCAGGCGCTCGATCAATGcgctggcgctgg                                  <  4:231263/147‑1 (MQ=255)                                                                                                                    cTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGCAGGCGCTCGATCAATGCGCTGGCGCTtgg                                 <  4:97389/148‑1 (MQ=255)                                                                                                                     tGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGCAGGCGCTCGATCAATGCGCTGGCGCTggg                                 >  8:37360/1‑147 (MQ=255)                                                                                                                      gggggCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGCAGGCGCTCGATCAATGCGCTGGCGCTGGGCAg                              >  5:208667/1‑149 (MQ=255)                                                                                                                       ggggCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGCAGGCGCTCGATCAATGCGCTGGCGCTGGGCa                               <  4:195617/147‑1 (MQ=255)                                                                                                                                   gggggcgaGTGTACAACTTTGCCGTTCCATGTGGTGCCTTTCCAGTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGCAGGCGCTCGATCAATGCGCTGGCGCTGGGCAGGTCatggatgg                   <  1:256825/141‑1 (MQ=255)                                                                                                                                         tAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGCAGGCGCTCGATCAATGCGCTGGCGCTGGGCAGGTCATGGATGGCACCGACg           >  7:44895/1‑149 (MQ=255)                                                                                                                                          aGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGCAGGCGCTCGATCAATGCGCTGGCGCTGGGCAGGTCATGGATGGCACCGACGc          >  4:215260/1‑149 (MQ=255)                                                                                                                                              tACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGCAGGCGCTCGATCAATGCGCTGGCGCTGGGCAGGTCATGGATGGCACCGACGCCCTg      >  6:123044/1‑149 (MQ=255)                                                                                                                                              tACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGCAGGCGCTCGATCAATGCGCTGGCGCTGGGCAGGTCATGGATGGCACCGACGCCCTg      <  6:33932/149‑1 (MQ=255)                                                                                                                                              tACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGCAGGCGCTCGATCAATGCGCTGGCGCTGGGCAGGTCATGGATGGCACCGACGCCCTg      >  8:132224/1‑149 (MQ=255)                                                                                                                                              tACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGCAGGCGCTCGATCAATGCGCTGGCGCTGGGCAGGTCATGGATGGCACCGACGCCCTg      >  3:108379/1‑149 (MQ=255)                                                                                                                                                cAACTTGCCCGTGCAAGGTGTTTCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGCAGGCGCTCGATCAATGCGCTGGCGCTGGGCAGGTCATGGATGGCACCGACTCCCTgg     <  6:208667/148‑1 (MQ=255)                                                                                                                                                cAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGCAGGCGCTCGATCAATGCGCTGGCGCTGGGCAGGTCATGGATGGCACCGACGCCCTGGc    >  5:295213/1‑149 (MQ=255)                                                                                                                                                cAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGCAGGCGCTCGATCAATGCGCTGGCGCTGGGCAGGTCATGGATGGCACCGACGCCCTGGc    <  5:123044/149‑1 (MQ=255)                                                                                                                                                 aaCTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGCAGGCGCTCGATCAATGCGCTGGCGCTGGGCAGGTCATGGATGGCACCGACGCCCTGGcc   <  8:263433/149‑1 (MQ=255)                                                                                                                                                  aCTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGCAGGCGCTCGATCAATGCGCTGGCGCTGGGCAGGTCATGGATGGCACCGACGCCCTGGccc  >  8:214047/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                      GTTGGCGTGCCACCAGGATCAGATTGCGCCCGCGCCGTGCCAGCGCTTCCGCCAAGGCCAGGCCCAGGCCGCTGGAAGCACCAGTGATCATGGCGTAACGGGTCATGCAAGGCTCCTGGGGGCGGAAGAGGGCGCCTAGTGTACAACTTGCCCGTGCAAGGTGTTGCCTTT‑TGCTACAAGGCGCGGGCGCGCAGAACCTGGCAGTGCTCCAGGGCTTGCCGGTACTCGCCATGCAGGCGCTCGATCAATGCGCTGGCGCTGGGCAGGTCATGGATGGCACCGACGCCCTGGCCC  >  NC_002947/516886‑517179

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