Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,786,944	A→G	8.0%	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	8.0%	172.8 / 5.2	62	T259A (ACC→GCC)	tsf	elongation factor Ts
Reads supporting (aligned to +/- strand):  ref base A (28/29);  new base G (3/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 = 8.84e-01

GAAAACATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCA  >  NC_002947/1786800‑1787092                                                                                                                                                 |                                                                                                                                                      gAAAACATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTa                                                                                                                                                   >  1:235518/1‑149 (MQ=255)   aaaCATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACt                                                                                                                                                 <  3:188791/149‑1 (MQ=255)      cATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTc                                                                                                                                               >  4:34831/1‑148 (MQ=255)        tCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCaaa                                                                                                                                            <  8:102274/149‑1 (MQ=255)        tCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCaaa                                                                                                                                            <  2:235518/149‑1 (MQ=255)           ttGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCGCCTACtttgaag                                                                                                                                           >  8:51973/1‑142 (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                                                                                                                              <  2:166375/148‑1 (MQ=255)                       tCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCAt                                                                                                                              <  3:34831/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                                                                                                                  <  5:49836/141‑1 (MQ=255)                                          gTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGt                                                                                                                  >  6:49836/1‑141 (MQ=255)                                            tCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCaa                                                                                                                                             <  1:61502/112‑1 (MQ=255)                                            tCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCaa                                                                                                                                             >  2:61502/1‑112 (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                                                                                                                                  >  2:240812/1‑50 (MQ=255)                                                                                                                     gAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAg                                                                                                                                  <  1:240812/50‑1 (MQ=255)                                                                                                                     gAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTg                              >  7:109835/1‑149 (MQ=255)                                                                                                                              tGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTc                       >  7:69571/1‑147 (MQ=255)                                                                                                                              tGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTc                       <  8:69571/147‑1 (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            >  4:178472/1‑142 (MQ=255)                                                                                                                                              ttCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGc            <  3:178472/142‑1 (MQ=255)                                                                                                                                                 aCCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTgggg  <  6:201720/149‑4 (MQ=255)                                                                                                                                                 |                                                                                                                                                      GAAAACATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCA  >  NC_002947/1786800‑1787092

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