Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,962,512	T→A	6.2%	L237Q (CTG→CAG)	PP_3492 →	acyl‑CoA dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,962,512	0	T	A	6.2%	193.2 / 2.8	64	L237Q (CTG→CAG)	PP_3492	acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (28/32);  new base A (2/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.75e-01

CCACCGACAACCCTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATT  >  NC_002947/3962365‑3962648                                                                                                                                                      |                                                                                                                                           ccACCGACAACCCTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTg                                                                                                                                           >  6:5409/1‑149 (MQ=255)  cACCGACTACCCTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGt                                                                                                                                          <  2:126071/149‑1 (MQ=255)  cACCGACAACCCTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGt                                                                                                                                          <  6:379535/149‑1 (MQ=255)   aCCGACAACCCTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTc                                                                                                                                         >  5:194718/1‑149 (MQ=255)       aCAACCCTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTc                                                                                                                                         >  1:164959/1‑145 (MQ=255)       aCAACCCTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTc                                                                                                                                         <  2:164959/145‑1 (MQ=255)        cAACCCTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAAcc                                                                                                                                    <  7:260838/149‑1 (MQ=255)          aCCCTGGGTTCTAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCAGATACCTGCGCGGTCACCTTCGACAACTCCGCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAAcc                                                                                                                                    <  1:111039/149‑1 (MQ=255)          aCCCTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTc                                                                                                                                  <  6:28354/149‑1 (MQ=255)           cccTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGGCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTc                                                                                                                                  >  1:302924/1‑148 (MQ=255)           cccTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTc                                                                                                                                         >  8:78194/1‑141 (MQ=255)           cccTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTc                                                                                                                                         <  7:78194/141‑1 (MQ=255)            ccTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCACCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGa                                                                                                                                <  3:183315/149‑1 (MQ=255)             cTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATAACTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCg                                                                                                                                 >  1:282192/1‑147 (MQ=255)              tGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTCCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATACTCGGTGAGCGTGGCAAGGGCCTCGCCATTGCCCTGTCCAACCTCGaat                                                                                                                              >  4:5263/1‑148 (MQ=255)                ggTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAg                                                                                                                              <  5:5409/147‑1 (MQ=255)                        gTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCAGTCCAACCTCGAAGGCGGCCGTAt                                                                                                                    >  7:159675/1‑149 (MQ=255)                             ccGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCa                                                                                                              >  2:282408/1‑149 (MQ=255)                                cAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTg                                                                                                           <  2:208080/149‑1 (MQ=255)                                       cccatgATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCa                                                                                                    <  3:5263/144‑1 (MQ=255)                                          aaGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGc                                                                                                 >  2:245640/1‑149 (MQ=255)                                          aaGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGc                                                                                                 >  8:308966/1‑149 (MQ=255)                                            gATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCt                                                                                               >  1:185095/1‑149 (MQ=255)                                                ggCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTgg                                                                                            <  2:185095/148‑1 (MQ=255)                                                    ttCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATc                                                                                       <  2:302924/149‑1 (MQ=255)                                                       gcgcCTCCGATATCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATc                                                                                                                   <  7:255044/119‑1 (MQ=255)                                                       gcgcCTCCGATATCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATc                                                                                                                   >  8:255044/1‑119 (MQ=255)                                                       gagcCTCCGATATCTGCGCGGTCACCTTCGAAAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCTGCCGTATc                                                                                                                   >  2:254062/3‑119 (MQ=255)                                                       gagcCTACGATATCTGCGCTGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCGTTGCCATTGCCCTGTCCAACCTAGAATGCGGCCGTATc                                                                                                                   <  1:254062/117‑1 (MQ=255)                                                        cgcCTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGccc                                                                                   <  2:282192/149‑1 (MQ=255)                                                          ccTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCcgc                                                                                 >  5:363619/1‑149 (MQ=255)                                                            tCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGc                                                                                     <  8:126848/143‑1 (MQ=255)                                                            tCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGc                                                                                     >  7:126848/1‑143 (MQ=255)                                                              cGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGcgccg                                                                             >  7:131617/1‑149 (MQ=255)                                                              cGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGcgcc                                                                              <  7:222168/148‑1 (MQ=255)                                                                 tACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCAt                                                                          <  6:194718/149‑1 (MQ=255)                                                                    cTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCGGCATTc                                                                        >  6:283057/1‑148 (MQ=255)                                                                       cgcgGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGaa                                                                     <  4:67645/148‑1 (MQ=255)                                                                         cgGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑gcg                                                                 <  5:174036/149‑1 (MQ=255)                                                                              aCCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑gcggcg                                                              <  6:283939/147‑1 (MQ=255)                                                                              aCCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑gcggcg                                                              >  5:283939/1‑147 (MQ=255)                                                                              aCCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGgcgc                                                             <  1:245640/148‑1 (MQ=255)                                                                               ccTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGC‑TTCGAATGCGgcgcgg                                                           <  7:107259/149‑3 (MQ=255)                                                                                cTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATCTGGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTg                                                           <  1:282408/149‑1 (MQ=255)                                                                                         cT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTcgc                                                 >  6:360144/1‑149 (MQ=255)                                                                                                gCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGACTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGAc                                            <  1:92367/149‑1 (MQ=255)                                                                                                   ttCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCAGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCg                                          >  6:243718/1‑148 (MQ=255)                                                                                                           gccAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCAGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCATACAGtt                                 <  8:159675/149‑1 (MQ=255)                                                                                                                  tCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCaa                           >  3:272400/1‑148 (MQ=255)                                                                                                                  tCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCaa                           <  4:272400/148‑1 (MQ=255)                                                                                                                  tCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCAAg                          >  5:22082/1‑149 (MQ=255)                                                                                                                    cTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCAAGc                         <  7:308966/148‑1 (MQ=255)                                                                                                                       ggTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCTGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCGAGCCGAt                     >  8:56281/1‑149 (MQ=255)                                                                                                                           aGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCCTTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCAAGCCGATCAAc                 >  3:78237/1‑149 (MQ=255)                                                                                                                              gTGGCAAGGGCCTGGCCATTGCCCAGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCAAGCCGATCAACGaa              <  5:243718/149‑1 (MQ=255)                                                                                                                               tGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCAAGCCGATCAACGaa              >  3:256039/1‑148 (MQ=255)                                                                                                                                 gCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCAAGCCGATCAACGAacac           >  3:70696/1‑149 (MQ=255)                                                                                                                                 gCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCAAGCCGATCAACGAacac           >  1:323024/1‑149 (MQ=255)                                                                                                                                   aaGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCAAGCCGATCAACGAACACCa         <  8:131617/149‑1 (MQ=255)                                                                                                                                       gCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCAAGCCGATCAACGAACACCAGAGc     <  6:229937/149‑1 (MQ=255)                                                                                                                                          tGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCAAGCCGATCAACGAACACCAGAGCAtt  >  2:68217/1‑149 (MQ=255)                                                                                                                                            gCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCAAGCCGATCAACGAACACCAGAGCAtt  <  4:256039/147‑1 (MQ=255)                                                                                                                                               aTTGCCCTGTCCAATCTCGTAGGCGGCCGTATC‑GGCATTGCCGCCCAGTCCTTGGGCATCGCCCGCGCAGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCATa                                      <  6:38215/108‑1 (MQ=255)                                                                                                                                               aTTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GAGGAGCTGGTCTACTCGCGTGACCGCATa                                      >  5:38215/1‑108 (MQ=255)                                                                                                                                                      |                                                                                                                                           CCACCGACAACCCTGGGTTCAAGGTTGACCGCAGCGAACACAAGATGGGCATTCGCGCCTCCGATACCTGCGCGGTCACCTTCGACAACT‑‑GCCGCATTCCTGCCGCCAATATCCTCGGTGAGCGTGGCAAGGGCCTGGCCATTGCCCTGTCCAACCTCGAAGGCGGCCGTATC‑GGCATTGCCGCCCAGGCCTTGGGCATCGCCCGCGCCGCATTCGAA‑GCGGCGCTGGTCTACTCGCGTGACCGCATACAGTTTGGCAAGCCGATCAACGAACACCAGAGCATT  >  NC_002947/3962365‑3962648

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