Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,920,270	A→G	11.5%	L499P (CTG→CCG)	acs ←	acyl‑CoA synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,920,270	0	A	G	11.5%	154.0 / 10.8	61	L499P (CTG→CCG)	acs	acyl‑CoA synthetase
Reads supporting (aligned to +/- strand):  ref base A (26/28);  new base G (3/4);  total (29/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.82e-01

CGCATCCGGAATTTCTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCATGT  >  NC_002947/3920126‑3920411                                                                                                                                                 |                                                                                                                                               cgcATCCGGAATTTCTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGAACATCTACATCACTGACAGCATTCCCCGGGTGCAGCCg                                                                                                                                            >  7:333870/1‑149 (MQ=255)    aTCCGGAATTTCTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACCGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCAACCTCACTCACAGCATGCCCCCAGTGCAGCCGCAc                                                                                                                                         >  4:241502/1‑149 (MQ=255)         gAATTTCTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAACACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAgg                                                                                                                                    >  6:272328/1‑149 (MQ=255)         gAATTTCTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCGGCTGCACCCggg                                                                                                                                    <  5:315253/149‑1 (MQ=255)           aTTTCTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCTGGGTGCGGCTGCACCCgggg                                                                                                                                   >  4:265059/1‑147 (MQ=255)           aTTTCTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGc                                                                                                                                   >  3:148557/1‑148 (MQ=255)              tCTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCTGGGTGCGGCTGCACCCGGGGCat                                                                                                                                <  3:265059/148‑2 (MQ=255)               cTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGa                                                                                                                              >  5:131945/1‑149 (MQ=255)                    ccTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATctct                                                                                                                         >  4:82123/1‑149 (MQ=255)                     cTTGGCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCtt                                                                                                                        <  4:102690/149‑1 (MQ=255)                     cTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCtt                                                                                                                        >  7:83235/1‑149 (MQ=255)                          ccgtcaccgtcaTCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCg                                                                                                                   >  5:388783/1‑149 (MQ=255)                          ccgtcaccgtcaTCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCg                                                                                                                   >  8:221614/1‑149 (MQ=255)                                ccgtcaTCGCGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTtg                                                                                                             <  8:223198/149‑1 (MQ=255)                                 cgtcaTCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTtgc                                                                                                            >  8:2548/1‑149 (MQ=255)                                     aTCTGGAACTAGTCGACAAAACGGAAATACCGGTTAAACATGAAGTGCGCAATCCTTGCCCTCGCCCACTCACTCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTtgctgc                                                                                                         <  8:333870/148‑1 (MQ=255)                                       cGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAg‑g                                                                                                     >  8:107668/1‑149 (MQ=255)                                          gAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCa                                                                                                  <  3:82123/149‑1 (MQ=255)                                             cTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGc                                                                                               <  6:137641/149‑1 (MQ=255)                                              tcgtcgACAAAACGGAAATACCGGGACACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGc                                                                                               >  1:85642/1‑148 (MQ=255)                                                gtcgACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCa                                                                                             >  2:284486/1‑148 (MQ=255)                                                gtcgACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCa                                                                                             >  6:261146/1‑148 (MQ=255)                                                    aCAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCGGCTGCACCCggggc                                                                                                                                  >  2:131853/1‑106 (MQ=255)                                                    aCAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCGGCTGCACCCggggc                                                                                                                                  <  1:131853/108‑3 (MQ=255)                                                       aaaCGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTg                                                                                     <  7:221614/149‑1 (MQ=255)                                                           ggAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGcac                                                                                  >  5:181632/1‑148 (MQ=255)                                                           ggAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACCCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGcaca                                                                                 >  6:185599/1‑149 (MQ=255)                                                                   cGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTtgc                                                                                                            <  7:25040/115‑1 (MQ=255)                                                                   cGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTtgc                                                                                                            >  8:25040/1‑115 (MQ=255)                                                                           cTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCCCCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGACCCCGCCgggt                                                                 >  2:295698/1‑149 (MQ=255)                                                                           cTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCgggt                                                                 <  1:4887/149‑1 (MQ=255)                                                                                    cgcAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTgaag                                                         >  1:191536/1‑148 (MQ=255)                                                                                     gcAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCCCCTGCACATCGGCCACCGCCGGGTGGGTgaagaa                                                       <  3:357530/149‑1 (MQ=255)                                                                                     gcAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTgaagaa                                                       <  5:243615/149‑1 (MQ=255)                                                                                      cAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTgaagaag                                                      <  6:181632/149‑1 (MQ=255)                                                                                        aTCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTgaagaagaa                                                    <  2:191536/149‑1 (MQ=255)                                                                                         tCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAGCGGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTgaagaagaa                                                    <  1:284486/149‑1 (MQ=255)                                                                                         tCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAAc                                                   <  4:265845/149‑1 (MQ=255)                                                                                               cccTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAActcctc                                              <  8:83235/148‑1 (MQ=255)                                                                                                ccTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCTATCACCTGCACATCGGCCACCGCCgggtgg                                                               >  6:136848/1‑130 (MQ=255)                                                                                                ccTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCTATCACCTGCACATCGGCCACCGCCgggtgg                                                               <  5:136848/130‑1 (MQ=255)                                                                                                    gCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGc                                          <  7:107668/147‑1 (MQ=255)                                                                                                    gCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGc                                          <  7:2548/147‑1 (MQ=255)                                                                                                            aCGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCGACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCggg                                 <  5:272328/148‑1 (MQ=255)                                                                                                            aCGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCgg                                  <  5:261146/147‑1 (MQ=255)                                                                                                              gCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTa                               >  6:158109/1‑148 (MQ=255)                                                                                                                aaCTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCAc                                                                         <  4:208835/104‑1 (MQ=255)                                                                                                                aaCTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCAc                                                                         >  3:208835/1‑104 (MQ=255)                                                                                                                aaCTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCa                                            >  4:369441/1‑133 (MQ=255)                                                                                                                aaCTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCa                                            <  3:369441/133‑1 (MQ=255)                                                                                                                       ccTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTtctc                      <  6:131945/148‑1 (MQ=255)                                                                                                                       ccTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTcg                     <  3:241502/149‑1 (MQ=255)                                                                                                                           aCTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTcgccgc                 >  7:57097/1‑149 (MQ=255)                                                                                                                            cTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTcgccgcc                <  2:85642/149‑1 (MQ=255)                                                                                                                               aCAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTcgccgccgc              >  3:349218/1‑148 (MQ=255)                                                                                                                                cAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGa            <  6:388783/149‑1 (MQ=255)                                                                                                                                       gcctgggTGCGGCTGCACCCGGGGCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCa     >  5:59310/5‑149 (MQ=21)                                                                                                                                       ccccGGGTGCAGCCGCACCCAGGCGTCGAGCTCTTCACCGTCCTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCa     <  1:295698/149‑1 (MQ=255)                                                                                                                                        cctgggTGCGGCTGCACCCGGGGCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCAt    <  6:59310/146‑1 (MQ=21)                                                                                                                                        cccGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCAt    <  5:185599/149‑1 (MQ=255)                                                                                                                                          cGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCATGt  >  7:223173/1‑149 (MQ=255)                                                                                                                                                 aGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAAt                            <  6:157272/116‑1 (MQ=255)                                                                                                                                                 aGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAAt                            >  5:157272/1‑116 (MQ=255)                                                                                                                                                 |                                                                                                                                               CGCATCCGGAATTTCTGCACCTTGCCCGTCACCGTCATCGGGAACTCGTCGACAAAACGGAAATACCGGGGCACCTTGAAGTGCGCAATCCTTGCCCTCGCCCACTCACGCAACTCTACCTCACTGACAGCATGCCCCGGGTGCAGCCGCACCCAGGCCACGATCTCTTCACCGTACTTGCTGCAG‑GGCACGCCAATCACCTGCACATCGGCCACCGCCGGGTGGGTGAAGAAGAACTCCTCCAGCTCGCGCGGGTAAATGTTCTCGCCGCCGCGAATGATCATGT  >  NC_002947/3920126‑3920411

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