Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,794,649	A→T	5.7%	L602Q (CTG→CAG)	fcs ←	feruloyl‑CoA‑synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,794,649	0	A	T	5.7%	219.9 / 2.6	70	L602Q (CTG→CAG)	fcs	feruloyl‑CoA‑synthetase
Reads supporting (aligned to +/- strand):  ref base A (33/33);  new base T (2/2);  total (35/35)
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.19e-01

GACCCAGGGCGTGCGCACCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGC  >  NC_002947/3794502‑3794788                                                                                                                                                     |                                                                                                                                             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                                                                                                         <  6:227416/149‑1 (MQ=255)                                    tgtCCTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGa                                                                                                         <  3:24711/149‑1 (MQ=255)                                       ccTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGAttt                                                                                                      >  7:378852/1‑149 (MQ=255)                                        cTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑Cgg                                                                                                                                       >  7:258297/1‑115 (MQ=255)                                        cTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑Cgg                                                                                                                                       <  8:258297/115‑1 (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                                                                                                         >  6:321406/1‑127 (MQ=255)                                                          cccGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGa                                                                                                         <  5:321406/127‑1 (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                       >  5:2341/1‑148 (MQ=255)                                                                                                                        gcgcCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGtt                       >  1:164337/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              <  4:148739/149‑1 (MQ=255)                                                                                                                                aCCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTc              <  5:183510/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)                                                                                                                                                     |                                                                                                                                             GACCCAGGGCGTGCGCACCATCTGGACGATGGCCGTGTCCTCGAAGACGGTGCACAGCCCGGCCGC‑ACTCAAGGCCGCACCTTGGCGTGCAATGCGTCGGGGTCTTCGCCACGGTACAGCGCCTCGACCTGAGCGGCGCGCCGCTGCAGCA‑CGGCGCGCTGATTGATCGAGCCCTTGTCGGTGATTTCACCGGCGTCGATCGACGGCGGCTCGGCCAGCAGCGACAGCCATTCGATACGGCTGGCGTTGCCTTGGGCATCGCGGTTCAAGCGCTCCAGCCAGTCAGC  >  NC_002947/3794502‑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: ‑