Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,097,463	G→C	8.1%	S137T (AGT→ACT)	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,463	0	G	C	8.1%	188.4 / 6.3	62	S137T (AGT→ACT)	gpx	hydroperoxy fatty acid reductase Gpx1
Reads supporting (aligned to +/- strand):  ref base G (34/23);  new base C (3/2);  total (37/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

CGTGAGTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAGTCGTTC  >  NC_002947/2097317‑2097610                                                                                                                                                   |                                                                                                                                                     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                                                                                                 <  1:158981/149‑1 (MQ=255)                                                   aCCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCg                                                                                                 >  2:9156/1‑149 (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                                                                              <  3:208092/148‑1 (MQ=255)                                                                       aaGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCa                                                                              <  2:71414/148‑1 (MQ=255)                                                                                 ccccGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAgg                                                                   >  4:275674/1‑149 (MQ=255)                                                                                 ccccGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAgg                                                                   >  6:121685/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                                                    >  5:289548/1‑149 (MQ=255)                                                                                                gTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGa                                                    >  3:115546/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                                                                   >  8:284413/1‑112 (MQ=255)                                                                                                                      cTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAgg                                                                   <  7:284413/112‑1 (MQ=255)                                                                                                                       ttCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCgcgg                             >  5:114434/1‑149 (MQ=255)                                                                                                                         cACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAg                                                                                                              >  6:93375/1‑67 (MQ=255)                                                                                                                         cACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAg                                                                                                              <  5:93375/67‑1 (MQ=255)                                                                                                                          aCCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCgcggcg                           >  3:308507/1‑148 (MQ=255)                                                                                                                            cAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGaa                        >  6:231829/1‑149 (MQ=255)                                                                                                                            cAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGaa                        >  8:288922/1‑149 (MQ=255)                                                                                                                            cAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGaa                        >  1:173396/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                 <  3:100144/142‑1 (MQ=255)                                                                                                                                          aCCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTg                 >  4:100144/1‑142 (MQ=255)                                                                                                                                            ccGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAg        >  3:137251/1‑149 (MQ=255)                                                                                                                                                  aGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATgg                                                        >  4:76402/1‑95 (MQ=255)                                                                                                                                                  aGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAGTCGTTc  <  6:117964/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CGTGAGTTTTCCGCTGTTTCGCAAGGTCGAGGTCAACGGCCCGGGCGCTCACCCTCTGTTCGTCGAACTCAAGCAACGCGCCCCGGGTATTTTGGGTAGCCAGAAGATCAAGTGGAACTTCACCAAGTTTCTGGTCGACCCGGCCAGTGGCCAGGTCAAGCGCTATGCCCCCACCACCAAGCCACAGGCC‑CTGGAGGCGGACATCGAGCGTTTGCTCAGCCGCTGAAGGGCACCCAATGGTCGACCAGGGCCAGTAGCTCTTCGCGGCGGAATGGCTTGGCCAGGTAGTCGTTC  >  NC_002947/2097317‑2097610

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