Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	581,758	G→A	5.5%	A350T (GCC→ACC)	selA →	selenocysteine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	581,758	0	G	A	5.5%	222.3 / 2.1	73	A350T (GCC→ACC)	selA	selenocysteine synthase
Reads supporting (aligned to +/- strand):  ref base G (38/31);  new base A (2/2);  total (40/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.44e-01

AAGCCGGGATCATCGTCGGCCACAAGGACTTGATCACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTG  >  NC_002947/581618‑581905                                                                                                                                             |                                                                                                                                                    aaGCCGGGATCATCGTCGGCCACAAGGACTTGATCACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAgcgc                                                                                                                                             <  6:136977/149‑1 (MQ=255)  aGCCGGGATCATCGTCGGCCACAAGGACTTGATCACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGAGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGcc                                                                                                                                            >  8:117890/1‑149 (MQ=255)          tcatcGTCGGCCACAAGGACTTGATCACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAgc                                                                                                                                    <  6:275201/149‑1 (MQ=255)           catcGTCGGCCACAAGGACTTGATCACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAgcc                                                                                                                                   >  2:125935/1‑149 (MQ=255)             tcgtcgGCCACAAGGACTTGATCACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAgcctg                                                                                                                                 <  1:4470/149‑1 (MQ=255)             tcgtcgGCCACAAGGACTTGATCACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAgcctg                                                                                                                                 >  8:161653/1‑149 (MQ=255)              cgtcgGCCACAAGGACTTGATCACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAgcctgc                                                                                                                                >  8:281845/1‑149 (MQ=255)                tcgGCCACAAGGACTTGATCACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAgcctgcgc                                                                                                                              <  7:281845/149‑1 (MQ=255)                     cacaAGGACTTGATCACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCctgct                                                                                                                         <  1:125935/149‑1 (MQ=255)                        aaGGACTTGATCACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGAc                                                                                                                      >  7:111662/1‑149 (MQ=255)                        aaGGACTTGATCACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGAc                                                                                                                      >  8:327005/1‑149 (MQ=255)                        aaGGACTTGATCACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGAc                                                                                                                      >  7:318732/1‑149 (MQ=255)                           gACTTGATCACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCg                                                                                                                   >  6:244890/1‑149 (MQ=255)                              ttGATCACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAg                                                                                                                >  7:234915/1‑149 (MQ=255)                                  tCACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAg                                                                                                            >  7:4910/1‑149 (MQ=255)                                  tCACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGACCGAGCGCCTGCCCAGCCTGCGCctgctg                                                                                                                        >  6:263434/1‑137 (MQ=255)                                  tCACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGACCGAGCGCCTGCCCAGCCTGCGCctgctg                                                                                                                        <  5:263434/137‑1 (MQ=255)                                   cACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAgg                                                                                                           >  3:58055/1‑149 (MQ=255)                                   cACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAgg                                                                                                           >  5:127375/1‑149 (MQ=255)                                   cACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAgg                                                                                                           >  5:17871/1‑149 (MQ=255)                                     ccAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGATGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCATGcc                                                                                                         >  6:46558/1‑149 (MQ=255)                                             aaGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGAccc                                                                                                                    <  3:279329/130‑1 (MQ=255)                                             aaGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGAccc                                                                                                                    >  4:279329/1‑130 (MQ=255)                                             aaGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATcca                                                                                                 <  7:327005/149‑1 (MQ=255)                                              aGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTATCGCAACCCCGATCGCCTGGCCGAgcgc                                                                                                                                             <  1:191133/104‑1 (MQ=255)                                              aGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAgcgc                                                                                                                                             >  2:191133/1‑104 (MQ=255)                                              aGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATcgag                                                                                                <  4:52269/149‑1 (MQ=255)                                              aGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATccag                                                                                                >  3:52269/1‑149 (MQ=255)                                              aGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATccag                                                                                                <  2:22832/149‑1 (MQ=255)                                              aGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATccag                                                                                                >  1:241970/1‑149 (MQ=255)                                              aGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATccag                                                                                                >  1:22832/1‑149 (MQ=255)                                                aaaaaTCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATccaggc                                                                                              >  4:121603/1‑149 (MQ=255)                                                   aaTCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATccaggccca                                                                                           <  4:147055/149‑1 (MQ=255)                                                    aTCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGACCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATccaggcccag                                                                                          >  4:281562/1‑149 (MQ=255)                                                     tacGCTCAAGCGTGCGTTGCGCGTCGACAAGATCAACCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACACGCAGCCAGGCCGAGATcccggcgccag                                                                                         >  4:19913/3‑144 (MQ=255)                                                     tCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATccaggcccagg                                                                                         >  3:227453/1‑149 (MQ=255)                                                         cTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGCGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGa                                                                                     <  7:117890/149‑1 (MQ=255)                                                            aaGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAgc                                                                                   <  4:58055/148‑1 (MQ=255)                                                                  gcgtTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAgcgcctggc                                                                            >  1:288743/1‑149 (MQ=255)                                                                       gcgcGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGc                                                                                                                                         <  4:177761/83‑1 (MQ=255)                                                                       gcgcGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGc                                                                                                                                         >  3:177761/1‑83 (MQ=255)                                                                               aCAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAg                                                               >  1:180697/1‑149 (MQ=255)                                                                                cAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAgg                                                              >  4:206139/1‑149 (MQ=255)                                                                                cAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAgg                                                              >  7:16067/1‑149 (MQ=255)                                                                                 aaGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGgc                                                             <  8:16067/149‑1 (MQ=255)                                                                                 aaGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGgc                                                             <  8:234915/149‑1 (MQ=255)                                                                                     tCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGAGCTcaggcgccaggcgctc                                                         <  7:161653/149‑3 (MQ=255)                                                                                         ccTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCg                                                     <  8:4910/149‑1 (MQ=255)                                                                                             gctgcACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGt                                                   <  7:119583/147‑1 (MQ=255)                                                                                             gctgcACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGACCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGt                                                   <  3:281562/147‑1 (MQ=255)                                                                                              ctgcACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATcca                                                                                                 >  3:103462/1‑100 (MQ=255)                                                                                              ctgcACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATcca                                                                                                 <  4:103462/100‑1 (MQ=255)                                                                                              ctgcACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCTAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGa                                                 <  8:318732/148‑1 (MQ=255)                                                                                              ctgcACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGa                                                 <  5:244890/148‑1 (MQ=255)                                                                                              ctgcACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGa                                                 <  3:206139/148‑1 (MQ=255)                                                                                                 cACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGc                                               <  8:111662/147‑1 (MQ=255)                                                                                                        aGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCAAGGCGCGTCTCGCTGAGCAATGGGCa                                       >  8:331451/1‑148 (MQ=255)                                                                                                            gTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAg                                  >  3:84539/1‑149 (MQ=255)                                                                                                            gTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAg                                  >  2:59789/1‑149 (MQ=255)                                                                                                             tGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGc                                 <  6:17871/149‑1 (MQ=255)                                                                                                             tGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGc                                 <  6:127375/149‑1 (MQ=255)                                                                                                                  gcgcTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGc                                        <  4:269345/137‑1 (MQ=255)                                                                                                                  gcgcTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGc                                        >  3:269345/1‑137 (MQ=255)                                                                                                                      tGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGCGCGCCTGGCGCCTGAg                                                                     <  2:358450/104‑1 (MQ=255)                                                                                                                      tGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGCGCGCCTGGCGCCTGAg                                                                     >  1:358450/1‑104 (MQ=255)                                                                                                                           cGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTTGAGCCCGCGtt                   <  2:288743/149‑1 (MQ=255)                                                                                                                           cGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGtt                   >  7:76732/1‑149 (MQ=255)                                                                                                                            gCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTg                  >  1:351191/1‑149 (MQ=255)                                                                                                                             cAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTgg                 <  3:121603/149‑1 (MQ=255)                                                                                                                              aaCCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTggg                >  8:255207/1‑149 (MQ=255)                                                                                                                                ccccGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGt                    >  3:199611/1‑143 (MQ=255)                                                                                                                                ccccGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGt                    <  4:199611/143‑1 (MQ=255)                                                                                                                                    gATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCCAGCGCCTGGCGCCTGAGCTCACGGCGCGTCACGGAGAGCAATGGGCAATCACCGTAGAGCCCGCGATGGGCATGAt          >  7:302859/1‑149 (MQ=255)                                                                                                                                            ggCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTg  >  2:62119/1‑149 (MQ=255)                                                                                                                                             |                                                                                                                                                    AAGCCGGGATCATCGTCGGCCACAAGGACTTGATCACCAGGATCAAGAAAAATCCGCTCAAGCGTGCGTTGCGCGTCGACAAGATCACCCTTGCTGCACTCGAAGCGGTGCTGGCGCTGTACCGCAACCCCGATCGCCTGGCCGAGCGCCTGCCCAGCCTGCGCCTGCTGACCCGCAGCCAGGCCGAGATCCAGGCCCAGGCCGAGCGCCTGGCGCCTGAGCTCCAGGCGCGTCTCGGTGAGCAATGGGCAATCAGCGTAGAGCCCGCGTTGGGCATGATCGGCAGTG  >  NC_002947/581618‑581905

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