Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,258,516	Δ1 bp	6.3%	coding (180/1005 nt)	asd →	aspartate‑semialdehyde dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,258,516	0	C	.	6.3%	206.2 / 6.9	63	coding (180/1005 nt)	asd	aspartate‑semialdehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (30/29);  new base . (2/2);  total (32/31)
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.00e-01

GCCACCGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAA  >  NC_002947/2258370‑2258663                                                                                                                                                   |                                                                                                                                                     gcCACCGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAg                                                                                                                                                    >  7:117242/1‑149 (MQ=255)            cGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCg                                                                                                                                        <  7:266552/149‑1 (MQ=255)             gTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGc                                                                                                                                       >  5:155358/1‑149 (MQ=255)              tCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGcc                                                                                                                                      <  6:155358/149‑1 (MQ=255)              tCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGcc                                                                                                                                      <  4:90920/149‑1 (MQ=255)                tgtgtAGCGTTGGTGCAGATACTCGAAGAAGTGGCATTCCCGGTTGCCAGTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCTA‑AGCTGTCGACTTCGCGCa                                                                                                                                    <  2:93172/148‑1 (MQ=255)                  tGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAgg                                                                                                                                  >  6:32732/1‑149 (MQ=255)                  tGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAgg                                                                                                                                  >  4:236383/1‑149 (MQ=255)                   gAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGTCAAGAAACTGAAGTTGCGCGAAGTCGA‑AGCTGTCGACTTCgagatag                                                                                                                                  >  5:187861/1‑142 (MQ=255)                   gAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGTTGCGCGAAGTCGA‑AGCTGTCGACTTc                                                                                                                                         <  6:187861/141‑1 (MQ=255)                            gTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAg                                                                                                                             >  4:291592/1‑144 (MQ=255)                            gTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAg                                                                                                                             <  3:291592/144‑1 (MQ=255)                             tGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGcc                                                                                                                       <  3:165673/149‑1 (MQ=255)                              gCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCct                                                                                                                      >  4:64904/1‑149 (MQ=255)                                    ccTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCg                                                                                                                >  1:118925/1‑149 (MQ=255)                                      tCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGcc                                                                                                              >  3:167546/1‑149 (MQ=255)                                      tCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGc                                                                                                               >  7:301439/1‑148 (MQ=255)                                                  cATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCgccgcc                                                                                                  <  8:117242/149‑1 (MQ=255)                                                  cATTCCCGGTTGCCACTCTGCACCTACTGTCCAGCATGGAATCGGCGGGAAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCTCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCgccgcc                                                                                                  <  7:257676/149‑1 (MQ=255)                                                    ttCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTAGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGt                                                                                                >  3:326181/1‑149 (MQ=255)                                                     tCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGt                                                                                                >  4:128258/1‑148 (MQ=255)                                                       ccGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAg                                                                                             >  2:307219/1‑149 (MQ=255)                                                         ggTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGcc                                                                                           <  5:4770/149‑1 (MQ=255)                                                              ccACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGc                                                                                      >  6:4118/1‑149 (MQ=255)                                                               cACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCt                                                                                     >  2:27615/1‑149 (MQ=255)                                                                 ctctGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCttt                                                                                   <  5:256371/149‑1 (MQ=255)                                                                  tctGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTg                                                                                  >  5:141016/1‑149 (MQ=255)                                                                  tctGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTg                                                                                  >  3:211744/1‑149 (MQ=255)                                                                  tctGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTg                                                                                  >  2:34374/1‑149 (MQ=255)                                                                    tGCACCTGCTGGCCGGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGc                                                                                                                            >  6:61074/1‑105 (MQ=255)                                                                    tGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGc                                                                                                                            <  5:61074/105‑1 (MQ=255)                                                                     gCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCg                                                                               >  2:314446/1‑149 (MQ=255)                                                                     gCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCg                                                                               >  2:319300/1‑149 (MQ=255)                                                                      cACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCgg                                                                              >  8:257323/1‑149 (MQ=255)                                                                       aCCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGc                                                                             <  3:64904/149‑1 (MQ=255)                                                                          tgctgGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAg                                                                          >  3:90802/1‑149 (MQ=255)                                                                                        taaacggcggGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCgg                                                             <  2:294721/144‑1 (MQ=255)                                                                                        gAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCgg                                                             <  1:314446/148‑1 (MQ=255)                                                                                        gAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCgg                                                             <  1:319300/148‑1 (MQ=255)                                                                                         aaTCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGc                                                                                                                                       <  2:240354/73‑1 (MQ=255)                                                                                         aaTCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGATGCGCGAAGTCGACAGCT‑TCGACTTCGc                                                                                                                                       >  1:240354/1‑73 (MQ=255)                                                                                                gggCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTCCCGGCAAGGCGCTGCAGGCGGGCTGCACgg                                                    <  1:27615/149‑1 (MQ=255)                                                                                                 ggCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGt                                                   <  3:236383/149‑1 (MQ=255)                                                                                                    agcagcGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCAt                                                >  7:140082/1‑149 (MQ=255)                                                                                                       agcGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGa                                             <  2:118925/149‑1 (MQ=255)                                                                                                       agcGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGa                                             <  3:128258/149‑1 (MQ=255)                                                                                                        gcGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGAc                                            <  1:307219/149‑1 (MQ=255)                                                                                                        gcGTGATGTTCGCTAGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGAc                                            <  4:326181/149‑1 (MQ=255)                                                                                                          gTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCt                                          <  1:34374/149‑1 (MQ=255)                                                                                                              tGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTcc                                      >  7:283876/1‑149 (MQ=255)                                                                                                                    cTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCggcggc                                <  5:4118/149‑1 (MQ=255)                                                                                                                        cAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGt                                        >  1:91604/1‑137 (MQ=255)                                                                                                                        cAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGt                                        <  2:91604/137‑1 (MQ=255)                                                                                                                         aaGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGga                           <  8:301439/149‑1 (MQ=255)                                                                                                                         aaGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGga                           >  5:255270/1‑149 (MQ=255)                                                                                                                         aaGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGga                           >  4:181576/1‑149 (MQ=255)                                                                                                                                tGAAGGTGCGCGAAGTCGACAGCT‑TCGACATCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgcc                    >  1:197978/1‑149 (MQ=255)                                                                                                                                        ggcgAAGTCGA‑AGCTGTCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCCCCGGCGCCGCCGTCAGCCGCAGCTTCGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccctggccct            >  4:243535/2‑149 (MQ=255)                                                                                                                                        cgcgAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccctggccc             <  2:197978/148‑1 (MQ=255)                                                                                                                                         gcgAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccctggccctg           >  5:21049/1‑149 (MQ=255)                                                                                                                                         gcgAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACgccctggccctg           <  8:283876/149‑1 (MQ=255)                                                                                                                                           gAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGt         <  5:32732/149‑1 (MQ=255)                                                                                                                                                  aCAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGaa  <  6:284351/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     GCCACCGGCAGCGTCGGTGAAGCCTTGGTGCAGATCCTCGAAGAACTGGCATTCCCGGTTGCCACTCTGCACCTGCTGGCCAGCATGGAATCGGCGGGCAGCAGCGTGATGTTCGCTGGCAAGAAACTGAAGGTGCGCGAAGTCGACAGCT‑TCGACTTCGCCCAGGTCAAGCTGGCCTTCTTCGCCACCGGCGCCGCCGTCAGCCGCAGCTTTGCCGGCAAGGCGCTGCAGGCGGGCTGCACGGTCATCGACCTGTCCGGCGGCCTGGACGACGCCCTGGCCCTGGTGCCGGAA  >  NC_002947/2258370‑2258663

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