Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,786,944	A→G	5.5%	T259A (ACC→GCC)	tsf →	elongation factor Ts

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,786,944	0	A	G	5.5%	269.0 / 3.8	91	T259A (ACC→GCC)	tsf	elongation factor Ts
Reads supporting (aligned to +/- strand):  ref base A (43/43);  new base G (3/2);  total (46/45)
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.21e-01

CCGGAAAACATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTG  >  NC_002947/1786797‑1787089                                                                                                                                                    |                                                                                                                                                     ccGGAAAACATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCAc                                                                                                                                                     <  3:307592/149‑1 (MQ=255)   ggAAAACATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCAcc                                                                                                                                                    >  2:250889/1‑148 (MQ=255)   ggAAAACATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCAcc                                                                                                                                                    >  5:219007/1‑148 (MQ=255)   ggAAAACATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGATTCCTTCAcc                                                                                                                                                    >  7:59727/1‑148 (MQ=255)     aaaaCATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTAc                                                                                                                                                 <  8:287640/149‑1 (MQ=255)         cATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTc                                                                                                                                              >  4:168593/1‑148 (MQ=255)         cATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCa                                                                                                                                             <  2:58331/149‑1 (MQ=255)           tCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCaaa                                                                                                                                           <  4:4299/149‑1 (MQ=255)           tCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCaaa                                                                                                                                           <  5:318904/149‑1 (MQ=255)             gTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCGCCTACtttgaag                                                                                                                                          >  5:98020/1‑143 (MQ=255)                gagaACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAg                                                                                                                                       <  1:141448/148‑1 (MQ=255)                 agaACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGATGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAgg                                                                                                                                      <  5:278318/148‑1 (MQ=255)                   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                                                                                                                                  <  6:218609/131‑1 (MQ=255)                                      tCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGaa                                                                                                                                  >  5:218609/1‑131 (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                                                                                                           <  4:90761/149‑1 (MQ=255)                                         ccAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTcgac                                                                                                           >  5:319523/1‑149 (MQ=255)                                         ccAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTcgac                                                                                                           <  1:238840/149‑1 (MQ=255)                                          cAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGt                                                                                                               >  1:193865/1‑144 (MQ=255)                                          cAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGt                                                                                                               <  2:193865/144‑1 (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                                                                                                                   >  3:54803/1‑127 (MQ=255)                                                       gAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGc                                                                                                                   <  4:54803/127‑1 (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                                                                                       >  2:224943/1‑149 (MQ=255)                                                             tCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGc                                                                                       >  4:209921/1‑149 (MQ=255)                                                                 tGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCg                                                                                     >  5:175439/1‑147 (MQ=255)                                                                 tGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCg                                                                                     <  6:175439/147‑1 (MQ=255)                                                                   aaaGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCt                                                                                   >  1:263077/1‑147 (MQ=255)                                                                   aaaGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCt                                                                                   >  6:244043/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:319019/149‑1 (MQ=255)                                                                                     aTGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCa                                                               <  1:250889/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                                             >  6:324880/1‑133 (MQ=255)                                                                                                                      aaGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCt                                             <  1:137688/133‑1 (MQ=255)                                                                                                                      aaGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCt                                             >  2:137688/1‑133 (MQ=255)                                                                                                                      aaGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCt                                             <  5:324880/133‑1 (MQ=255)                                                                                                                      aaGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGc                             >  3:8393/1‑149 (MQ=255)                                                                                                                        gAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTg                           >  5:271822/1‑149 (MQ=255)                                                                                                                        gAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTg                           >  2:71974/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)                                                                                                                                                    |                                                                                                                                                     CCGGAAAACATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAA‑‑GGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGCAGTAAGACAG‑CCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTG  >  NC_002947/1786797‑1787089

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