Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,522,239	T→G	7.4%	F30V (TTC→GTC) ‡	idh →	isocitrate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,522,239	0	T	G	7.4%	206.6 / 7.4	68	F30V (TTC→GTC) ‡	idh	isocitrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base T (30/33);  new base G (3/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 6.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.09e-01

CGCCTACCCTGCGCACCACGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTACCACCCCTGAA  >  NC_002947/4522094‑4522385                                                                                                                                                  |                                                                                                                                                   cgcCTACCCTGCGCACCACGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCa                                                                                                                                                 <  1:15318/149‑1 (MQ=255)  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:357052/149‑1 (MQ=255)                   cGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCg                                                                                                                               <  3:313872/149‑1 (MQ=255)                       aCTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGt                                                                                                                           <  7:163130/149‑1 (MQ=255)                       aCTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGt                                                                                                                           >  7:362581/1‑149 (MQ=255)                        cTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTc                                                                                                                          >  6:35100/1‑149 (MQ=255)                        cTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTc                                                                                                                          <  2:238590/149‑1 (MQ=255)                         tCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTc                                                                                                                          >  6:145347/1‑148 (MQ=255)                             gCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGaaa                                                                                                                      >  4:319553/1‑148 (MQ=255)                                cGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTcg                                                                                                                  >  4:66158/1‑149 (MQ=255)                                cGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTcg                                                                                                                  >  5:177353/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                                                                                   >  6:109927/1‑149 (MQ=255)                                                               ccACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTgcc                                                                                   <  3:66158/149‑1 (MQ=255)                                                                 aCCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTgc                                                                                    >  2:200284/1‑146 (MQ=255)                                                                 aCCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTgc                                                                                    <  1:200284/146‑1 (MQ=255)                                                                  cccGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTgccgcc                                                                                <  5:362996/149‑1 (MQ=255)                                                                     gTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTc                                                                             <  5:120150/149‑1 (MQ=255)                                                                         cAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCg                                                                          <  1:86701/148‑1 (MQ=255)                                                                          aaGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTgc                                                                                    <  4:6304/137‑1 (MQ=255)                                                                          aaGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTgc                                                                                    >  3:6304/1‑137 (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                         <  4:168381/149‑1 (MQ=255)                                                                                                                         cTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTgg                         <  7:244570/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)                                                                                                                                                  |                                                                                                                                                   CGCCTACCCTGCGCACCACGAGACTCGAGCACGCTCAGCACACAGAGAGTTAACCCGCATGCCCACCCGTTCCAAGATCATCTATACCTTCACCGACGAAGCCCCCGCCCTCGCCACTTACTCGCTGCTGCCGATCATCGAAGCGTTCACCGCTTCGGCTGACATCGCTGTCGAAACTCGCGACATTTCCCTGGCTGGCCGTATCCTTGCCGCCTTCCCGGAGCAACTGGGCGCAGAGAAGCAAGTAGGCGATCACCTGGCGGAACTGGGCCAGCTGGCTACCACCCCTGAA  >  NC_002947/4522094‑4522385

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