Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,741,283	C→G	intergenic (+474/+278)	gltA → / ← yeiW	citrate synthase/putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,741,283	0	C	G	86.4%	179.4 / 17.4	66	intergenic (+474/+278)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (6/3);  new base G (27/30);  total (33/33)
Fisher's exact test for biased strand distribution p-value = 4.75e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

GTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CC‑TCC‑GC‑CAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCT  >  NC_002947/4741155‑4741406                                                                                                                                      |                                                                                                                              gTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑TTC‑TG‑CGACTTGCATTCGTGaatccag                                                                                                          >  1:2438691/1‑145 (MQ=255)            aTCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑TTC‑TG‑CGACTTGCATTCGTGAA‑TCCAGGCACAGTcccc                                                                                              >  2:351622/1‑150 (MQ=255)             tCCGCCTCCTTTATGGTCGTAATGGGCAGCACGGAATCACCGTGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑TTC‑TG‑CGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCg                                                                                             <  2:1069787/150‑1 (MQ=255)               cGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑TTCTGC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGa                                                                                           <  1:2449091/150‑1 (MQ=255)                gCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCT‑TCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGac                                                                                          >  1:2714321/1‑150 (MQ=255)                gCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACGGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑TTC‑TG‑CGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGac                                                                                          >  1:2662072/1‑150 (MQ=255)                gCCTCCTTGATGGTCGCAACGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGac                                                                                          >  2:2432392/1‑150 (MQ=255)                  ctcctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGacac                                                                                        >  1:2092425/1‑150 (MQ=255)                    cctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGt                                                                                      >  2:2528252/1‑150 (MQ=255)                      ttGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTTGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑TTC‑TG‑CGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGtcc                                                                                    >  2:1872002/1‑150 (MQ=255)                      ttGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGtcc                                                                                    <  1:1119175/150‑1 (MQ=255)                      ttGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑TTCTGC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGtcc                                                                                    >  2:441132/1‑150 (MQ=255)                              cGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAg                                                                            <  1:2293204/150‑1 (MQ=255)                                      tcagcaCGGAATCACGCGGCTTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑TTC‑TG‑CGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGa                                                                    <  1:1093165/149‑1 (MQ=255)                                      gcagcaCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGa                                                                    <  2:724115/150‑1 (MQ=255)                                        agcaCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGAcc                                                                  <  2:2096657/150‑1 (MQ=255)                                         gcaCGGAATCGCGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCt                                                                 >  2:2555691/1‑150 (MQ=255)                                          caCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCtt                                                                >  2:2615227/1‑150 (MQ=255)                                          caCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCtt                                                                >  1:1666752/1‑150 (MQ=255)                                           aCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg                                                               <  1:2615227/150‑1 (MQ=255)                                           aCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCCCAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg                                                               <  2:2644124/150‑1 (MQ=255)                                            cGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑TTCTGC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTgg                                                              >  1:2284367/1‑150 (MQ=255)                                             ggAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGc                                                             <  1:862602/150‑1 (MQ=255)                                              gAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGcc                                                            >  1:1953998/1‑150 (MQ=255)                                                  cACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcg                                                        <  2:2308743/150‑1 (MQ=255)                                                      cggcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGgcgc                                                    <  2:2759002/150‑1 (MQ=255)                                                        gcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCtt                                                                >  1:1696775/1‑136 (MQ=255)                                                        gcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCtt                                                                <  2:1696775/136‑1 (MQ=255)                                                          gTCACTGTCCGTGGCGTATTGAAATGGTTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTg                                                >  2:1507924/1‑150 (MQ=255)                                                           tCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGa                                               >  1:2361738/1‑150 (MQ=255)                                                                  ccATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GCTTCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGt                                        <  2:814480/150‑1 (MQ=255)                                                                       gCGTATTGAACTTGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                                   <  1:1872002/150‑1 (MQ=255)                                                                        cGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACAcc                                  <  2:2092425/150‑1 (MQ=255)                                                                          tATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc                                >  1:1869456/1‑150 (MQ=255)                                                                          tATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc                                <  2:2438691/150‑1 (MQ=255)                                                                           aTTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgcg                               >  1:695135/1‑150 (MQ=255)                                                                            ttGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGa                              <  2:1092841/150‑1 (MQ=255)                                                                            ttGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGc                                                             <  2:1108615/119‑1 (MQ=255)                                                                            ttGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CA‑GC‑TTCTGC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGc                                                             >  1:1108615/1‑119 (MQ=255)                                                                                  tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGcc                        >  2:2104616/1‑150 (MQ=255)                                                                                  tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGcc                        <  1:1693565/150‑1 (MQ=255)                                                                                  tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑GCTTCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGcc                        >  1:541695/1‑150 (MQ=255)                                                                                    ggTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGcc                                                   <  1:1435114/121‑1 (MQ=255)                                                                                    ggTATAGGTCCAGCTCTCGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGcc                                                   >  2:1435114/1‑121 (MQ=255)                                                                                      tataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg                                                               <  1:1262766/107‑1 (MQ=255)                                                                                      tataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg                                                               >  2:1262766/1‑107 (MQ=255)                                                                                       ataGGTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTa                   >  2:2435824/1‑150 (MQ=255)                                                                                          ggTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACtt                <  2:1541124/150‑1 (MQ=255)                                                                                          ggTCCAGCTCTGGTCACCTTCGTGCTTGAG‑CAGCT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACtt                <  1:2528252/150‑1 (MQ=255)                                                                                           gTCCAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTg               <  1:2555691/150‑1 (MQ=255)                                                                                            tCCAGCGCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTgg              >  2:2575053/1‑150 (MQ=255)                                                                                             ccAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGa             >  1:2648729/1‑150 (MQ=255)                                                                                             ccAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc                         >  1:2775361/1‑138 (MQ=255)                                                                                             ccAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc                         <  2:2775361/138‑1 (MQ=255)                                                                                              cAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACtt                >  2:2447173/1‑146 (MQ=255)                                                                                              cAGCTCTGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGGGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACtt                <  1:2447173/146‑1 (MQ=255)                                                                                                 ctctGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCg         >  1:420651/1‑150 (MQ=255)                                                                                                  tctGGTCACCTTCGTGCTTGAG‑CA‑GC‑TTC‑TG‑CGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCgg        <  1:2143581/150‑1 (MQ=255)                                                                                                   ctGGTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGc       <  2:1869456/150‑1 (MQ=255)                                                                                                      gTCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCAc    <  1:2185067/150‑1 (MQ=255)                                                                                                      gTCACCTTCGTGCTTGAG‑CA‑GC‑TTA‑TCGCGAC‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCAc    >  1:73442/1‑150 (MQ=255)                                                                                                       tCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc                         <  1:1729139/128‑1 (MQ=255)                                                                                                       tCACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc                                <  2:1320055/121‑1 (MQ=255)                                                                                                       tCACCTTCGTGCTTGAG‑CA‑GC‑TTCTGC‑GACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc                                >  1:1320055/1‑121 (MQ=255)                                                                                                       tCAACTTCGTGCTTGAG‑CA‑GC‑TTC‑TG‑CGACTTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc                         >  2:1729139/1‑128 (MQ=255)                                                                                                        cACCTTCGTGCTTGAGCAG‑CT‑TCT‑GC‑GACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCt  <  2:2193086/150‑1 (MQ=255)                                                                                                                                      |                                                                                                                              GTCAGGCTGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCTTGAG‑AA‑CC‑TCC‑GC‑CAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCT  >  NC_002947/4741155‑4741406

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