Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,794,645	Δ1 bp	5.3%	coding (1809/1884 nt)	fcs ←	feruloyl‑CoA‑synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,794,645	0	C	.	5.3%	250.1 / 6.5	75	coding (1809/1884 nt)	fcs	feruloyl‑CoA‑synthetase
Reads supporting (aligned to +/- strand):  ref base C (36/35);  new base . (2/2);  total (38/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

TCGACCCAGGGCGTGCGCACCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGC  >  NC_002947/3794500‑3794788                                                                                                                                                   |                                                                                                                                                 tCGACCCAGGGCGTGCGCACCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTgc                                                                                                                                               <  3:178185/149‑1 (MQ=255) tCGACCCAGGGCGTGCGCACCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTgc                                                                                                                                               <  3:172064/149‑1 (MQ=255)   gACCCAGGGCGTGCGCACCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTgcag                                                                                                                                             >  1:151134/1‑149 (MQ=255)   gACCCAGGGCGTGCGCACCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTgca                                                                                                                                              <  6:324681/148‑1 (MQ=255)     cccAGGGCGTGCGCACCATCTGTACGATGGCCGTGTCCTCGAAGACCGTGCACAGCCCGGCCGC‑ACTCAACGCCGCACCTTGGCGTGAAATGAGTCGGTGTCTTCTCCACGGTACATCGCCTCTACCTGAGCGGCGCGCCGCTccagca                                                                                                                                           >  6:374679/1‑149 (MQ=255)     cccAGGGCGTGCGCACCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTgcagca                                                                                                                                           >  4:116866/1‑149 (MQ=255)      ccAGGGCGTGCGCACCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑c                                                                                                                                         >  5:227416/1‑149 (MQ=255)       cAGGGCGTGCGCACCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑Cg                                                                                                                                        <  3:116866/149‑1 (MQ=255)             gTGCGCACCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGgcgcgc                                                                                                                                  <  6:365998/149‑1 (MQ=255)              tGCGCACCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCt                                                                                                                                 <  5:149592/149‑1 (MQ=255)                         tGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAg                                                                                                                      >  6:231632/1‑149 (MQ=255)                             cGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCt                                                                                                                  >  3:148739/1‑149 (MQ=255)                              gATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCtt                                                                                                                 >  5:123729/1‑149 (MQ=255)                                tGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGt                                                                                                               <  6:123729/149‑1 (MQ=255)                                  gCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCg                                                                                                             >  4:24711/1‑149 (MQ=255)                                      tgtCCTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGa                                                                                                         <  3:24711/149‑1 (MQ=255)                                      tgtCCTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGa                                                                                                         <  6:227416/149‑1 (MQ=255)                                         ccTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGAttt                                                                                                      >  7:378852/1‑149 (MQ=255)                                          cTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑Cgg                                                                                                                                       <  8:258297/115‑1 (MQ=255)                                          cTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑Cgg                                                                                                                                       >  7:258297/1‑115 (MQ=255)                                          cTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATc                                                                                                                         <  3:292080/129‑1 (MQ=255)                                          cTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATc                                                                                                                         >  4:292080/1‑129 (MQ=255)                                           tCGAAGACGGTGCACAGCCCGGCCGCGACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGt                                                                                                           >  1:213628/1‑143 (MQ=255)                                           tCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGt                                                                                                           <  2:213628/142‑1 (MQ=255)                                           tCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCa                                                                                                    >  6:95493/1‑149 (MQ=255)                                             gAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCAc                                                                                                   <  4:149298/148‑1 (MQ=255)                                              aaGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCg                                                                                                 <  6:407583/149‑1 (MQ=255)                                                      gCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGtcgat                                                                                         >  1:64311/1‑149 (MQ=255)                                                       cacaGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGtcgatc                                                                                        <  2:64311/149‑1 (MQ=255)                                                        acaGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGCCCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGtcgatc                                                                                        >  6:91043/1‑148 (MQ=255)                                                        acaGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGtcgat                                                                                         >  6:311653/1‑147 (MQ=255)                                                        acaGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGtcgat                                                                                         <  5:311653/147‑1 (MQ=255)                                                            cccGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGgcgcgc                                                                                                                                  >  1:392849/1‑102 (MQ=255)                                                            cccGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGgcgcgc                                                                                                                                  <  2:392849/102‑1 (MQ=255)                                                            cccGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGa                                                                                                         <  5:321406/127‑1 (MQ=255)                                                            cccGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGa                                                                                                         >  6:321406/1‑127 (MQ=255)                                                                  cGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTc                                                                             <  3:30006/149‑1 (MQ=255)                                                                    c‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCgg                                                                           >  6:191805/1‑149 (MQ=255)                                                                        tCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGc                                                                          >  4:250740/1‑147 (MQ=255)                                                                        tCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGc                                                                          >  7:316409/1‑147 (MQ=255)                                                                        tCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCca                                                                        >  8:73107/1‑149 (MQ=255)                                                                              ccGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCg                                                                  >  4:121507/1‑149 (MQ=255)                                                                                gCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGa                                                                 >  3:374688/1‑148 (MQ=255)                                                                                  aCCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCTGCGtcgatc                                                                                        <  6:84328/123‑1 (MQ=255)                                                                                  aCCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAg                                                              >  8:110732/1‑149 (MQ=255)                                                                                  aCCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGAGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGtcgatc                                                                                        >  5:84328/1‑123 (MQ=255)                                                                                    cTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGc                                                             >  6:236368/1‑148 (MQ=255)                                                                                        gCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTc                                                        <  5:95493/149‑1 (MQ=255)                                                                                                  cGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTg                                              <  6:61095/149‑1 (MQ=255)                                                                                                   gTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTgg                                             <  7:42667/149‑1 (MQ=255)                                                                                                    tCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGc                                            >  5:286390/1‑149 (MQ=255)                                                                                                    tCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGc                                            <  3:121507/149‑1 (MQ=255)                                                                                                       gggTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGtt                                         <  6:286390/149‑1 (MQ=255)                                                                                                        ggTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTg                                        <  5:91043/149‑1 (MQ=255)                                                                                                          tCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGcc                                      >  5:300390/1‑149 (MQ=255)                                                                                                            ttCGCCACGGTGCAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCtt                                    >  7:238718/1‑149 (MQ=255)                                                                                                              cGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTTATAGAGCCCTTGTCGGTGATTTCACCGGCGTCGATAGACGGCGGCTCGGCCAGCAGCCACAGCCATTCGATACTGCTGGCGTTGCATTg                                   >  4:86878/1‑148 (MQ=255)                                                                                                              cGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCtt                                    <  3:86878/147‑1 (MQ=255)                                                                                                                    ggTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATcg                            <  1:160716/149‑1 (MQ=255)                                                                                                                        cAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGt                        <  8:238718/149‑1 (MQ=255)                                                                                                                          gcgcCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGtt                       >  1:164337/1‑148 (MQ=255)                                                                                                                          gcgcCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGtt                       >  5:2341/1‑148 (MQ=255)                                                                                                                            gcCTCGCCCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCaa                    <  5:192986/149‑1 (MQ=255)                                                                                                                            gcCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCaa                    <  5:236368/149‑1 (MQ=255)                                                                                                                             ccTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGc                                                                          <  3:90330/94‑1 (MQ=255)                                                                                                                             ccTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGc                                                                          >  4:90330/1‑94 (MQ=255)                                                                                                                                  aCCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTc              <  5:183510/149‑1 (MQ=255)                                                                                                                                  aCCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTc              <  4:148739/149‑1 (MQ=255)                                                                                                                                   ccTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTc              >  2:46098/1‑148 (MQ=255)                                                                                                                                       atcagcgcgcCG‑TGCTGCAGCGGCGCGCCGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCgg                         <  8:347430/129‑1 (MQ=255)                                                                                                                                        tcagcgcgcCG‑TGCTGCAGCGGCGCGCCGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGAc                                                                                     >  7:316168/4‑72 (MQ=21)                                                                                                                                        tcagcgcgcCG‑TGCTGCAGCGGCGCGCCGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGAc                                                                                     <  8:316168/69‑1 (MQ=21)                                                                                                                                        tcagcgcgcCG‑TGCTGCAGCGGCGCGCCGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTccagc          >  8:343967/4‑147 (MQ=255)                                                                                                                                        gcggcgCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTccagcca        <  8:378852/149‑1 (MQ=255)                                                                                                                                             gcgcCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAg   >  8:123958/1‑149 (MQ=255)                                                                                                                                              cgccgcTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGc  <  4:374688/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                 TCGACCCAGGGCGTGCGCACCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGC  >  NC_002947/3794500‑3794788

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