Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,284,643	C→T	8.9%	N16N (AAC→AAT)	hbd →	3‑hydroxybutyryl‑CoA dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,284,643	0	C	T	8.9%	188.0 / 6.3	67	N16N (AAC→AAT)	hbd	3‑hydroxybutyryl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (29/32);  new base T (2/4);  total (31/36)
Fisher's exact test for biased strand distribution p-value = 6.78e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.86e-01

AGCGCATCCAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGCACTCTCGACGCCGACAAGG  >  NC_002947/4284500‑4284785                                                                                                                                                 |                                                                                                                                                aGCGCATCCAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCAt                                                                                                                                            >  8:32084/1‑149 (MQ=255)  gcgcATCCAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCAtt                                                                                                                                           >  4:137278/1‑149 (MQ=255)       tCCAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAATGCCGTTGCCCa                                                                                                                                      <  7:100104/149‑1 (MQ=255)         cAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAATGCCGTTGCCCatg                                                                                                                                    <  8:83306/149‑3 (MQ=255)          aGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCACACCTGGGCAATGCCGTTGCCCatg                                                                                                                                    <  1:154752/149‑3 (MQ=255)                gCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTgcgc                                                                                                                             >  2:87673/1‑149 (MQ=255)                      cccTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGc                                                                                                                       >  1:308402/1‑149 (MQ=255)                       ccTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGcc                                                                                                                      >  3:149036/1‑149 (MQ=255)                             cACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTAc                                                                                                                >  7:104309/1‑149 (MQ=255)                                    tgtATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGtgc                                                                                                         >  1:292102/1‑149 (MQ=255)                                        tCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGtgct‑gc                                                                                                     <  5:67942/148‑1 (MQ=255)                                        tCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGtgct‑gc                                                                                                     <  7:32084/148‑1 (MQ=255)                                           gcggTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCtggt                                                                                                 <  5:196641/149‑1 (MQ=255)                                           gcggTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCtggt                                                                                                 <  4:149036/149‑1 (MQ=255)                                            cggTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCtggtg                                                                                                <  3:76005/149‑1 (MQ=255)                                                 gCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCg                                                                                                                     <  2:100925/124‑1 (MQ=255)                                                 gCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCg                                                                                                                     >  1:100925/1‑124 (MQ=255)                                                  ccAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAATGCCGTTGCCCatg                                                                                                                                    <  5:45855/108‑3 (MQ=255)                                                  ccAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAATGCCGTTGCCCatg                                                                                                                                    >  6:45855/1‑106 (MQ=255)                                                   cAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGttt                                                                                         >  3:39403/1‑149 (MQ=255)                                                   cAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGtt                                                                                          >  4:32978/1‑148 (MQ=255)                                                               gtcgtGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAgcgc                                                                             <  3:137278/149‑1 (MQ=255)                                                                  gtgtTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTc                                                                           <  2:4444/148‑1 (MQ=255)                                                                     ttCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAgc                                                                       >  5:289263/1‑149 (MQ=255)                                                                          gcgcGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGgt                                                                  <  3:32978/149‑1 (MQ=255)                                                                            gcgTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGgtgt                                                                <  2:292102/149‑1 (MQ=255)                                                                            gcgTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGgtgt                                                                <  8:104309/149‑1 (MQ=255)                                                                              gTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTgg                                                              >  2:185402/1‑149 (MQ=255)                                                                              gTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTgg                                                              <  1:185402/149‑1 (MQ=255)                                                                               tCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGc                                                             <  2:153878/149‑1 (MQ=255)                                                                                      gaggCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACg                                                                                   <  5:22796/120‑1 (MQ=255)                                                                                      gaggCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACg                                                                                   >  6:22796/1‑120 (MQ=255)                                                                                       aggCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGtgtg                                                               <  5:70676/139‑1 (MQ=255)                                                                                       aggCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGtgtg                                                               >  6:70676/1‑139 (MQ=255)                                                                                       aggCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGa                                                     <  4:39403/149‑1 (MQ=255)                                                                                        ggCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCTCGCTGGTGGATGTTTCCGACGCAGTGCTCGAGCGCGGTGTGGCCACCTTGa                                                     <  1:87673/149‑1 (MQ=255)                                                                                           tGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGtgc                                                                                                         <  8:174257/94‑1 (MQ=255)                                                                                           tGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAgcg                                                                      >  7:274280/1‑128 (MQ=255)                                                                                           tGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAgcg                                                                      <  8:274280/128‑1 (MQ=255)                                                                                               acaTGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGa                                                     >  4:309152/1‑141 (MQ=255)                                                                                               acaTGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGa                                                     <  3:309152/141‑1 (MQ=255)                                                                                                        ttGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAgcgc                                    >  6:86513/1‑149 (MQ=255)                                                                                                        ttGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAgc                                      <  1:156100/147‑1 (MQ=255)                                                                                                        ttGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAgc                                      >  2:156100/1‑147 (MQ=255)                                                                                                         tGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGACATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAgcgc                                    >  1:188516/1‑148 (MQ=255)                                                                                                           aaaaGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCa                                  <  2:308402/148‑1 (MQ=255)                                                                                                           aaaaGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAg                                 >  5:276895/1‑149 (MQ=255)                                                                                                             aaGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGACATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGt                               <  2:188516/149‑1 (MQ=255)                                                                                                                aTCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAg                            >  2:178425/1‑149 (MQ=255)                                                                                                                   gCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGc                           >  6:315734/1‑147 (MQ=255)                                                                                                                    ccGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTc                                                                           <  4:207996/98‑1 (MQ=255)                                                                                                                     cGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTTCCCAGGTGTGCGCCATTGGCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACTCAGCGCTCGAGCGCGTTGTGGCCACCTTTAGCAAGAACCTCGAGCGCCAGGTCAGCAAgg                       <  4:327302/149‑1 (MQ=255)                                                                                                                     cGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAgg                       >  5:267002/1‑149 (MQ=255)                                                                                                                     cGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAgg                       >  8:276952/1‑149 (MQ=255)                                                                                                                            ggCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGCActct                <  6:276895/149‑1 (MQ=255)                                                                                                                                cgGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGCACTCTCGAc            <  7:276952/149‑1 (MQ=255)                                                                                                                                 gggCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGCACTCTCGACg           >  4:88896/1‑149 (MQ=255)                                                                                                                                  ggCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGCACTCTCGACGc          <  5:86513/149‑1 (MQ=255)                                                                                                                                      c‑cATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGCACTCTCGACGCCGa       >  2:57006/1‑148 (MQ=255)                                                                                                                                      c‑cATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGCACTCTCGACGCCGAc      <  6:289263/149‑1 (MQ=255)                                                                                                                                         ctGGGCAATGCCGTTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGCACTCTCGACGCCGAc      >  3:88115/2‑147 (MQ=255)                                                                                                                                         aTGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAg                        >  2:122976/1‑129 (MQ=255)                                                                                                                                         aTGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAg                        <  1:122976/129‑1 (MQ=255)                                                                                                                                         aTGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGCAct                  <  6:249849/135‑1 (MQ=255)                                                                                                                                         aTGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGCAct                  >  5:249849/1‑135 (MQ=255)                                                                                                                                         aTGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGCACTCTCGACGCCGAc      >  7:133419/1‑147 (MQ=255)                                                                                                                                            ggCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGCACTCTCGACGCCGACAAgg  <  1:178425/148‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                AGCGCATCCAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGCACTCTCGACGCCGACAAGG  >  NC_002947/4284500‑4284785

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