Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,893,846	Δ1 bp	6.2%	coding (1213/1572 nt)	yegH ←	putative transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,893,846	0	C	.	6.2%	258.1 / 10.5	80	coding (1213/1572 nt)	yegH	putative transporter
Reads supporting (aligned to +/- strand):  ref base C (41/34);  new base . (2/3);  total (43/37)
Fisher's exact test for biased strand distribution p-value = 6.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.83e-01

CTGAGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAG  >  NC_002947/3893701‑3893987                                                                                                                                                   |                                                                                                                                              cTGAGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGtt                                                                                                                                            <  1:6487/149‑1 (MQ=255)  tGAGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATACGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGACATGGTAAAC‑AGCCCGCTCAAGTCACCGAAATCGtt                                                                                                                                            >  7:253618/1‑148 (MQ=255)    aGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGtt                                                                                                                                            <  2:228948/146‑1 (MQ=255)    aGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGtt                                                                                                                                            >  1:228948/1‑146 (MQ=255)      gTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCa                                                                                                                                       >  2:35555/1‑149 (MQ=255)       ttCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACGAGTTCGGTGACTTCACCGGGCT‑GTTGACCat                                                                                                                                      <  3:212457/149‑2 (MQ=37)        tCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAT‑AGCCCGGTGAAGTCACCGAACTCGTTGACCAc                                                                                                                                      >  1:346940/1‑148 (MQ=255)        tCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCAc                                                                                                                                      >  1:264688/1‑148 (MQ=255)          aGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACGAGTTCGGTGACTTCACCGGGCT‑GTTGACCATGa                                                                                                                                    <  2:363340/148‑1 (MQ=37)          aGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGaa                                                                                                                                   <  2:99897/149‑1 (MQ=255)             gcACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGc                                                                                                                                <  2:240139/149‑1 (MQ=255)              cACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACGAGTTCGGTGACTTCACCGGGCT‑GTTGACCat                                                                                                                                      <  3:100155/142‑2 (MQ=37)              cACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAACGAGTTCGGTGACTTCACCGGGCT‑GTTGACCat                                                                                                                                      >  4:100155/1‑141 (MQ=37)                ccGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAat                                                                                                                             <  4:240520/149‑1 (MQ=255)                  gCTGACCACAAAGCCTTCCGCCTCCTGAACGATGCCCTGGGCTTCAACCTCACTGGCATGAGGCAGCTCACCGGGAATGGACTCAAGGATGTCGGTCATGGTGTTC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGTAGGCAatat                                                                                                                           <  8:39677/149‑1 (MQ=255)                  gCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAT‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAatat                                                                                                                           <  2:346940/149‑1 (MQ=255)                           aaaGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATTGGTCGACt                                                                                                                  >  3:183640/1‑149 (MQ=255)                           aaaGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGGCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACt                                                                                                                  >  3:267471/1‑149 (MQ=255)                                 ttCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACcgc                                                                                                            >  7:286899/1‑149 (MQ=255)                                     cccTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCAt                                                                                                        >  4:361595/1‑149 (MQ=255)                                      cctcctGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCAtt                                                                                                       >  7:240782/1‑149 (MQ=255)                                      cctcctGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCAtt                                                                                                       <  2:264688/149‑1 (MQ=255)                                      cctcctGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCAtt                                                                                                       <  8:62886/149‑1 (MQ=255)                                       ctcctGAACGATGCCAGGGCCTTCAACATCAGTGGCATCCGGCAGCGCACCGGCAATCGAGTCAAGCATGTCGGTAATGGTCAAC‑AGCCCGGTGAAGACACCGAACTCGTTGACCACGCAGGCAAGATGGGTCGCCTGACCGCGCAttt                                                                                                      >  1:279723/1‑149 (MQ=255)                                          ctGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCt                                                                                                   >  2:344419/1‑149 (MQ=255)                                          ctGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGAATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCt                                                                                                   >  7:235419/1‑149 (MQ=255)                                            gAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTc                                                                                                  >  2:73142/1‑148 (MQ=255)                                                gATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAggg                                                                                             >  4:231071/1‑149 (MQ=255)                                                    cccGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCAtt                                                                                         <  1:226834/149‑1 (MQ=255)                                                        ggCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAgg                                                                                     <  8:240782/149‑1 (MQ=255)                                                           cTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCAtt                                                                                         >  6:156485/1‑142 (MQ=255)                                                           cTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCAtt                                                                                         <  5:156485/142‑1 (MQ=255)                                                                   tCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGATGATGCTGAAGct                                                                          >  3:250107/1‑149 (MQ=255)                                                                      cTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCtctc                                                                       <  8:235419/149‑1 (MQ=255)                                                                      cTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCtctc                                                                       <  1:344419/149‑1 (MQ=255)                                                                       tGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTc                                                                                                  <  4:49309/121‑1 (MQ=255)                                                                       tGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTc                                                                                                  >  3:49309/1‑121 (MQ=255)                                                                       tGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCtctc                                                                       >  6:109809/1‑148 (MQ=255)                                                                         gCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTg                                                                               <  5:215686/138‑1 (MQ=255)                                                                         gCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTg                                                                               >  6:215686/1‑138 (MQ=255)                                                                           aTCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATTGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCcagca                                                                  <  4:183640/149‑1 (MQ=255)                                                                            tCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGa                                                                                                                    >  4:244546/1‑98 (MQ=255)                                                                            tCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGa                                                                                                                    <  3:244546/98‑1 (MQ=255)                                                                               ggCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGACGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGCAGCTCACCAGCAGGtt                                                              >  2:315850/1‑149 (MQ=255)                                                                                       aCCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGgcg                                                      >  6:269094/1‑149 (MQ=255)                                                                                         cGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGcc                                                    >  8:245320/1‑149 (MQ=255)                                                                                             aaTCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCggg                                                >  6:37961/1‑149 (MQ=255)                                                                                                 gACTCAAGGATGTCGGTCATGGTTAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCAt                                            >  3:122675/1‑149 (MQ=255)                                                                                                      aaGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACTACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCttt                                       >  1:176576/1‑149 (MQ=255)                                                                                                      aaGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCttt                                       >  1:305179/1‑149 (MQ=255)                                                                                                      aaGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCttt                                       >  6:150813/1‑149 (MQ=255)                                                                                                        ggATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTc                                      >  4:356649/1‑148 (MQ=255)                                                                                                            gTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGt                                 >  8:175336/1‑149 (MQ=255)                                                                                                                ttcgtGGTCAACGAGTTCGGTGACTTCACCGGGCT‑GTTGACCATGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCGAGCAGGTTCAACGGCGCCCGGGCCATGCTTGCCAGGTCgcg                             >  4:264668/5‑147 (MQ=25)                                                                                                                gTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCggg                             <  1:35555/149‑1 (MQ=255)                                                                                                                 tCATCGTCAAC‑AGCCCCGTGAAGTCACCGACCTCGTTGACCACGAAGGCAATCTGGGTCGACTGACCGCGCATTTGCGGCTGTGCATTTAGCCTCCTGCAGCTCTCCAGCAGGCTCAACGTCGCCCGGGCCATGCTTTCCAGGTCGggc                            <  2:279723/149‑1 (MQ=255)                                                                                                                    tGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGggctgg                         <  6:236880/149‑1 (MQ=255)                                                                                                                      gTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGggctggct                       >  2:162017/1‑149 (MQ=255)                                                                                                                         aaC‑AGCCCGGTGAAGTCAACGAACTCGTTGACAACGAAGGCAATATTTGTCGACTGACCGCGCATTTGCTCAAGGGCATTGAGGATGCTGATGCTATCCAGCAGGTTAAAAGGAGCCCTGGCCATGCTTTCCAGGTCGGGCTTGCTGcc                    >  2:233898/1‑149 (MQ=255)                                                                                                                          aC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGccc                   <  7:245320/149‑1 (MQ=255)                                                                                                                          aC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGccc                   <  8:286899/149‑1 (MQ=255)                                                                                                                           c‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGccc                   >  5:359673/1‑148 (MQ=255)                                                                                                                             aGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCga                 <  1:162017/149‑1 (MQ=255)                                                                                                                               cccGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTCCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCAGCCAGCCCgacc               <  3:231071/149‑3 (MQ=255)                                                                                                                                ccgtttaaGTCACCGAACTCGTTGACAACGAAGGCAATATGGGGGGACTGACCGCGCATTTGCTCCAGGGCATTGTGGATGCTGAAGCTCTCCTGCAGGTTCAACTGCGACCGGGCCATTCTTTCCAGGTCGGGCTGTCTGCCCGagag              <  8:253618/143‑1 (MQ=255)                                                                                                                                ccGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGagag              <  6:359673/149‑1 (MQ=255)                                                                                                                                  ggTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagca            <  7:175336/149‑1 (MQ=255)                                                                                                                                  ggTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagca            <  3:361595/149‑1 (MQ=255)                                                                                                                                       aGTCACCGAACTCGTTGACTACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttcc       <  2:176576/149‑1 (MQ=255)                                                                                                                                        gTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttcct      >  4:296480/1‑149 (MQ=255)                                                                                                                                         tCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttcctt     <  5:269094/149‑1 (MQ=255)                                                                                                                                         tCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttcctt     <  3:356649/149‑1 (MQ=255)                                                                                                                                           aCCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCggg                             >  2:339510/1‑123 (MQ=255)                                                                                                                                           aCCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCggg                             <  1:339510/123‑1 (MQ=255)                                                                                                                                           aCCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttcctt     >  1:90083/1‑147 (MQ=255)                                                                                                                                           aCCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttcctt     >  6:64489/1‑147 (MQ=255)                                                                                                                                             cGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCGTTCCAGGTCGGGCTGGCTGCACGAGagccgttgcttcag  >  8:343228/1‑148 (MQ=255)                                                                                                                                              gAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGagcagttccttcag  >  5:279622/1‑147 (MQ=255)                                                                                                                                                  tCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTg                                                                               <  7:297115/66‑1 (MQ=255)                                                                                                                                                  tCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTg                                                                               >  8:297115/1‑66 (MQ=255)                                                                                                                                                   |                                                                                                                                              CTGAGGTTCAGCGCACCGCTGACCACAAAGCCTTCCCCCTCCTGAACGATGCCCGGGCCTTCAACCTCACTGGCATCCGGCAGCTCACCGGCAATCGACTCAAGGATGTCGGTCATGGTCAAC‑AGCCCGGTGAAGTCACCGAACTCGTTGACCACGAAGGCAATATGGGTCGACTGACCGCGCATTTGCTCCAGGGCATTGAGGATGCTGAAGCTCTCCAGCAGGTTCAACGGCGCCCGGGCCATGCTTTCCAGGTCGGGCTGGCTGCCCGAGAGCAGTTCCTTCAG  >  NC_002947/3893701‑3893987

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: ‑