Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,097,457	C→T	8.2%	P135L (CCG→CTG)	gpx →	hydroperoxy fatty acid reductase Gpx1

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,097,457	0	C	T	8.2%	169.1 / 4.1	61	P135L (CCG→CTG)	gpx	hydroperoxy fatty acid reductase Gpx1
Reads supporting (aligned to +/- strand):  ref base C (32/24);  new base T (3/2);  total (35/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.01e-01

ATTTCGGCGTGAGTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAG  >  NC_002947/2097310‑2097604                                                                                                                                                    |                                                                                                                                                     aTTTCGGCGTGAGTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCg                                                                                                                                                     <  4:252073/149‑1 (MQ=255) aTTTCGGCGTGAGTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCg                                                                                                                                                     <  4:247702/149‑1 (MQ=255)        cGTGAGTTTTCCGCTGTTTCGCAAGGTCGAGGTGAACGGCCCGGGCGCTCACCCTCTGTTCCTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGTAGATCAAGTGGAAGTTCTCCAAGTTTCTGGTCGACCCGGCCAGTg                                                                                                                                              <  5:256810/149‑1 (MQ=255)          tGAGTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGc                                                                                                                                            >  4:103242/1‑149 (MQ=255)             gTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAg                                                                                                                                         <  2:116663/149‑1 (MQ=255)              ttttCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCTGGCCACTGGCCggg                                                                                                                                        <  1:338522/149‑1 (MQ=255)                ttCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGt                                                                                                                                       <  3:217877/148‑1 (MQ=255)                 tCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAACTGGAACTTCACAAACTTTCTGGTCGACCCGGCCAGTGGCCAGGTc                                                                                                                                      >  3:109349/1‑148 (MQ=255)                             cAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGccc                                                                                                                         >  8:238876/1‑149 (MQ=255)                              aaGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGcccc                                                                                                                        >  8:46521/1‑149 (MQ=255)                                 gTCGAGGGCAACGGCCCGGGCGCTCACGCTCTGTGCGTCGACCTCAGGCACCGCGCCCCGGGTAGTTTGGGTAGGCGGAAGATCAGGTGGAGCCTCACCAAGTTTGTGGTAGACCCGGCGCGTGGCCAGGTGAAGCGCTATGccccgcg                                                                                                                     >  2:93632/1‑146 (MQ=255)                                  tCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCccac                                                                                                                     >  3:120659/1‑148 (MQ=255)                                    gAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCccaccac                                                                                                                  >  5:186111/1‑149 (MQ=255)                                       gTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCccaccacca                                                                                                                >  2:158981/1‑148 (MQ=255)                                       gTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCccaccacc                                                                                                                 >  5:117964/1‑147 (MQ=255)                                               cccGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAgg                                                                                                       <  6:298134/149‑1 (MQ=255)                                                      gcTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGa                                                                                               <  1:83271/149‑1 (MQ=255)                                                        tCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGCGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAgg                                                                                             >  1:71414/1‑149 (MQ=255)                                                          aCCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCg                                                                                           >  2:9156/1‑149 (MQ=255)                                                          aCCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCg                                                                                           <  1:158981/149‑1 (MQ=255)                                                            ccTCAGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGa                                                                                         >  4:208092/1‑149 (MQ=255)                                                                  ttCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATc                                                                                     <  7:238876/147‑1 (MQ=255)                                                                  ttCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGa                                                                                   <  6:186111/149‑1 (MQ=255)                                                                        gAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGttt                                                                             >  7:175863/1‑149 (MQ=255)                                                                            tCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGc                                                                           <  8:270738/147‑1 (MQ=255)                                                                            tCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCta                                                                         >  7:270738/1‑148 (MQ=255)                                                                              aaGCAACGCGCCCCGGGTTTTTTGGGTAGCCAGAAGATCAAGTGTAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCATGTCAAGCGCTTTGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCa                                                                        <  4:109349/148‑1 (MQ=255)                                                                              aaGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCa                                                                        <  2:71414/148‑1 (MQ=255)                                                                              aaGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCa                                                                        <  3:208092/148‑1 (MQ=255)                                                                                        ccccGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAgg                                                             >  6:121685/1‑149 (MQ=255)                                                                                        ccccGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAgg                                                             >  4:275674/1‑149 (MQ=255)                                                                                               tATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGACAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGccacc                                                       >  7:296208/1‑144 (MQ=255)                                                                                                aTTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCaa                                                     >  1:153128/1‑149 (MQ=255)                                                                                                     gggTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGTGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGt                                                 <  3:103242/148‑1 (MQ=255)                                                                                                     gggTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGt                                                 <  1:9156/148‑1 (MQ=255)                                                                                                       gTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGa                                              >  3:115546/1‑149 (MQ=255)                                                                                                       gTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGa                                              >  5:289548/1‑149 (MQ=255)                                                                                                          gCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCa                                           <  8:175863/149‑1 (MQ=255)                                                                                                          gCCAGAAGATCAAGTGGAACCTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCCACTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGAc                                             <  1:93632/148‑1 (MQ=255)                                                                                                               aagaTCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGc                                       >  2:18693/1‑148 (MQ=255)                                                                                                                  aTCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGt                                                 >  8:196764/1‑135 (MQ=255)                                                                                                                             cTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAgg                                                             <  7:284413/112‑1 (MQ=255)                                                                                                                             cTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAgg                                                             >  8:284413/1‑112 (MQ=255)                                                                                                                              ttCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCgcgg                       >  5:114434/1‑149 (MQ=255)                                                                                                                                cACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAg                                                                                                        <  5:93375/67‑1 (MQ=255)                                                                                                                                cACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAg                                                                                                        >  6:93375/1‑67 (MQ=255)                                                                                                                                 aCCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCgcggcg                     >  3:308507/1‑148 (MQ=255)                                                                                                                                   cAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGaa                  >  6:231829/1‑149 (MQ=255)                                                                                                                                   cAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGaa                  >  1:173396/1‑149 (MQ=255)                                                                                                                                   cAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGaa                  >  8:288922/1‑149 (MQ=255)                                                                                                                                         tCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGc              <  7:46521/147‑1 (MQ=255)                                                                                                                                             cttgACCTGGCCACTGGCCGGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGc         >  3:213017/4‑148 (MQ=255)                                                                                                                                              ttgACCTGGCCACTGGCCGGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGc         >  6:244035/3‑147 (MQ=255)                                                                                                                                              tCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGc         >  6:44863/1‑147 (MQ=255)                                                                                                                                              tCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGc         >  4:187170/1‑147 (MQ=255)                                                                                                                                               tgACCTGGCCACTGGCCGGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAg      >  1:57512/2‑149 (MQ=255)                                                                                                                                                gACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAgg     <  5:231829/149‑1 (MQ=255)                                                                                                                                                gACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAg      <  2:153128/148‑1 (MQ=255)                                                                                                                                                 aCCTGGCCACTGGCCGGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAgg     <  5:244035/148‑1 (MQ=255)                                                                                                                                                 aCCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTg           >  4:100144/1‑142 (MQ=255)                                                                                                                                                 aCCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTg           <  3:100144/142‑1 (MQ=255)                                                                                                                                                   ccGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAg  >  3:137251/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                     ATTTCGGCGTGAGTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAG  >  NC_002947/2097310‑2097604

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑