Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,786,952	C→T	6.1%	F261F (TTC→TTT)	tsf →	elongation factor Ts

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,786,952	0	C	T	6.1%	193.1 / 2.2	65	F261F (TTC→TTT)	tsf	elongation factor Ts
Reads supporting (aligned to +/- strand):  ref base C (30/31);  new base T (2/2);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.17e-01

CATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGG  >  NC_002947/1786805‑1787099                                                                                                                                                    |                                                                                                                                                     cATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTc                                                                                                                                                      >  4:34831/1‑148 (MQ=255)   tCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCaaa                                                                                                                                                   <  8:102274/149‑1 (MQ=255)   tCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCaaa                                                                                                                                                   <  2:235518/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                                                                                                                         <  5:49836/141‑1 (MQ=255)                                     gTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGt                                                                                                                         >  6:49836/1‑141 (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                                                                                 >  2:11782/1‑139 (MQ=255)                                                                               gAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCgctgct                                                                                 <  1:11782/139‑1 (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                   >  4:178472/1‑142 (MQ=255)                                                                                                                                         ttCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGc                   <  3:178472/142‑1 (MQ=255)                                                                                                                                            aCCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTgggg         <  6:201720/149‑4 (MQ=255)                                                                                                                                                 cTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACg    >  5:52613/1‑149 (MQ=255)                                                                                                                                                  ttCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑a                                                                        >  4:158655/1‑80 (MQ=255)                                                                                                                                                  ttCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑a                                                                        <  3:158655/80‑1 (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)                                                                                                                                                    |                                                                                                                                                     CATCGTTGAGAACATGATCAACGGTCGTATCACCAAGTTCAAGGCCGAAGCCTCGCTGAAAGAGCAAGCCTTCGTCATGAACCCGGAAATCAAGGTTGGCGAGCTGGCCAAGAAAGCCGGTGCAGAAATCGTTTCCTTCACCTACTTCAAAGTAGGCGAAGGCATCGAGAAGCCAGTCGACGACTTCGCTGCTGAAGTTGCCGCTCAGGTCGCTGCTGCCAAGC‑AGTAAGACAGCCCCGTCTGTCGCCCCAAAGAGGCTGCCCGCTCACGCGCGCAGCCTCTTTGTCAAAACGGG  >  NC_002947/1786805‑1787099

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