Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,522,243	C→A	7.2%	T31N (ACC→AAC)	idh →	isocitrate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,522,243	0	C	A	7.2%	213.8 / 3.0	67	T31N (ACC→AAC)	idh	isocitrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base C (31/31);  new base A (3/2);  total (34/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCCTACCCTGCGCACCACGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTACCACCCCTGAAGCCAA  >  NC_002947/4522095‑4522390                                                                                                                                                     |                                                                                                                                                    gcCTACCCTGCGCACCACGATACTCTAGCAAGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCAc                                                                                                                                                     <  6:428519/149‑1 (MQ=255)      cccTGCGCACCACGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCtt                                                                                                                                                <  6:69384/149‑1 (MQ=255)        cTGCGCACCACGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGGTGAACGCTTCg                                                                                                                                              <  3:319297/149‑1 (MQ=255)           cgcACCACGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCGCGCCCTCGCCACTAACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCt                                                                                                                                           <  5:360724/149‑1 (MQ=255)           cgcACCACGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGGTGAACGCTTCga                                                                                                                                             <  7:24106/147‑2 (MQ=255)             caccacGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGa                                                                                                                                         <  7:42575/149‑1 (MQ=255)              accacGAGACTCGAGCACGCTCATCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTc                                                                                                                                               <  4:335750/142‑1 (MQ=255)              accacGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTc                                                                                                                                               >  3:335750/1‑142 (MQ=255)                 acGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACAt                                                                                                                                      >  5:422992/1‑148 (MQ=255)                  cGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCg                                                                                                                                    <  3:313872/149‑1 (MQ=255)                  cGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCg                                                                                                                                    <  3:357052/149‑1 (MQ=255)                      aCTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGt                                                                                                                                <  7:163130/149‑1 (MQ=255)                      aCTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGt                                                                                                                                >  7:362581/1‑149 (MQ=255)                       cTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTc                                                                                                                               <  2:238590/149‑1 (MQ=255)                       cTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTc                                                                                                                               >  6:35100/1‑149 (MQ=255)                        tCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTc                                                                                                                               >  6:145347/1‑148 (MQ=255)                            gCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGaaa                                                                                                                           >  4:319553/1‑148 (MQ=255)                               cGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTcg                                                                                                                       >  5:177353/1‑149 (MQ=255)                               cGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTcg                                                                                                                       >  4:66158/1‑149 (MQ=255)                                gCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATGGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTcgc                                                                                                                      <  1:304040/149‑1 (MQ=255)                                gCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTcgc                                                                                                                      <  4:211381/149‑1 (MQ=255)                                    aGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACa                                                                                                                  <  4:392365/149‑1 (MQ=255)                                       acacaGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACAttt                                                                                                               <  3:41886/149‑1 (MQ=255)                                         acaGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTAACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTcc                                                                                                             >  4:80200/1‑149 (MQ=255)                                           agagagTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTACCGACATTTCCct                                                                                                           >  4:10828/1‑149 (MQ=255)                                           agagagTTAACCCGCAGGCCCACCCGGTCCAAGATCATCTATACCTCCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCct                                                                                                           <  2:254267/149‑1 (MQ=255)                                                gTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCctggctg                                                                                                      <  6:177353/149‑1 (MQ=255)                                                     cccGCATGCCCACCCGCTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCg                                                                                                  >  3:168381/1‑148 (MQ=255)                                                        gCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATc                                                                                              <  3:319553/149‑1 (MQ=255)                                                             cccACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCGATCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTgc                                                                                         <  1:363651/149‑1 (MQ=255)                                                             cccACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTgc                                                                                         <  8:317908/149‑1 (MQ=255)                                                              ccACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTgcc                                                                                        <  3:66158/149‑1 (MQ=255)                                                              ccACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTgcc                                                                                        >  6:109927/1‑149 (MQ=255)                                                                aCCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTgc                                                                                         <  1:200284/146‑1 (MQ=255)                                                                aCCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTgc                                                                                         >  2:200284/1‑146 (MQ=255)                                                                 cccGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTgccgcc                                                                                     <  5:362996/149‑1 (MQ=255)                                                                    gTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTc                                                                                  <  5:120150/149‑1 (MQ=255)                                                                        cAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCg                                                                               <  1:86701/148‑1 (MQ=255)                                                                         aaGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTgc                                                                                         >  3:6304/1‑137 (MQ=255)                                                                         aaGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTgc                                                                                         <  4:6304/137‑1 (MQ=255)                                                                             tcatcTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCaa                                                                         >  7:221625/1‑149 (MQ=255)                                                                             tcatcTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCaa                                                                         >  2:5996/1‑149 (MQ=255)                                                                                tcTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTg                                                                      >  3:202029/1‑149 (MQ=255)                                                                                   ataCCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGgc                                                                   >  8:124943/1‑149 (MQ=255)                                                                                    taCCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTg                                                                      <  5:145347/145‑1 (MQ=255)                                                                                       cTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGgcg                                                                  >  8:244570/1‑146 (MQ=255)                                                                                       cTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGgcg                                                                  >  7:26167/1‑146 (MQ=255)                                                                                          cACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCagag                                                             >  4:402208/1‑148 (MQ=255)                                                                                            ccGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAg                                                          >  6:74191/1‑149 (MQ=255)                                                                                             cgacgaAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAg                                                          >  6:7471/1‑148 (MQ=255)                                                                                                   aGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAg                                                   >  4:68775/1‑149 (MQ=255)                                                                                                     cccccGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGc                                                 >  2:12619/1‑149 (MQ=255)                                                                                                                gCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTgg                                       <  3:402208/148‑1 (MQ=255)                                                                                                                    cTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGa                                   >  3:81932/1‑148 (MQ=255)                                                                                                                       aCTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTg                               >  7:182584/1‑149 (MQ=255)                                                                                                                        cTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTgg                              <  7:244570/149‑1 (MQ=255)                                                                                                                        cTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTgg                              <  4:168381/149‑1 (MQ=255)                                                                                                                           gctgctGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGc                            <  8:26167/148‑1 (MQ=255)                                                                                                                              gctgcCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGc                        >  4:163819/1‑149 (MQ=255)                                                                                                                                    gATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTa                  <  8:221625/149‑1 (MQ=255)                                                                                                                                    gATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTACCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGAAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTa                  >  1:336021/1‑149 (MQ=255)                                                                                                                                     atcatcGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTac                 <  5:109927/149‑1 (MQ=255)                                                                                                                                       agccGAAGCGGTGAACGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCTGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTacca               >  1:100533/4‑149 (MQ=255)                                                                                                                                        gccGAAGCGGTGAACGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTaccac              >  6:376633/3‑149 (MQ=255)                                                                                                                                          cGAAGCGGTGAACGCTTCGGCTGACATCGCTGTCGAAACTCGCGACAGTTCCCTGGCTGGCCGTATCCTTGCCGGCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGGGGAACTGGGCCAGCTGGCTACCAccc            >  4:5778/1‑149 (MQ=255)                                                                                                                                            aaGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTACCACCCCt          >  4:277672/1‑149 (MQ=255)                                                                                                                                               cGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTACCACCCCTGaa       <  5:35100/149‑1 (MQ=255)                                                                                                                                                    aCCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTACCACCCCTGAAGCCaa  >  4:24908/1‑149 (MQ=255)                                                                                                                                                     ccGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGACGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTACCACCCCTGAAGCCaa  >  6:49877/1‑148 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GCCTACCCTGCGCACCACGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTACCACCCCTGAAGCCAA  >  NC_002947/4522095‑4522390

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