Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,786,953	A→G	6.4%	K262E (AAA→GAA)	tsf →	elongation factor Ts

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,786,953	0	A	G	6.4%	182.7 / 2.0	62	K262E (AAA→GAA)	tsf	elongation factor Ts
Reads supporting (aligned to +/- strand):  ref base A (29/29);  new base G (2/2);  total (31/31)
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.98e-01

TCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGA  >  NC_002947/1786807‑1787101                                                                                                                                                   |                                                                                                                                                      tCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCaaa                                                                                                                                                     <  2:235518/149‑1 (MQ=255) tCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCaaa                                                                                                                                                     <  8:102274/149‑1 (MQ=255)    ttGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAg                                                                                                                                                    >  2:145729/1‑147 (MQ=255)        gaACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCg                                                                                                                                              <  1:145729/149‑1 (MQ=255)         aaCATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCg                                                                                                                                              >  8:21997/1‑148 (MQ=255)          aCATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCGCCTACTTTGAAGTAGGTGaa                                                                                                                                            <  7:51973/149‑1 (MQ=255)          aCATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGaa                                                                                                                                            <  3:76820/149‑1 (MQ=255)           ccagATCAACGGGCGTATCAACCACGTCAACGCCGCACCCTCCCCTAAACATCAAGCCTTCCTCAATAACCCGGAAATCAACGGTGGCCATCTTGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAg                                                                                                                                           <  6:46011/146‑1 (MQ=255)           cATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAATCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAg                                                                                                                                           <  3:165990/149‑1 (MQ=255)            aTGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAgg                                                                                                                                          <  5:37892/149‑1 (MQ=255)                tCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCAt                                                                                                                                       <  3:34831/148‑1 (MQ=255)                tCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCAt                                                                                                                                       <  2:166375/148‑1 (MQ=255)                     ggTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGaa                                                                                                                                 <  7:21997/149‑1 (MQ=255)                         gTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGc                                                                                                                               >  8:191090/1‑147 (MQ=255)                            tCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTc                                                                                                                          <  7:287718/149‑1 (MQ=255)                              aCCAAGTTCAAGGCCTAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTcga                                                                                                                        <  7:191090/149‑1 (MQ=255)                              aCCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTcga                                                                                                                        >  6:197637/1‑149 (MQ=255)                               ccAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTcgac                                                                                                                       >  4:23751/1‑149 (MQ=255)                                   gTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGt                                                                                                                           >  6:49836/1‑141 (MQ=255)                                   gTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGt                                                                                                                           <  5:49836/141‑1 (MQ=255)                                     tCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCaa                                                                                                                                                      >  2:61502/1‑112 (MQ=255)                                     tCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCaa                                                                                                                                                      <  1:61502/112‑1 (MQ=255)                                       aaGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCgc                                                                                                               >  7:80497/1‑149 (MQ=255)                                                gCCTCGCTGAAAGAGGAAGCCTTCGTCATGAACCCGGAAATCAAGGTCGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGt                                                                                                      <  7:50260/149‑1 (MQ=255)                                                   tCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGc                                                                                                   >  5:1838/1‑149 (MQ=255)                                                     gCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCg                                                                                                 >  2:147779/1‑149 (MQ=255)                                                         aaaGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCt                                                                                               <  3:277173/147‑1 (MQ=255)                                                         aaaGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTc                                                                                              <  6:1838/148‑1 (MQ=255)                                                              gCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTcga                                                                                                                        >  7:72366/1‑117 (MQ=255)                                                              gCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTCGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTcga                                                                                                                        <  8:72366/117‑1 (MQ=255)                                                                  gCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCgctgc                                                                                    >  7:158821/1‑149 (MQ=255)                                                                          cATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAg                                                                             >  7:140980/1‑148 (MQ=255)                                                                             gAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCgctgct                                                                                   <  1:11782/139‑1 (MQ=255)                                                                             gAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCgctgct                                                                                   >  2:11782/1‑139 (MQ=255)                                                                               aCCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTaa                                                                      <  3:23751/149‑1 (MQ=255)                                                                                 ccGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGa                                                                    >  3:112558/1‑149 (MQ=255)                                                                                  cGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGa                                                                    >  4:140777/1‑148 (MQ=255)                                                                                       aTCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTAACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTAGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTG‑CACGCTAGTAAGACAGccc                                                              >  1:260540/1‑149 (MQ=255)                                                                                            ggTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCt                                                         <  8:140980/149‑1 (MQ=255)                                                                                              ttGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGt                                                       >  8:277032/1‑149 (MQ=255)                                                                                                       ttgCCAAGAAAGCCGGTGCATAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCaaa                                              <  3:5987/147‑1 (MQ=255)                                                                                                           cAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAgg                                          >  5:227275/1‑149 (MQ=255)                                                                                                              gAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAg                                                                                                                                           <  1:240812/50‑1 (MQ=255)                                                                                                              gAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAg                                                                                                                                           >  2:240812/1‑50 (MQ=255)                                                                                                              gAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTg                                       >  7:109835/1‑149 (MQ=255)                                                                                                                       tGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTc                                <  8:69571/147‑1 (MQ=255)                                                                                                                       tGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTc                                >  7:69571/1‑147 (MQ=255)                                                                                                                          aGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCAcgcg                           >  1:192554/1‑149 (MQ=255)                                                                                                                            aaaTCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCt                                                         <  3:276968/117‑1 (MQ=255)                                                                                                                            aaaTCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCACGTCt                                                         >  4:276968/1‑117 (MQ=255)                                                                                                                             aaTCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACgcgcgc                        <  1:147779/149‑1 (MQ=255)                                                                                                                                 gatgccTTCGCCTACTTTGAAGTAGGTGAAGGAAACGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGcc                    >  6:20265/5‑149 (MQ=38)                                                                                                                                 gTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGc                     <  5:197637/148‑1 (MQ=255)                                                                                                                                  tttCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCct                   >  5:201720/1‑149 (MQ=255)                                                                                                                                     ccTTCGCCTACTTTGAAGTAGGTGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTc                                >  6:289932/1‑133 (MQ=255)                                                                                                                                     ccTTCGCCTACTTTGAAGTAGGTGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTc                                <  5:289932/133‑1 (MQ=255)                                                                                                                                       ttCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGc                     <  3:178472/142‑1 (MQ=255)                                                                                                                                       ttCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGc                     >  4:178472/1‑142 (MQ=255)                                                                                                                                          aCCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTgggg           <  6:201720/149‑4 (MQ=255)                                                                                                                                               cTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACg      >  5:52613/1‑149 (MQ=255)                                                                                                                                                ttCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑a                                                                          <  3:158655/80‑1 (MQ=255)                                                                                                                                                ttCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑a                                                                          >  4:158655/1‑80 (MQ=255)                                                                                                                                                ttCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACgg     >  3:9738/1‑149 (MQ=255)                                                                                                                                                ttCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCACAAAGAGGCTGCGCGCGTGAGCGGGCAGCctc                  >  2:30693/1‑136 (MQ=38)                                                                                                                                                ttCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCACAAAGAGGCTGCGCGCGTGAGCGGGCAGCctc                  <  1:30693/136‑1 (MQ=38)                                                                                                                                                  cAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACggg    <  4:112558/148‑1 (MQ=255)                                                                                                                                                   aaaGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGgga  <  6:227275/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                      TCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGGGA  >  NC_002947/1786807‑1787101

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