Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,134,274	C→T	5.9%	S252S (TCG→TCA)	PP_3638 ←	acyl‑CoA dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,134,274	0	C	T	5.9%	204.1 / 2.2	67	S252S (TCG→TCA)	PP_3638	acyl‑CoA dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (30/33);  new base T (2/2);  total (32/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 = 7.52e-01

GTCGGCAGCACTTTCAGCCAGGATGATGCGCAGCTTGAACAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTACCCGTGGCAGGAATTCGCGTGCATTTTCATGCAGCAACCAGTCGCGGGC  >  NC_002947/4134127‑4134420                                                                                                                                                    |                                                                                                                                                   gTCGGCAGCACTTTCAGCCAGGATGATGCGCAGCTTGAACAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACg                                                                                                                                                   <  4:133653/149‑1 (MQ=255)     gcagcaCTTTCAGCCAGGATGATGCGCAGCTTGAACAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTCAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGATTCGTCGAGGAACGAgg                                                                                                                                                >  1:182691/1‑148 (MQ=255)           cTTTCAGCCAGGATGATGCGCAGCTTGAACAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGcc                                                                                                                                         <  3:218436/149‑1 (MQ=255)            tttCAGCCAGGATGATGCGCAGCTTGAACAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAATGAGGCTCgg                                                                                                                                           <  5:6763/146‑1 (MQ=255)            tttCAGCCAGGATGATGCGCAGCTTGAACAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAATGAGGCTCgg                                                                                                                                           <  5:1484/146‑1 (MQ=255)            tttCAGCCAGGATGATGCGCAGCTTGAACAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAATGAGGCTCgg                                                                                                                                           >  6:1484/1‑146 (MQ=255)            tttCAGCCAGGATGATGCGCAGCTTGAACAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAATGAGGCTCgg                                                                                                                                           >  6:6763/1‑146 (MQ=255)                aGCCAGGATGATGCGCAGCTTGAACAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGc                                                                                                                                     <  2:217656/148‑1 (MQ=255)                 gCCAGGATGATGCGCAGCTTGAACAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTg                                                                                                                                      >  5:152502/1‑146 (MQ=255)                  ccAGGATGATGCGCAGCTTGAACAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGa                                                                                                                                  <  8:59491/149‑1 (MQ=255)                  ccAGGATGATGCGCAGCGTGAACAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAg                                                                                                                                                 >  5:274928/1‑134 (MQ=255)                  ccAGGATGATGCGCAGCCTGAACAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAg                                                                                                                                                 <  6:274928/134‑1 (MQ=255)                   cAGGATGATGCGCAGCTTGAACAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTc                                                                                                                                             >  3:270964/1‑137 (MQ=255)                   cAGGATGATGCGCAGCTTGAACAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTc                                                                                                                                             <  4:270964/137‑1 (MQ=255)                             cgcAGCTTGAACAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGaa                                                                                                                        <  8:199173/148‑1 (MQ=255)                              gcAGCTTGAACAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAAcc                                                                                                                      >  7:101664/1‑149 (MQ=255)                                aGCTTGAACAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACctc                                                                                                                    <  2:239094/149‑1 (MQ=255)                                  cTTGAACAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAAc                                                                                                                       <  3:114946/144‑1 (MQ=255)                                  cTTGAACAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAAc                                                                                                                       >  4:114946/1‑144 (MQ=255)                                    tGAACAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTcc                                                                                                                >  7:251510/1‑149 (MQ=255)                                     gAACAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGc                                                                                                                                               >  3:226332/1‑117 (MQ=255)                                     gAACAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGc                                                                                                                                               <  4:226332/117‑1 (MQ=255)                                      aaCAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAg                                                                                                              >  7:175109/1‑149 (MQ=255)                                        cAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAgcg                                                                                                            <  2:109609/149‑1 (MQ=255)                                         aaGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACTAGGCTCGGCCTTGc                                                                                                                                     <  5:173571/123‑1 (MQ=255)                                           ggCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCtc                                                                                                         >  5:170881/1‑149 (MQ=255)                                              cgccgGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCtt                                                                                                                  >  5:95976/1‑137 (MQ=255)                                              cgccgGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCtt                                                                                                                  <  6:95976/137‑1 (MQ=255)                                                ccgGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTcgc                                                                                                    <  8:101664/149‑1 (MQ=255)                                                       ctgctgGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAg                                                                                               <  3:137078/147‑1 (MQ=255)                                                          ctgGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAg                                                                                                                                                 >  7:258377/1‑94 (MQ=255)                                                          ctgGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAg                                                                                                                                                 <  8:258377/94‑1 (MQ=255)                                                             gAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATg                                                                                        <  6:334812/148‑1 (MQ=255)                                                                     cGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCca                                                                               >  3:27589/1‑149 (MQ=255)                                                                      gTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCca                                                                               >  3:69672/1‑148 (MQ=255)                                                                                  aTGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAGCGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCa                                                                  <  1:194303/149‑1 (MQ=255)                                                                                     cTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAgg                                                                <  4:274794/148‑1 (MQ=255)                                                                                           ttCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCg                                                         >  3:22644/1‑149 (MQ=255)                                                                                                   cGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTAcc                                                 >  5:29115/1‑149 (MQ=255)                                                                                                    ggTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCgg                                                        >  8:263120/1‑141 (MQ=255)                                                                                                    ggTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCgg                                                        <  7:263120/141‑1 (MQ=255)                                                                                                    ggTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTAccc                                                >  2:280178/1‑149 (MQ=255)                                                                                                      tgtTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTACCCGt                                              >  4:82114/1‑149 (MQ=255)                                                                                                             cAGGCTTTCGCTAAGCACCTGACGTGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTACCCGTGGCAGGa                                       <  4:229040/149‑1 (MQ=255)                                                                                                             cAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTACCCGTGGCAgg                                        <  4:22644/148‑1 (MQ=255)                                                                                                               ggCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTACCCGTGGCAGGAAt                                     >  6:70480/1‑149 (MQ=255)                                                                                                                 cGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGCTGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTACCCGTGGCAGGAATTc                                   <  8:251510/149‑1 (MQ=255)                                                                                                                 cGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTACCCGTGGCAGGAATTc                                   <  8:175109/149‑1 (MQ=255)                                                                                                                  gTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTACCCGTGGCAGGAATTcg                                  >  2:360517/1‑149 (MQ=255)                                                                                                                   ttCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTACCCGTGGCAGGAATTcgc                                 <  6:170881/149‑1 (MQ=255)                                                                                                                      gCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTACCCGTGGCAGGAATTCGCGt                               <  3:82114/148‑1 (MQ=255)                                                                                                                        tAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTACCCGTGGCAGGAATTCGCGTGCa                            >  6:344443/1‑149 (MQ=255)                                                                                                                              ccTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTACCCGTGGCAGGAATTCGCGTGCATTTTCa                      >  3:147739/1‑149 (MQ=255)                                                                                                                                  aCGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTACCCGTGGCAGGAATTcgc                                 >  8:121190/1‑134 (MQ=255)                                                                                                                                  aCGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTACCCGTGGCAGGAATTcgc                                 <  7:121190/134‑1 (MQ=255)                                                                                                                                  aCGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTACCCGTGGCAGGAATTCGCGTGCATTTTCATg                    <  6:152502/147‑1 (MQ=255)                                                                                                                                    gggCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCgg                                                        <  5:311116/109‑1 (MQ=255)                                                                                                                                    gggCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCgg                                                        >  6:311116/1‑109 (MQ=255)                                                                                                                                    gggCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTACCCGTGGCAGGAATTCGCGTGCATTTTCATgcag                 >  1:69167/1‑148 (MQ=255)                                                                                                                                            tcgAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTACCCGTGGCAGGAATTCGCGTGCATTTTCATGCAGCAACCAGTc        >  2:299721/1‑149 (MQ=255)                                                                                                                                            tcgAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTACCCGTGGCAGGAATTCGCGTGCATTTTCATGCAGCAACCAGTc        <  5:171389/149‑1 (MQ=255)                                                                                                                                              gAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTACCCGTGGCAGGAATTCGCGTGCATTTTCATGCAGCAACCAGTcgc      >  2:371666/1‑149 (MQ=255)                                                                                                                                                ggAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTACCCGTGGCAGGAATTCGCGTGCATTTTCATGTAGCAACCAGTCGCgg    <  1:280178/149‑1 (MQ=255)                                                                                                                                                 gAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTACCCGTGGCAGGAATTCGCGTGCATTTTCATGCAGCAACCAGTCGCggg   >  5:59809/1‑149 (MQ=255)                                                                                                                                                  aaCGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTACCCGTGGCAGGAATTCGCGTGCATTTTCATGCAGCAACCAGTCGCGGGc  <  4:27589/149‑1 (MQ=255)                                                                                                                                                   aCGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTACCCGTGGCAGGAATTcgc                                 <  6:348449/117‑1 (MQ=255)                                                                                                                                                   aCGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTACCCGTGGCAGGAATTcgc                                 >  5:348449/1‑117 (MQ=255)                                                                                                                                                    |                                                                                                                                                   GTCGGCAGCACTTTCAGCCAGGATGATGCGCAGCTTGAACAAGGCCGCCGGCTGCTGCTGGAAGCGACCGTCGAGCAAGCCATGCTTGAGTTCGCTGACGGTGTTTTCCAGGCGTTCGCTAAGCACCTGACGGGCTTCGTCGAGGAACGAGGCTCGGCCTTGCAGATGCGCCTGAACCTCTTCCAGCGCTCGTCGCGCCAGGCCAATGGCCATGCCACATTGCAGCGCCAGGAATGCCGGGCGTACCCGTGGCAGGAATTCGCGTGCATTTTCATGCAGCAACCAGTCGCGGGC  >  NC_002947/4134127‑4134420

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