Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,631,093	T→C	6.5%	intergenic (+27/‑154)	PP_1430 → / → lepA	DegP‑like serine endoprotease/elongation factor 4

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,631,093	0	T	C	6.5%	173.1 / 2.1	61	intergenic (+27/‑154)	PP_1430/lepA	DegP‑like serine endoprotease/elongation factor 4
Reads supporting (aligned to +/- strand):  ref base T (26/31);  new base C (2/2);  total (28/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

ATCGGTTCGGCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCTTTTCGAGTGT  >  NC_002947/1630947‑1631240                                                                                                                                                    |                                                                                                                                                    aTCGGTTCGGCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTc                                                                                                                                                   >  3:167114/1‑149 (MQ=255)    ggTTCGGCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCgg                                                                                                                                                 <  8:262032/148‑1 (MQ=255)         ggCCAAGGCATTCACCGAAATCGCCACCGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGccc                                                                                                                                           <  8:331740/149‑1 (MQ=255)          gCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCCGAAGCTGCCCt                                                                                                                                          >  4:220803/1‑149 (MQ=255)            cAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCttt                                                                                                                                        <  5:22025/149‑1 (MQ=255)             aaGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCCGAAGCTGCCCtttc                                                                                                                                       <  3:220803/149‑2 (MQ=255)               ggAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTc                                                                                                                                     <  7:248202/149‑1 (MQ=255)                gAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCCGAAGCTGCCCtttcgcc                                                                                                                                    <  2:38571/149‑5 (MQ=255)                 aaTTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCAt                                                                                                                                   <  4:21730/149‑1 (MQ=255)                  aTTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGcc                                                                                                                                            >  6:106264/1‑139 (MQ=255)                  aTTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGcc                                                                                                                                            <  5:106264/139‑1 (MQ=255)                        cGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGccc                                                                                                                                           >  8:259330/1‑134 (MQ=255)                        cGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGccc                                                                                                                                           <  7:259330/134‑1 (MQ=255)                        cGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCtt                                                                                                                            >  4:83531/1‑149 (MQ=255)                             tCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACaa                                                                                                                       >  3:31169/1‑149 (MQ=255)                                ccAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTc                                                                                                                           >  4:60284/1‑142 (MQ=255)                                ccAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTc                                                                                                                           <  3:60284/142‑1 (MQ=255)                                 cAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGc                                                                                                                   >  1:238327/1‑149 (MQ=255)                                  aaGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGcc                                                                                                                  <  4:31169/149‑1 (MQ=255)                                                gaaTCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCtgt                                                                                                    <  6:42997/149‑1 (MQ=255)                                                 actCGTTCCGTGGCGATGCGCGGACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TATGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCt                                                                                                      <  2:130011/144‑1 (MQ=255)                                                 aaTCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TATGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATCAATCTTCACAATTGCCGGAATGATGGcg                                                                                                      >  1:130011/1‑145 (MQ=255)                                                   tCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTa                                                                                                 <  3:83531/149‑1 (MQ=255)                                                      ttCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAAGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAg                                                                                              >  4:248157/1‑149 (MQ=255)                                                        ccGTGTCGATGCGCGTACTGCGTCAGGGCCGGTCCAGCTGCATCACATTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGt                                                                                             <  6:342248/148‑1 (MQ=255)                                                             tCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGAc                                                                                        <  7:168302/148‑1 (MQ=255)                                                                aTGCGCGTACTGCGTCAGGGCCGTGCCAGCTTAATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAATGTGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACagag                                                                                    <  3:21632/149‑1 (MQ=255)                                                                   cgcgTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGa                                                                                  <  6:121396/148‑1 (MQ=255)                                                                       tACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAAtctc                                                                             <  4:167114/149‑1 (MQ=255)                                                                        aCTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAAGTAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCtgt                                                                                                    >  1:260423/1‑126 (MQ=255)                                                                        aCTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATGTTCACATTTGCCGGAATGTTTGCGGGGTACAGTCGGAGAGAGGaagatcc                                                                            >  4:208273/1‑144 (MQ=255)                                                                        aCTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCtgt                                                                                                    <  2:260423/125‑1 (MQ=255)                                                                                    ccGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGTAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCcagca                                                                <  2:238327/149‑1 (MQ=255)                                                                                             cTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCt                                                                                                      >  6:41119/1‑102 (MQ=255)                                                                                             cTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCt                                                                                                      <  5:41119/102‑1 (MQ=255)                                                                                                  tcaCTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAAt                                                  >  4:7006/1‑149 (MQ=255)                                                                                                  tcaCTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAAt                                                  >  2:223046/1‑149 (MQ=255)                                                                                                      tttCTAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCcagca                                                                >  4:47754/1‑131 (MQ=255)                                                                                                      tttCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCcagca                                                                <  3:47754/131‑1 (MQ=255)                                                                                                      tttCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTccc                                              >  1:272522/1‑149 (MQ=255)                                                                                                         cAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCgg                                            <  6:131295/148‑1 (MQ=255)                                                                                                                 cTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCcagca                                                                <  4:99605/120‑1 (MQ=255)                                                                                                                 cTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCcagca                                                                >  3:99605/1‑120 (MQ=255)                                                                                                                  tGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTcg                                  >  7:154668/1‑149 (MQ=255)                                                                                                                     a‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTcggcg                               >  7:132627/1‑149 (MQ=255)                                                                                                                        aaGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCggcggg                             >  6:192294/1‑149 (MQ=255)                                                                                                                           cAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTa                                                      >  4:118593/1‑121 (MQ=255)                                                                                                                           cAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTa                                                      <  3:118593/121‑1 (MQ=255)                                                                                                                              tttCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTccc                                                                           >  7:306442/1‑97 (MQ=255)                                                                                                                              tttCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTccc                                                                           <  8:306442/97‑1 (MQ=255)                                                                                                                                  gCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTc                   >  1:286386/1‑149 (MQ=255)                                                                                                                                   cAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGTCGTCCGGCTCg                  >  8:280122/1‑149 (MQ=255)                                                                                                                                   cAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCg                  >  4:153094/1‑149 (MQ=255)                                                                                                                                       gCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGc              <  8:132627/149‑1 (MQ=255)                                                                                                                                         gAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCt            >  4:104231/1‑149 (MQ=255)                                                                                                                                         aaaaGGGCAGCCGAAGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCt            >  3:135883/2‑149 (MQ=255)                                                                                                                                            aGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCtttt         <  2:272522/149‑1 (MQ=255)                                                                                                                                               gCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCTTTTCGa      <  5:192294/149‑1 (MQ=255)                                                                                                                                               gCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCTTTTCGa      <  3:7006/149‑1 (MQ=255)                                                                                                                                               gCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCTTTTCGa      <  8:154668/149‑1 (MQ=255)                                                                                                                                                   cTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTACGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCTTTTCGAgtgt  <  2:286386/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                    ATCGGTTCGGCCAAGGAATTCACCGAAATCGCCAAGGAACTGCCGAAGAATCGTTCCGTGTCGATGCGCGTACTGCGTCAGGGCCGTGCCAGCTTCATCACTTTCAAGCTTGCTGAA‑TAAGCAGTTTCGCAAGGCGAAAGGGCAGCTTCGGCTGCCCTTTTTCATGAATCTTCACAATTGCCGGAATGATTGCTGTGTACAGTCGGACAGAGGAATCTCCCATAACCCAGCAGCTTAAGGTACAATTCCCGGCTATTTTTCGGCGGGCGTCCGGCTCGCAGCCTTTTCGAGTGT  >  NC_002947/1630947‑1631240

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