Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,521,977	G→T	7.3%	intergenic (‑240/‑175)	icd ← / → idh	NADP(+)‑specific isocitrate dehydrogenase/isocitrate dehydrogenase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,521,977	0	G	T	7.3%	212.6 / 6.2	68	intergenic (‑240/‑175)	icd/idh	NADP(+)‑specific isocitrate dehydrogenase/isocitrate dehydrogenase
Reads supporting (aligned to +/- strand):  ref base G (33/30);  new base T (3/2);  total (36/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.75e-01

GCATTGCTGCTTAGTCTTTGGTCTTATAAATGGGCCGATATCACTTGGCCTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTTCCGCTCATGGCGCTGCCAGGGCGATGCGCCTACCCTGCGCACCACGAGACTCGAG  >  NC_002947/4521835‑4522122                                                                                                                                               |                                                                                                                                                   gcATTGCTGCTTAGTCTTTGGTCTTATAAATGGGCCGATATCACTTGGCCTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGc                                                                                                                                              <  7:161046/149‑1 (MQ=255)    ttGCTGCTTAGTCTTTGGTCTTATAAATGGGCCGATATCACTTGGCCTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAAATCCAGGAAAATGCGGACCGCCCTAGGCCATTCATGGccta                                                                                                                                           <  4:275777/149‑3 (MQ=255)              tgcTTTGGTCTTAAAAATGGGCCGATATCACTTGGCCTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑ggg                                                                                                                                <  8:175460/147‑1 (MQ=255)               tCTTTGGTCTTATAAATGGGCCGATATCACTTGGCCTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGgt                                                                                                                               <  6:94134/149‑1 (MQ=255)                     gTCTTATAAATGGGCCGATATCACTTGGCCTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCACGACCGAAGGGTCATCGGGGCGAATCCCTGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTAc                                                                                                                         <  3:130269/149‑1 (MQ=255)                      tCTTATAAATGGGCCGATATCACTTGGCCTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTAcc                                                                                                                        >  6:191867/1‑149 (MQ=255)                          ataAATGGGCCGATATCACTTGGCCTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGc                                                                                                                      <  3:120508/147‑1 (MQ=255)                          ataAATGGGCCGATATCACTTGGCCTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAg                                                                                                                    >  5:326184/1‑149 (MQ=255)                           taAATGGGCCGATATCACTTGGCCTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAg                                                                                                                    <  8:21487/148‑1 (MQ=255)                           taAATGGGCCGATATCACTTGGCCTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGc                                                                                                                   >  6:22854/1‑149 (MQ=255)                           taAATGGGCCGATATCACTTGGCCTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGc                                                                                                                   >  6:364702/1‑149 (MQ=255)                            aaaTGGGCCGATATCACTTGGCCTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTa                                                                                                                          >  7:160588/1‑141 (MQ=255)                            aaaTGGGCCGATATCACTTGGCCTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTa                                                                                                                          <  8:160588/141‑1 (MQ=255)                             aaTGGGCCGATATCACTTGGCCTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCt                                                                                                                  >  2:136801/1‑148 (MQ=255)                             aaTGGGCCGATATCACTTGGCCTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCt                                                                                                                  >  8:84863/1‑148 (MQ=255)                                   ccGATATCACTTGGCCTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACt                                                                                                           >  5:29610/1‑149 (MQ=255)                                    cGATATCACTTGGCCTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACtt                                                                                                          <  7:49501/149‑1 (MQ=255)                                           ccTTGGCCTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTc                                                                                                   <  5:280902/148‑1 (MQ=255)                                            cTTGGCCTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATaaa                                                                                                                                        >  2:358188/1‑112 (MQ=255)                                            cTTGGCCTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATaaa                                                                                                                                        <  1:358188/112‑1 (MQ=255)                                                gCCTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATaaa                                                                                                                                        <  7:181503/108‑1 (MQ=255)                                                gCCTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATaaa                                                                                                                                        >  8:181503/1‑108 (MQ=255)                                                gCCTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACt                                                                                               <  6:59049/148‑1 (MQ=255)                                                 ccTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGa                                                                                             <  6:378414/149‑1 (MQ=255)                                                  cTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCGCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGAt                                                                                            <  5:4655/149‑1 (MQ=255)                                                   ttGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATg                                                                                           >  5:74364/1‑149 (MQ=255)                                                       tGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATc                                                                                       >  3:125748/1‑149 (MQ=255)                                                              cccAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCAt                                                                                 >  7:9887/1‑148 (MQ=255)                                                               ccAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCa                                                                               <  5:191867/149‑1 (MQ=255)                                                                cAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCAc                                                                              >  4:204467/1‑149 (MQ=255)                                                                  aaTGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCg                                                                            >  1:312272/1‑149 (MQ=255)                                                                      gTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTc                                                                         <  2:312272/148‑1 (MQ=255)                                                                         tGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCa                                                                      <  1:136801/148‑1 (MQ=255)                                                                          ggTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAg                                                                     <  5:345005/148‑1 (MQ=255)                                                                            tataCTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCCGCct                                                                  >  7:335747/1‑149 (MQ=255)                                                                              taCTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCctct                                                                <  7:84863/149‑1 (MQ=255)                                                                                cTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCtctc                                                               <  6:29610/148‑1 (MQ=255)                                                                                 tCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCg                                                              >  7:2021/1‑148 (MQ=255)                                                                                 tCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGa                                                             >  2:214912/1‑149 (MQ=255)                                                                                 tCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGa                                                             >  3:358375/1‑149 (MQ=255)                                                                                          aCCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGGCGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTTCCGCt                                                    >  8:240393/1‑149 (MQ=255)                                                                                            cGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTTCCGCTCa                                                  >  5:138955/1‑149 (MQ=255)                                                                                            cGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTTCCGCTCa                                                  <  4:125748/149‑1 (MQ=255)                                                                                              aaGGGTCATCGGGGCGAATCCCACGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACATTCCGCTCATg                                                <  6:326184/149‑1 (MQ=255)                                                                                                   tCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTTCCGCTCATGGCGCt                                           >  2:71823/1‑149 (MQ=255)                                                                                                          ggCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTTCCGCTCATGGCGCTGCCAggg                                    <  8:2021/149‑1 (MQ=255)                                                                                                          ggCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTTCCGCTCATGGCGCTGCCAggg                                    <  3:204467/149‑1 (MQ=255)                                                                                                          ggCGAATCCCAGGAAAATGCGGACCGCCCTACGCCATTAATGTCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCTCTCGCAGACTCTGGACTTTCCGCTCATGGCGCTGCCAggg                                    <  7:311554/149‑1 (MQ=255)                                                                                                              aaTCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTTCCGCTCATGGCGCTGCCAGGGCGAt                                <  5:22854/149‑1 (MQ=255)                                                                                                              aaTCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTTCCGCTCATGGCGCTGCCAGGGCGAt                                <  5:364702/149‑1 (MQ=255)                                                                                                               aTCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTTCCGCTCATGGCGCTGCCAGGGCGATg                               >  3:338386/1‑149 (MQ=255)                                                                                                               aTCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCCCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTACCGCTCATGGCGCTGCCAGGGCGATg                               >  1:386550/1‑149 (MQ=255)                                                                                                                  ccAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTTCCGCTCATGGCGCTGCCAGGGCGATgcgc                            >  2:263614/1‑149 (MQ=255)                                                                                                                      gAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTTCCGCTCATGGCGCTGCCAGGGCGATGTGCCTAc                        <  4:358375/149‑1 (MQ=255)                                                                                                                      gAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTTCCGCTCATGGCGCTGCCAGGGCGATGCGCCTAc                        >  8:267619/1‑149 (MQ=255)                                                                                                                      gAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GCGTTGTTAACGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTTCCGCTCATGGCGCTGCCAGGGCGATGCGCCTAc                        >  5:101918/1‑149 (MQ=255)                                                                                                                              ggACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTTCCGCTCATGGCGCTGCCAGGGCGATGCGCCTACCCTGCGca                >  5:286386/1‑149 (MQ=255)                                                                                                                                aCCGCCCTAGGCGATGAATGGCCATAATGCCT‑GGGTTGTTAGCGCAGCTCAGCTCTTGTCGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTTCCGCTCATGGCGCTGCCAGGGCGATGCGCCTACCCTGCGcacc              <  2:386550/149‑1 (MQ=255)                                                                                                                                aCCGCCCTAGGCCATGAATGGCCATAAAGCCTGGGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGCCTTTCCGCTCATGGCGCTGCCAGGGCGATGCGCCTACCCTGCGcac               >  1:46973/1‑149 (MQ=255)                                                                                                                                    tttatggCCATTCATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTTCCGCTCATGGCGCTGCCAGGGCGATGCGCCTACCCTGCGcaccac            >  2:285591/6‑147 (MQ=255)                                                                                                                                    cccTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTTCCGCTCATGGCGCTGCCAGGGCGATGCGCCTACCCTGCGCACCACga          <  8:335747/149‑1 (MQ=255)                                                                                                                                     ttatggCCATTCATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTTCCGCTCATGGCGCTGCCAGGGCGATGCGCCTACCCTGCGcaccac            >  3:6290/5‑146 (MQ=255)                                                                                                                                       atggCCATTCATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTTCCGCTCATGGCGCTGCCAGGGCGATGCGCCTACCCTGCGCACCACGAGAc       >  2:37662/3‑149 (MQ=255)                                                                                                                                          gCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTTCCGCTCATGGCGCTGCCAGGGCGATGCGCCTACCCTGCGCACCACGAGACTc     >  4:231602/1‑148 (MQ=255)                                                                                                                                           ccctGAATCGCCATAAAGCCT‑TGGTGCTTACCGCAGCTCAGCACGTGACGCTCGACTGATGCATCCACCATCTCCCCTCGCAGCCTCGCGACTTCCCGCTCATGGCGCTGCCAGTGCCATGCGCCTACCCTGCGCACCACGAGCCTCGa   <  6:101918/146‑1 (MQ=255)                                                                                                                                            cATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTTCCGCTCATGGCGCTGCCAGGGCGATGCGCCTACCCTGCGCACCACGAGACTCGAg  >  3:310967/1‑149 (MQ=255)                                                                                                                                            cATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTTCCGCTAATGGCGCTGCCAGGGCGATGCGCCTACCCTGCGCACCACGAGACTCGAg  >  3:125985/1‑149 (MQ=255)                                                                                                                                             aTGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTTCCGCTCATGGCGCTGCCAGGGCGATGCGCCTACCCTGCGCACCACGAGACTCGAg  >  4:89026/1‑148 (MQ=255)                                                                                                                                               |                                                                                                                                                   GCATTGCTGCTTAGTCTTTGGTCTTATAAATGGGCCGATATCACTTGGCCTTGCTGATTAGCCCAAATGGTTTGGTATACTCCGCCGCGACCGAAGGGTCATCGGGGCGAATCCCAGGAAAATGCGGACCGCCCTAGGCCATGAATGGCCATAAAGCCT‑GGGTTGTTACCGCAGCTCAGCACTTGACGCTCGACTGATGCATCCACCATCACCGCTCGCAGCCTCTCGACTTTCCGCTCATGGCGCTGCCAGGGCGATGCGCCTACCCTGCGCACCACGAGACTCGAG  >  NC_002947/4521835‑4522122

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑