Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,786,961	C→T	5.4%	G264G (GGC→GGT)	tsf →	elongation factor Ts

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,786,961	0	C	T	5.4%	279.0 / 3.4	92	G264G (GGC→GGT)	tsf	elongation factor Ts
Reads supporting (aligned to +/- strand):  ref base C (47/40);  new base T (2/3);  total (49/43)
Fisher's exact test for biased strand distribution p-value = 6.62e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.09e-01

AACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTA  >  NC_002947/1786815‑1787109                                                                                                                                                   |                                                                                                                                                        aaCATGATCAACGGTCGCATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCg                                                                                                                                                        >  8:298717/1‑148 (MQ=255)   cagtaTCCACGGGCGTATCACCAAGTTCCAGGCCGAAGCCGCGCTTAAAGAGCAACCCTTCGTCATGAACCCGTAAATCAAGGTTGGCGAGCTGGGCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑g                                                                                                                                                   <  8:211536/145‑1 (MQ=255)   cATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCGCCTACTTTGAAGTAGGTGAA‑‑g                                                                                                                                                   <  6:113947/149‑1 (MQ=255)   cATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑g                                                                                                                                                   <  7:103361/149‑1 (MQ=255)   cATGATCAACGGTCGCATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑g                                                                                                                                                   <  7:298717/149‑1 (MQ=255)    aTGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑gg                                                                                                                                                  <  7:127799/149‑1 (MQ=255)     tGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCGCCTACTTTGAAGTAGGTGAA‑‑gga                                                                                                                                                 <  6:98020/149‑2 (MQ=255)        tCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATc                                                                                                                                              <  3:25178/149‑1 (MQ=255)             ggTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGaa                                                                                                                                         <  4:151569/149‑1 (MQ=255)                    tCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGaa                                                                                                                                                      >  5:218609/1‑131 (MQ=255)                    tCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGaa                                                                                                                                                      <  6:218609/131‑1 (MQ=255)                    tCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTc                                                                                                                                  >  7:249455/1‑149 (MQ=255)                    tCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTc                                                                                                                                  >  7:33770/1‑149 (MQ=255)                    tCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATAAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑g                                                                                                                                                   >  6:237460/1‑132 (MQ=255)                    tCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATAAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑g                                                                                                                                                   <  5:237460/132‑1 (MQ=255)                       ccAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTcgac                                                                                                                               >  5:319523/1‑149 (MQ=255)                       ccAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTcgac                                                                                                                               <  1:238840/149‑1 (MQ=255)                       ccAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTcgac                                                                                                                               <  4:90761/149‑1 (MQ=255)                        cAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGt                                                                                                                                   <  2:193865/144‑1 (MQ=255)                        cAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGt                                                                                                                                   >  1:193865/1‑144 (MQ=255)                          aGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCgacgac                                                                                                                            >  1:158460/1‑149 (MQ=255)                           gTTCAAGGCCGAAGCCTCGCCGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACt                                                                                                                           >  1:16433/1‑149 (MQ=255)                               aaGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTcgacga                                                                                                                             >  3:288769/1‑143 (MQ=255)                               aaGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTcgacga                                                                                                                             <  4:288769/143‑1 (MQ=255)                               aaGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCg                                                                                                                        <  8:59727/148‑1 (MQ=255)                                    cGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCgctgct                                                                                                                   >  2:272423/1‑148 (MQ=255)                                     gAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGc                                                                                                                                       <  4:54803/127‑1 (MQ=255)                                     gAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGc                                                                                                                                       >  3:54803/1‑127 (MQ=255)                                      aaGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGc                                                                                                                                       >  6:3211/1‑126 (MQ=255)                                      aaGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGc                                                                                                                                       <  5:3211/126‑1 (MQ=255)                                          cTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGGA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTg                                                                                                            >  1:53086/1‑149 (MQ=255)                                           tCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGc                                                                                                           >  4:209921/1‑149 (MQ=255)                                           tCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGc                                                                                                           >  2:224943/1‑149 (MQ=255)                                               tGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCg                                                                                                         >  5:175439/1‑147 (MQ=255)                                               tGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCg                                                                                                         <  6:175439/147‑1 (MQ=255)                                                 aaaGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCt                                                                                                       >  6:244043/1‑147 (MQ=255)                                                 aaaGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCt                                                                                                       >  1:263077/1‑147 (MQ=255)                                                 aaaGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTc                                                                                                      <  3:45788/148‑1 (MQ=255)                                                       cAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCgc                                                                                               >  8:89767/1‑149 (MQ=255)                                                          gCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATTGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCgctgc                                                                                            >  4:287396/1‑149 (MQ=255)                                                           ccTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCgctgct                                                                                           <  6:219007/149‑1 (MQ=255)                                                                gTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGcc                                                                                        >  2:319019/1‑147 (MQ=255)                                                                 tCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTTCCGCTCAGGTCGCTGCTGCCAAg                                                                                     <  2:263077/149‑1 (MQ=255)                                                                 tCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAg                                                                                     <  7:167896/149‑1 (MQ=255)                                                                  cATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAg                                                                                     >  7:81351/1‑148 (MQ=255)                                                                   aTGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCa                                                                                   <  1:250889/149‑1 (MQ=255)                                                                   aTGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCa                                                                                   <  1:319019/149‑1 (MQ=255)                                                                    tGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAg                                                                                  <  8:33770/149‑1 (MQ=255)                                                                     gAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAg                                                                                  >  6:198571/1‑148 (MQ=255)                                                                      aaCCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGt                                                                                 <  6:319523/148‑1 (MQ=255)                                                                       aCCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTaa                                                                               <  2:158460/149‑1 (MQ=255)                                                                       aCCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCGGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGt                                                                                 >  8:208824/1‑147 (MQ=255)                                                                         ccGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGa                                                                             >  1:35072/1‑149 (MQ=255)                                                                                   aGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCgct                                                                                                                      <  7:40306/98‑1 (MQ=255)                                                                                   aGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCgct                                                                                                                      >  8:40306/1‑98 (MQ=255)                                                                                    ggTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCt                                                                 >  5:323932/1‑149 (MQ=255)                                                                                         gCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGc                                                            <  4:107887/149‑1 (MQ=255)                                                                                             gCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGcccc                                                         >  1:283577/1‑148 (MQ=255)                                                                                                 gCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAag                                                     <  2:35072/148‑1 (MQ=255)                                                                                                   cAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAgg                                                  >  1:146768/1‑149 (MQ=255)                                                                                                    aaGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCt                                                                 <  1:137688/133‑1 (MQ=255)                                                                                                    aaGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCt                                                                 <  5:324880/133‑1 (MQ=255)                                                                                                    aaGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCt                                                                 >  2:137688/1‑133 (MQ=255)                                                                                                    aaGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCt                                                                 >  6:324880/1‑133 (MQ=255)                                                                                                    aaGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGc                                                 >  3:8393/1‑149 (MQ=255)                                                                                                      gAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTg                                               >  2:71974/1‑149 (MQ=255)                                                                                                      gAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTg                                               >  5:271822/1‑149 (MQ=255)                                                                                                       aaaGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGc                                              <  5:244043/149‑1 (MQ=255)                                                                                                          gcgggTGCTCAAATCTTTGCATTCACCTACTTAAACGTAGGCTAA‑‑GGCATCGTTAAGCCAGTCTACTACTGGGCTGCGTAATTTGCCGCTCAGGTCGCTGCTGCAATGCTGT‑AGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGCCCGc                                          <  2:16433/146‑1 (MQ=255)                                                                                                               tGCAGAAATAGTTTCCTTCACCTACTTCAAAGTATGCGAATCTTCCATGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTgcgcgcgt                                        <  1:272423/149‑3 (MQ=255)                                                                                                                  aGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCAcgcg                                   >  1:214478/1‑149 (MQ=255)                                                                                                                   gAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCAcgcgc                                  >  1:325430/1‑149 (MQ=255)                                                                                                                         gTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGc                             <  8:249455/148‑1 (MQ=255)                                                                                                                              cTTCGCCTACTTTGAAGTAGGTGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCttt                       >  4:204002/1‑149 (MQ=255)                                                                                                                              cTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCttt                       >  8:212708/1‑149 (MQ=255)                                                                                                                              cTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCttt                       >  5:104295/1‑149 (MQ=255)                                                                                                                              cTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCtt                        <  4:8393/148‑1 (MQ=255)                                                                                                                               ttCGCCTACTTTGAAGTAGGTGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTg                      >  5:228809/1‑149 (MQ=255)                                                                                                                                 cACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAGCCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCttt                       <  1:71974/147‑1 (MQ=255)                                                                                                                                   ccTACTTCAAAGTAGGCGAA‑‑CGCATCGAGAACCCAGTCGACGACTTCGCAGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCACTAAGACAG‑CCCCGACTGTCGCCCCAAAGAGGCTGCCCGCTCACGCCAGCAGCCTCATTGTCa                   >  4:67913/1‑148 (MQ=255)                                                                                                                                    cTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCAcgc                                    >  7:16391/1‑130 (MQ=255)                                                                                                                                    cTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCAcgc                                    <  8:16391/130‑1 (MQ=255)                                                                                                                                    cTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCaa                  >  4:220812/1‑148 (MQ=255)                                                                                                                                       cTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGACTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACg              >  2:31594/1‑149 (MQ=255)                                                                                                                                          tgaaGTAGGTGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACgg             <  6:228809/145‑1 (MQ=255)                                                                                                                                            aaGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTg                      >  5:288834/1‑136 (MQ=255)                                                                                                                                            aaGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTg                      <  6:288834/136‑1 (MQ=255)                                                                                                                                             aGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGggag         >  5:66666/1‑148 (MQ=255)                                                                                                                                                 ggCGAA‑‑GGCATCGAGAAGCCAGTCGACGAGTTCGCTGCTGAAGTTGCCGCTGAGGTCGCTGCTGCCAAGCAGTAAGACTG‑CGCCGTCTGTCGCGCCAAAGAGGCTGCCCGCTCACGCGCGGAGCCTCTTTGTCAAAACGGGGAGGAGt     <  3:67913/148‑1 (MQ=255)                                                                                                                                                  gCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGttt   >  5:326002/1‑149 (MQ=255)                                                                                                                                                  gCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGttt   <  5:198571/149‑1 (MQ=255)                                                                                                                                                  gCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGtt    >  8:182600/1‑148 (MQ=255)                                                                                                                                                   cGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTta  <  1:224943/149‑1 (MQ=255)                                                                                                                                                   cGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTta  <  7:89767/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                        AACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGAGGAGTTTA  >  NC_002947/1786815‑1787109

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