Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,298,632	A→C	5.7%	R490R (CGT→CGG)	recD ←	multifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,298,632	0	A	C	5.7%	267.1 / 7.2	88	R490R (CGT→CGG)	recD	multifunctional ATP‑dependent DNA helicase/dsDNA/ssDNA exonuclease/ssDNA endonuclease
Reads supporting (aligned to +/- strand):  ref base A (37/46);  new base C (3/2);  total (40/48)
Fisher's exact test for biased strand distribution p-value = 6.56e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.05e-01

GCGCCACCCGCTCGTTGAGGCCTTCGACACCCCAGGCGCCCCGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGTTCGTGCTTGAGCGAAAGGCTGTGCACGTCGGCAGGCGGCGCGGCCAGCAG  >  NC_002947/5298484‑5298778                                                                                                                                                     |                                                                                                                                                     gcgcCACCCGCTCGTTGAGGCCTTCGACACCCCAGGCGCCCCGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTaa                                                                                                                                                      >  3:253405/1‑149 (MQ=255)   gcCACCCGCTCGTTGAGGCCTTCGACACCCCAGGCGCCCCGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACg                                                                                                                                                    <  6:191540/149‑1 (MQ=255)     cACCCGCGCGTTGAGGCCTTCGACACCCCAGGCGCCCCGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGc                                                                                                                                                  <  6:213644/149‑1 (MQ=255)        ccGCTCGTTGAGGCCTTCGACACCCCAGGCGCCCCGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGcc                                                                                                                                                 >  3:193425/1‑147 (MQ=255)         cGCTCGTTGAGGCCTTCGACACCCCAGGCGCCCCGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAg                                                                                                                                              <  1:225924/149‑1 (MQ=255)           cTCGTTGAGGCCTTCGACACCCCAGGCGCCCCGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGTTACCGGccggcac                                                                                                                                            <  2:312467/149‑6 (MQ=255)               ttGAGGCCTTCTACAGCCCAGGCGGCCCGGCGCACCGCTAACAGCAAGTGGAAATCGTCAAAACTGCGCAACACCTTGCCCGCCGTTTGCTCCCAGGCGGGGTAGTCATACGCGGTATCATGTGCCGGCCGGTAACGGCCAAGGGTgcg                                                                                                                                        <  3:343206/149‑1 (MQ=255)                    gCCTTCGACACCCCAGGCGCCCCGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGt                                                                                                                                   <  6:370589/149‑1 (MQ=255)                      cTTCGACACCCCAGGCGCCCCGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCTCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAg                                                                                                                                 >  8:306848/1‑149 (MQ=255)                            cacCCCAGGCGCCCCGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGACCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCgg                                                                                                                           >  6:9255/1‑149 (MQ=255)                             acCCCAGGCGCCCCGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGt                                                                                                                          <  6:356640/149‑1 (MQ=255)                              ccccAGGCGCCCCGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTa                                                                                                                         <  5:6223/149‑1 (MQ=255)                               cccAGGCGCCCCGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCACCACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGGATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAg                                                                                                                        >  7:272970/1‑149 (MQ=255)                               cccAGGCGCCCCGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAg                                                                                                                        >  6:287341/1‑149 (MQ=255)                                ccAGGCGCCCCGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGc                                                                                                                       >  6:255053/1‑149 (MQ=255)                                 cAGGCGCCCCGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGc‑c                                                                                                                     <  2:324763/149‑1 (MQ=255)                                  aGGCGCCCCGGCGCACCGCGCACAGCAACTGGAGATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGCTACCTGCgcaccctt                                                                                                                     <  6:60838/149‑8 (MQ=255)                                   ggCGCCCCGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑Ctt                                                                                                                   <  4:215571/149‑1 (MQ=255)                                      gcCCCGTCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCg                                                                                                                <  4:10653/149‑1 (MQ=255)                                      gcCCCGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCg                                                                                                                <  8:53418/149‑1 (MQ=255)                                      gcCCCGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCg                                                                                                                >  8:237707/1‑149 (MQ=255)                                      gcCCCGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCg                                                                                                                <  5:9255/149‑1 (MQ=255)                                       ccccGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCgg                                                                                                               >  2:289671/1‑149 (MQ=255)                                        cccGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAAATGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGt                                                                                                              >  2:100768/1‑149 (MQ=255)                                         ccGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTc                                                                                                             <  3:370804/149‑1 (MQ=255)                                          cGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACg                                                                                                                                                    <  7:282587/110‑1 (MQ=255)                                          cGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACg                                                                                                                                                    >  8:282587/1‑110 (MQ=255)                                             cgcACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATc                                                                                                         <  4:193425/149‑1 (MQ=255)                                                 ccGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCAt                                                                                                          >  4:381498/1‑144 (MQ=255)                                                 ccGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCAt                                                                                                          <  3:381498/144‑1 (MQ=255)                                                 ccGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTa                                                                                                     >  3:356751/1‑149 (MQ=255)                                                 ccGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTa                                                                                                     >  1:239984/1‑149 (MQ=255)                                                 ccGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCCTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTa                                                                                                     >  5:48615/1‑149 (MQ=255)                                                  cgcgCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTAc                                                                                                    >  2:10227/1‑149 (MQ=255)                                                    cgcACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACcg                                                                                                  >  7:112605/1‑149 (MQ=255)                                                     gcACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGTGTCGTCAAAGGAGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATAGGTACcg‑c                                                                                                >  3:185148/1‑149 (MQ=255)                                                       acaGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑Cgg                                                                                              <  5:287341/149‑1 (MQ=255)                                                           cAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGtt                                                                                            >  4:338703/1‑147 (MQ=255)                                                           cAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCa                                                                                          <  5:255053/149‑1 (MQ=255)                                                            aaCTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAg                                                                                         <  3:24274/149‑1 (MQ=255)                                                            aaCTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAg                                                                                         <  5:326173/149‑1 (MQ=255)                                                             aCTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAgg                                                                                        <  4:356751/149‑1 (MQ=255)                                                             aCTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAgg                                                                                        <  1:10227/149‑1 (MQ=255)                                                                ggAAATCATCAAAACTGCGCAACACATTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCg                                                                                     >  3:162231/1‑149 (MQ=255)                                                                 gAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGcc                                                                                      <  2:386879/147‑1 (MQ=255)                                                                 gAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGt                                                                                    >  5:288874/1‑149 (MQ=255)                                                                    aTCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGTTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGa                                                                                 <  7:306848/149‑1 (MQ=255)                                                                    aTCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGa                                                                                 <  3:338703/149‑1 (MQ=255)                                                                    aTCCTCAAAACTGCGCAACACCCTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGa                                                                                 <  6:48615/149‑1 (MQ=255)                                                                       cTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGt                                                                               <  7:237707/148‑1 (MQ=255)                                                                        tCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTaa                                                                             >  8:309361/1‑149 (MQ=255)                                                                        tCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTaa                                                                             >  5:243685/1‑149 (MQ=255)                                                                        tCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTaa                                                                             >  2:319441/1‑149 (MQ=255)                                                                          aaaaCTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTaa                                                                                                                                                      <  4:82349/76‑1 (MQ=255)                                                                          aaaaCTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTaa                                                                             <  7:309361/147‑1 (MQ=255)                                                                          aaaaCTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACa                                                                           <  6:288874/149‑1 (MQ=255)                                                                           aaaCTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTg                                                                                                                  <  1:13823/110‑1 (MQ=255)                                                                           aaaCTGCGCAACACCTTGCCCGCCCATTGCTCACAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTg                                                                                                                  >  2:13823/1‑110 (MQ=255)                                                                            aaCTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAgg                                                                         >  5:189643/1‑149 (MQ=255)                                                                               tGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCgg                                                                      >  2:248727/1‑149 (MQ=255)                                                                                     acacGTTGCCCGCCCATTGCTGCGAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAAc                                                                <  4:43301/149‑1 (MQ=255)                                                                                                 cccATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGt                                                    >  3:43669/1‑149 (MQ=255)                                                                                                   cATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCTGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGTTc                                                  <  8:112605/149‑1 (MQ=255)                                                                                                   cATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCTGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGTTc                                                  <  5:87861/149‑1 (MQ=255)                                                                                                     ttGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGGCGAGTAACAGGCGGTCGAACGCACGATCCTGTTctt                                                <  4:162231/149‑3 (MQ=255)                                                                                                      tGCTCCCAGGCGGGGTCGTCACCGGCGGTAGCAAGTCCCGGCCGGTAACGGCCAAGGCGGCGCAGGTAGCTGCGGTAGC‑CTTGCGCTCCATCGCTACCG‑CGGTTCAGGCCGCCGAGTAACAGGCGGTCGAACGCACGATCCTGCTCGTg                                               <  6:134367/149‑1 (MQ=255)                                                                                                       gCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCGCCGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGTTCGt                                                <  1:289671/148‑1 (MQ=255)                                                                                                       gCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGTTCGTg                                               <  8:278543/148‑1 (MQ=255)                                                                                                        cTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGTTCGTGc                                              <  1:248727/148‑1 (MQ=255)                                                                                                        cTCACAGGCGGGGGCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTAGAGTAACAGGTGGTCGAACGCACGATCCTGTTCGTGc                                              <  4:185148/148‑1 (MQ=255)                                                                                                              ggCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAg                                                                                         >  3:28128/1‑99 (MQ=255)                                                                                                              ggCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAg                                                                                         <  4:28128/99‑1 (MQ=255)                                                                                                                 ggggTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGTTCGTGCTTGAGCGaaa                                    <  4:43669/149‑1 (MQ=255)                                                                                                                  gggTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGTTCGTGCTTGAGCGAAAg                                   >  8:50658/1‑149 (MQ=255)                                                                                                                       gtcAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGTTCGTGCTTGAGCGAAAGGCtg                               <  7:50658/148‑1 (MQ=255)                                                                                                                        tcAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGTTCGTGCTTGAGCGAAAGGCtgtg                             >  8:366739/1‑149 (MQ=255)                                                                                                                          aaaGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGTTCGTGCTTGAGCGAAAGGCTGTGCa                           >  8:322728/1‑149 (MQ=255)                                                                                                                           aaGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGTTCGTGCTTGAGCGAAAGGCTGTGCAc                          <  1:100768/149‑1 (MQ=255)                                                                                                                                    tCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGTTCGTGCTTGAGCGAAAGGCTGTGCACGTCGGCAgg                 >  4:98091/1‑149 (MQ=255)                                                                                                                                      aaGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGTTCGTGc                                              >  7:353263/1‑118 (MQ=255)                                                                                                                                      aaGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGTTCGTGc                                              <  8:353263/118‑1 (MQ=255)                                                                                                                                      aaGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGTTCGTGCTTGAGCGAAAGGCTGTGCACGTCGGCAggcg               >  8:322769/1‑149 (MQ=255)                                                                                                                                         ccttggccgtTACCGGCCAAGGGTGCGCAGGTAGCTGCGGTAGCGCTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGt                                                                     >  1:81688/5‑93 (MQ=255)                                                                                                                                         ccttggccgtTACCGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGt                                                                     <  2:81688/88‑1 (MQ=255)                                                                                                                                         ccttggccgtTACCGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGTTCGTGCTTGAGCGAAAGGCTGTGCACGTCGGCAGgcggcg            >  3:328180/5‑149 (MQ=255)                                                                                                                                           ttggccgtTACCGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGTTCGTGCTTGAGCGAAAGGCTGTGCACGTCGGCAGgcggcgcg          >  8:8096/3‑149 (MQ=255)                                                                                                                                             ggccggTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCGATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGTTCGTGCTTGAGCGAAAGGCTGTGCACGTCGGCAGgcggcgcggc        <  6:189643/149‑1 (MQ=255)                                                                                                                                              gccggTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGTTCGTGCTTGAGCGAAAGGCTGTGCACGTCGGCAGgcggcgcggc        >  2:123948/1‑148 (MQ=255)                                                                                                                                                  gTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGTTCGTGCTTGAGCGAAAGGCTGTGCACGTCGGCAGGCGGCGCGGCcag     >  2:355703/1‑147 (MQ=255)                                                                                                                                                   tAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGCACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGTTCGTGCTTGAGCGAAAGGCTGTGCACGTCGGCAGGCGGCGCGGCcagcag  >  3:31891/1‑149 (MQ=255)                                                                                                                                                     aCGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGTTCGTGCTTGAGCGAAAGGCTGTGCACGTCGGCAGgcggcgc           <  3:157909/138‑1 (MQ=255)                                                                                                                                                     aCGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGTTCGTGCTTGAGCGAAAGGCTGTGCACGTCGGCAGgcggcgc           >  4:157909/1‑138 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GCGCCACCCGCTCGTTGAGGCCTTCGACACCCCAGGCGCCCCGGCGCACCGCGCACAGCAACTGGAAATCCTCAAAACTGCGCAACACCTTGCCCGCCCATTGCTCCCAGGCGGGGTCGTCAAAGGCGGTATCAAGTGCCGGCCGGTAACGGCCAAGGGTGCGCAGGTAGCTGCGGTAGC‑CTTGCGGTCCATCGGTACCG‑CGGTTCAGGCCGTCGAGTAACAGGCGGTCGAACGCACGATCCTGTTCGTGCTTGAGCGAAAGGCTGTGCACGTCGGCAGGCGGCGCGGCCAGCAG  >  NC_002947/5298484‑5298778

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