Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,284,645	G→C	9.0%	G17A (GGC→GCC)	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,645	0	G	C	9.0%	192.5 / 9.0	67	G17A (GGC→GCC)	hbd	3‑hydroxybutyryl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (28/31);  new base C (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 = 6.29e-01

AGCGCATCCAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGCACTCTCGACGCCGACAAGGCCGCAGC  >  NC_002947/4284500‑4284792                                                                                                                                                   |                                                                                                                                                     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                                                                       <  8:104309/149‑1 (MQ=255)                                                                            gcgTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGgtgt                                                                       <  2:292102/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                                                                                          >  6:22796/1‑120 (MQ=255)                                                                                      gaggCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACg                                                                                          <  5:22796/120‑1 (MQ=255)                                                                                       aggCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGtgtg                                                                      >  6:70676/1‑139 (MQ=255)                                                                                       aggCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGtgtg                                                                      <  5:70676/139‑1 (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                                             >  2:156100/1‑147 (MQ=255)                                                                                                        ttGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAgc                                             <  1:156100/147‑1 (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)                                                                                                                                                  ggCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGCACTCTCGACGCCGACAAGGCCGCAgc  <  1:57006/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     AGCGCATCCAGGGCCGCTACGCCCTGGCCACGATGTGTATCGGCGGTGGCCAAGGCATCGCCGTCGTGTTCGAGCGCGTCTGAGGGAGGCTGACACATGAGCATTGAAAAGATCGCCGTGATCGGCGCGGGCAC‑CATGGGCAACGGCATTGCCCAGGTGTGCGCCATTGCCGGCTACCAGGTGCT‑GCTGGTGGATGTTTCCGACGCAGCGCTCGAGCGCGGTGTGGCCACCTTGAGCAAGAACCTCGAGCGCCAGGTCAGCAAGGGCACTCTCGACGCCGACAAGGCCGCAGC  >  NC_002947/4284500‑4284792

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