Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	4,741,277	C→T	intergenic (+468/+284)	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,277	0	C	T	78.6%	171.2 / 33.7	85	intergenic (+468/+284)	gltA/yeiW	citrate synthase/putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (10/8);  new base T (39/27);  total (49/36)
Fisher's exact test for biased strand distribution p-value = 7.94e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑AA‑CC‑TCC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACT  >  NC_002947/4741162‑4741410                                                                                                                        |                                                                                                                                          tGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAgg                                                                                                      >  1:1301951/1‑150 (MQ=255)      tCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑CA‑GC‑TTCTGCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCg                                                                                                 >  2:1016724/1‑150 (MQ=255)       ccGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑CA‑GC‑TTC‑TGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCgg                                                                                                >  1:106653/1‑150 (MQ=255)        cGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGa                                                                                               >  2:1005339/1‑150 (MQ=255)        cGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑CA‑GC‑TTCTGCGACT‑TGCATTCGTGAA‑TCCAGCCACACGCCCCGGa                                                                                               >  1:1598542/1‑150 (MQ=255)         gCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGac                                                                                              >  1:363665/1‑150 (MQ=255)           ctcctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGacac                                                                                            <  1:3511815/150‑1 (MQ=255)              ctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtc                                                                                         <  2:2330049/150‑1 (MQ=255)              ctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACATTCCCCAGACACGtc                                                                                         >  1:1595217/1‑150 (MQ=255)              ctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGtc                                                                                         >  1:3279145/1‑150 (MQ=255)              ctTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑CA‑GCTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGtc                                                                                         >  1:2104103/1‑150 (MQ=255)                 gATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtcctc                                                                                      <  2:1387847/150‑1 (MQ=255)                 gATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGtcctc                                                                                      >  1:1119265/1‑150 (MQ=255)                  aTGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtcctcc                                                                                     <  1:1005339/150‑1 (MQ=255)                  aTGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtcctcc                                                                                     <  2:3178017/150‑1 (MQ=255)                   tGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑CA‑GCTTCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTcctcct                                                                                    <  1:1219842/150‑1 (MQ=255)                   tGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑CA‑GC‑TTC‑TGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTcctcct                                                                                    <  2:635949/150‑1 (MQ=255)                         tAATGTGCAGCACGGAATCACGCAGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAgcg                                                                              >  1:1846643/1‑150 (MQ=255)                          aaTGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAgcgc                                                                             >  1:1618197/1‑150 (MQ=255)                             gtgCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTc                                                                          <  1:1293543/150‑1 (MQ=255)                                 agcaCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGAcc                                                                      >  1:1757107/1‑150 (MQ=255)                                  gcaCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCt                                                                     >  2:2326959/1‑150 (MQ=255)                                   caCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑CAGCT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCCGCGCCTCGACCtt                                                                    >  1:267805/1‑150 (MQ=255)                                   caCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑CAGCT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCtt                                                                    <  2:1301951/150‑1 (MQ=255)                                    aCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg                                                                   <  2:2413069/150‑1 (MQ=255)                                    aCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg                                                                   >  2:1048076/1‑150 (MQ=255)                                    aCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg                                                                   <  2:2104103/150‑1 (MQ=255)                                    aCGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑CA‑GC‑TTCTGCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg                                                                   <  2:2292581/150‑1 (MQ=255)                                     cGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑CAGCT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTgg                                                                  >  1:171953/1‑150 (MQ=255)                                      ggAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑CAGCT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGc                                                                 >  1:1409817/1‑150 (MQ=255)                                          tCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑CAGCT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgc                                                             >  2:1399895/1‑150 (MQ=255)                                           cACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtcct                                                                                       >  1:751053/1‑123 (MQ=255)                                           cACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑CA‑GC‑TTCTGCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGtcct                                                                                       <  2:751053/123‑1 (MQ=255)                                           cACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑CA‑GC‑TTCTGCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcg                                                            <  1:1016724/150‑1 (MQ=255)                                            aCGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑CA‑GC‑TTC‑TGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcgg                                                           >  1:3530272/1‑150 (MQ=255)                                              gcggcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAgcggcg                                                         >  1:3150688/1‑150 (MQ=255)                                               cggcgTCACCGGTCATGGCGTATTGAACTGGGTATAGGTCCAGATCTGGTC‑CCGTCGGTCTGTGATT‑TT‑GG‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGgcgc                                                        <  1:2185839/150‑1 (MQ=255)                                               cggcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑CAGCT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGgcgc                                                        <  1:3177325/150‑1 (MQ=255)                                                 gcgTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCt                                                      <  2:1757107/150‑1 (MQ=255)                                                    tCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑CAGCT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCAGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGa                                                   >  1:1266592/1‑150 (MQ=255)                                                     cACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGAt                                                  <  1:3305146/150‑1 (MQ=255)                                                      aCCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGTCACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATc                                                 >  2:3589824/1‑150 (MQ=255)                                                        cGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTTGCCATCGGCCCCTTGATCAc                                               >  2:3491220/1‑150 (MQ=255)                                                        cGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCAc                                               >  2:3148582/1‑150 (MQ=255)                                                         ggCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                                              <  2:1119265/150‑1 (MQ=255)                                                         ggCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                                              <  2:363665/150‑1 (MQ=255)                                                         ggCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑CA‑GC‑TTC‑TGCGACTTGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACg                                              <  2:1060637/150‑1 (MQ=255)                                                          gCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACgg                                             <  2:478747/150‑1 (MQ=255)                                                           ccATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtcc                                                                                        >  2:2732321/1‑106 (MQ=255)                                                           ccATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGtcc                                                                                        <  1:2732321/106‑1 (MQ=255)                                                            cATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTc                                           >  1:1851826/1‑150 (MQ=255)                                                             aTGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                                          >  1:3108045/1‑150 (MQ=255)                                                             aTGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTca                                          >  2:1055578/1‑150 (MQ=255)                                                                gCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCacac                                       <  1:915968/150‑1 (MQ=255)                                                                 cGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg                                                                   <  2:2485962/121‑1 (MQ=255)                                                                 cGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑CA‑GC‑TTCTGCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTg                                                                   >  1:2485962/1‑121 (MQ=255)                                                                   tATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGcc                                                       >  1:1502312/1‑131 (MQ=255)                                                                   tATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGcc                                                       <  2:1502312/131‑1 (MQ=255)                                                                   tATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACcgc                                    >  2:883294/1‑150 (MQ=255)                                                                      tGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGAt                                 >  1:3211042/1‑150 (MQ=255)                                                                       gAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATa                                >  2:1661056/1‑150 (MQ=255)                                                                        aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg                               >  2:1606038/1‑150 (MQ=255)                                                                        aaCTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg                               >  2:2034046/1‑150 (MQ=255)                                                                          cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc                             >  1:577811/1‑150 (MQ=255)                                                                          cTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑CA‑GCTTCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGc                             <  1:3148582/150‑1 (MQ=255)                                                                           tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAg                               <  2:2843699/147‑1 (MQ=255)                                                                           tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACAGCGATAg                               >  1:2843699/1‑147 (MQ=255)                                                                           tGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑CA‑GC‑TTC‑TGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGcc                            >  1:943997/1‑150 (MQ=255)                                                                             ggTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCtt                          <  1:883294/150‑1 (MQ=255)                                                                              gTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTg                         >  1:129150/1‑150 (MQ=255)                                                                              gTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTg                         <  1:1661056/150‑1 (MQ=255)                                                                                ataGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTa                       <  1:696726/150‑1 (MQ=255)                                                                                 taGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑CAGCT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTAc                      <  1:1399895/150‑1 (MQ=255)                                                                                   ggTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACtt                    >  1:2709985/1‑150 (MQ=255)                                                                                    gTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑CAGCC‑TCT‑GAGTGT‑TGACCTACTGAA‑TCCAGGCACAGGCCCCTGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGGCACACCGCGATAGGCCTTGTACTTg                   >  1:1487704/1‑150 (MQ=255)                                                                                     tCCAGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGAATTCGTGAA‑TCAAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTgg                  >  1:1019844/1‑150 (MQ=255)                                                                                     tCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑CA‑GC‑TTC‑TGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTgg                  >  1:2318023/1‑150 (MQ=255)                                                                                        aGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCTGGACACGTCCTCCTGCAGCGCCTc                                                                          >  2:694671/1‑91 (MQ=255)                                                                                        aGCTCTGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCTGGACACGTCCTCCTGCAGCGCCTc                                                                          <  1:694671/91‑1 (MQ=255)                                                                                         gCTCTGGTCACCTTCGTGCT‑TGA‑G‑CA‑GC‑TTC‑TAGGACTTACATTCGTGAA‑TCCAGGCACCGGCCCCGAACGCGTCCTCCTACACCGCTGCTACCTTGGCCAGCGGCACCTTGATCACGGTCACACCGCGATAGGACGTGTACGTGGAccc              >  1:2009110/1‑150 (MQ=255)                                                                                            ctGGTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCTGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACGTGGCCCCGGc           <  2:943997/150‑1 (MQ=255)                                                                                            ctGGTCACCTTCGTGCT‑TGA‑G‑CA‑GC‑TTC‑TGCGACTTGCATTCGTGAA‑TCCAGGCACAGGCCCCGTACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGGCCCGGc           <  2:106653/150‑1 (MQ=255)                                                                                              ggTCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCa         >  1:2174325/1‑150 (MQ=255)                                                                                                tCACCTTCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAATCCAGGC‑ACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCAcc       <  2:1987745/150‑1 (MQ=255)                                                                                                     ttCGTGCT‑TGA‑GCAG‑CT‑TCT‑GCGACT‑TGCATTCGTGAA‑TCCAGGCACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACCTGGACCCCGCCACCTCAcc  >  2:1212577/1‑149 (MQ=255)                                                                                                                        |                                                                                                                                          TGCCATCCGCCTCCTTGATGGTCGTAATGTGCAGCACGGAATCACGCGGCGTCACCGGCCATGGCGTATTGAACTGGGTATAGGTCCAGCTCTGGTCACCTTCGTGCT‑TGA‑G‑AA‑CC‑TCC‑GCCAAC‑TGCATTCGTGAA‑CCAAGCAACAGGCCCCGGACACGTCCTCCTGCAGCGCCTCGACCTTGGCCAGCGGCGCCTTGATCACGGTCACACCGCGATAGGCCTTGTACTTGGACCCGGCCACCTCACT  >  NC_002947/4741162‑4741410

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