Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	626,556	T→G	6.1%	K213Q (AAG→CAG)	dedA ←	inner membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	626,556	0	T	G	6.1%	265.1 / 6.5	86	K213Q (AAG→CAG)	dedA	inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (39/40);  new base G (2/3);  total (42/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

GTCGCTGATGCCGTCCAGCATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCC  >  NC_002947/626411‑626700                                                                                                                                                   |                                                                                                                                                  gTCGCTGATGCCGTCCAGCATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGc                                                                                                                                                <  2:55456/149‑1 (MQ=255) gTCGCTGATGCCGTCCAGCATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGc                                                                                                                                                <  2:101183/149‑1 (MQ=255)     cTGATGCCGTCCAGCATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTa                                                                                                                                            <  6:194676/149‑1 (MQ=255)     cTGATGCCGTCCAGCATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTGGGCCGTa                                                                                                                                            >  5:50127/1‑149 (MQ=255)      tGATGCCGTCCAGCATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTAc                                                                                                                                           >  2:2759/1‑149 (MQ=255)      tGATGCCGTCCAGCATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGGTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTAc                                                                                                                                           >  6:40930/1‑149 (MQ=255)           ccGTCCAGCATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGccc                                                                                                                                      <  2:320595/149‑1 (MQ=255)           ccGTCCAGCATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGcc                                                                                                                                       <  1:223100/148‑1 (MQ=255)            cGTCCAGCATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCa                                                                                                                                     >  8:103541/1‑149 (MQ=255)            cGTCCAGCATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCa                                                                                                                                     <  5:374107/149‑1 (MQ=255)             gTCCAGCATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCAA‑g                                                                                                                                   <  3:164447/149‑1 (MQ=255)               acAGCAATGGCAGGCCGTACGCGAGCGG‑CTGCCAATGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑Ggc                                                                                                                                 >  4:371750/2‑149 (MQ=37)                cAGCATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑Ggcg                                                                                                                                >  4:354950/1‑149 (MQ=255)                cAGCATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑Ggcg                                                                                                                                >  3:183661/1‑149 (MQ=255)                 aGCATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑Ggcgc                                                                                                                               <  6:389735/149‑1 (MQ=255)                 aGCATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTGGGCCGTACGGCCAA‑Ggctc                                                                                                                               <  6:50127/149‑3 (MQ=255)                  gCATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTGGGCCGTACGGCCAA‑Ggctca                                                                                                                              <  6:194666/149‑4 (MQ=255)                   cATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTGGGCCGTACGGCCAA‑Ggctcac                                                                                                                             <  8:313346/149‑5 (MQ=255)                        gCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACgg                                                                                                                                         >  2:161254/1‑133 (MQ=255)                        gCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACgg                                                                                                                                         <  1:161254/133‑1 (MQ=255)                        gCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACg                                                                                                                          >  8:55242/1‑147 (MQ=255)                          a‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACgcag                                                                                                                       >  1:48441/1‑148 (MQ=255)                            gCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACgcagc                                                                                                                      >  3:26915/1‑148 (MQ=255)                               gCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGc                                                                                                                  >  2:163025/1‑149 (MQ=255)                                   gcgTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCga                                                                                                               <  1:2759/148‑1 (MQ=255)                                    cgTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCgaga                                                                                                             <  8:142268/149‑1 (MQ=255)                                           ccTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCg                                                                                                      <  5:235210/149‑1 (MQ=255)                                             tGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCgg                                                                                                     >  2:43103/1‑148 (MQ=255)                                             tGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCg                                                                                                      >  5:184298/1‑147 (MQ=255)                                             tGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCg                                                                                                      <  6:184298/147‑1 (MQ=255)                                               ccATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGc                                                                                                    <  1:163025/147‑1 (MQ=255)                                               ccATTGCTGTGGCGTGATCGGGTAACGCGACAGACTACGCTTGCGCCGCCAGCCGCTGAATGACCACATGCAGATCCTTCTTTACCTGTACAGTGAGCCTTAGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCAc                                                                                                  >  2:298895/1‑149 (MQ=255)                                                 aTTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTGTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCAc                                                                                                  >  1:123627/1‑147 (MQ=255)                                                        tgGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAg                                                                                                                   <  1:77925/123‑1 (MQ=255)                                                        tgGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAg                                                                                                                   >  2:77925/1‑123 (MQ=255)                                                        tgGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGAc                                                                                           <  6:172878/147‑1 (MQ=255)                                                          gCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAg                                                                                         <  4:153157/147‑1 (MQ=255)                                                          gCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAg                                                                                         >  3:153157/1‑147 (MQ=255)                                                           cGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGgat                                                                                      <  7:83419/149‑1 (MQ=255)                                                             tGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGgatga                                                                                    >  5:256165/1‑149 (MQ=255)                                                              gATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGgatgat                                                                                   >  3:177383/1‑149 (MQ=255)                                                                  gggTAACGCGCCATAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACTTGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGatgatgatg                                                                               <  2:112597/149‑1 (MQ=255)                                                                  gggTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACTTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGatgatgatg                                                                               >  1:294434/1‑149 (MQ=255)                                                                    gTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGcc                                                                             <  4:26915/149‑1 (MQ=255)                                                                      aaCGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCcca                                                                           >  7:69695/1‑149 (MQ=255)                                                                       aCGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCccac                                                                          <  3:371750/149‑1 (MQ=255)                                                                         gcgcCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCccacca                                                                        <  2:123627/149‑1 (MQ=255)                                                                           gcCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑Ggc                                                                                                                                 <  6:310814/90‑1 (MQ=255)                                                                           gcCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑Ggc                                                                                                                                 >  5:310814/1‑90 (MQ=255)                                                                             cAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCAt                                                                    >  7:211331/1‑149 (MQ=255)                                                                                agTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGACCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAg                                                                 >  5:347307/1‑149 (MQ=255)                                                                                 gTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGc                                                                <  7:257394/149‑1 (MQ=255)                                                                                     gCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCa                                                                                                <  7:70133/113‑1 (MQ=255)                                                                                     gCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCa                                                                                                >  8:70133/1‑113 (MQ=255)                                                                                      cTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAg                                                           <  2:48441/149‑1 (MQ=255)                                                                                         gcgcCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGGGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGt                                                         <  3:296346/148‑1 (MQ=255)                                                                                                cAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCtt                                                  >  5:82083/1‑148 (MQ=255)                                                                                                 aGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTtga                                                <  5:204009/149‑1 (MQ=255)                                                                                                   gcgcTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTtgat                                               >  5:70809/1‑148 (MQ=255)                                                                                                    cgcTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTtgatg                                              >  8:24219/1‑148 (MQ=255)                                                                                                        gAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCACACGAAGCAAGCCGACAATCATCGGCAAAAGCGACATGATGATCATGACCCCCACCATCAGCGACAGGTGATGCATCATGAAAgg                                         >  5:8311/1‑149 (MQ=255)                                                                                                         aaTGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCATGGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATgg                                         >  1:106321/1‑149 (MQ=255)                                                                                                               cacaTGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTCGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTg                                  >  5:113249/1‑149 (MQ=255)                                                                                                                acaTGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATAGCACGTTGc                                 >  8:237392/1‑149 (MQ=255)                                                                                                                acaTGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTCGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGc                                 <  6:113249/149‑1 (MQ=255)                                                                                                                  aTGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCg                               <  7:94654/149‑1 (MQ=255)                                                                                                                   tGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACGAGCGACAGGATGATGATGCCGACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGGCga                              <  6:256165/149‑1 (MQ=255)                                                                                                                      cgttcgttCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCgaag                            >  2:282737/1‑148 (MQ=255)                                                                                                                       gttcgttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCgaagaa                          <  7:103541/149‑1 (MQ=255)                                                                                                                        ttcgttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGaagaag                         <  1:282737/149‑1 (MQ=255)                                                                                                                          cgttcTTTACCTGTTCAGTGAGCCTTGGCCGTACTGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCCTCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGt                        <  8:200516/148‑1 (MQ=255)                                                                                                                            ttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGc                     <  4:183661/149‑1 (MQ=255)                                                                                                                            ttcTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGc                     <  1:43103/149‑1 (MQ=255)                                                                                                                                 tACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGccc                   <  4:177383/146‑1 (MQ=255)                                                                                                                                  aCCTGTTCAGTGATGCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAg                <  5:40930/148‑1 (MQ=255)                                                                                                                                    cTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGaagaag                         >  4:73206/1‑137 (MQ=255)                                                                                                                                    cTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGaagaag                         <  3:73206/137‑1 (MQ=255)                                                                                                                                    cTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGt              >  2:361153/1‑148 (MQ=255)                                                                                                                                       ttCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGAc           >  8:218540/1‑148 (MQ=255)                                                                                                                                         cAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCCGCGACAGGTGCTGCTTGATGAATAGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCa         <  7:237392/148‑1 (MQ=255)                                                                                                                                            tGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCCGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGAAGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGgca     >  2:190156/1‑148 (MQ=255)                                                                                                                                            tGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGgcc     >  8:339449/1‑149 (MQ=255)                                                                                                                                             gcgCCTGGGCCGTACGGCCAA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGgccg    >  8:144937/3‑149 (MQ=255)                                                                                                                                                ccTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGcc                                >  7:70098/1‑118 (MQ=255)                                                                                                                                                ccTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGcc                                <  8:70098/118‑1 (MQ=255)                                                                                                                                                ccTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGgccgcc  >  5:391927/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                  GTCGCTGATGCCGTCCAGCATTGGCA‑GCCGCTCGCGTACGGCCTGCCATTGCTGTGGCGTGATCGGGTAACGCGCCAGAGTACGCTTGCGCCGCCAGGCGCTGAATGACCACATGCCGTTCGTTCTTTACCTGTTCAGTGAGCCTTGGCCGTACGGCCCA‑GGCGCCCACGCAGCAAGCCGAGAATCATCGGCACCAGCGACAGGATGATGATGCCCACCACCATCAGCGACAGGTGCTGCTTGATGAATGGCACGTTGCCGAAGAAGTAGCCCAAGGTGACCAGGCCGCC  >  NC_002947/626411‑626700

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