Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	330,555	T→C	6.1%	W365R (TGG→CGG)	betT →	choline transporter of high affinity

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	330,555	0	T	C	6.1%	204.7 / 3.1	70	W365R (TGG→CGG)	betT	choline transporter of high affinity
Reads supporting (aligned to +/- strand):  ref base T (33/30);  new base C (2/2);  total (37/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 = 9.54e-01

TGAATAACTGGACGCTCTTCTTCTGGGCATGGTGGGTGGCATGGTCGCCGTTTGTCGGCTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTC  >  NC_000913/330409‑330698                                                                                                                                                   |                                                                                                                                                tGAATAACTGGACGCTCTTCTTCTGGGCATGGTGGGTGGCATGGTCGCCGTTTGTCGGCTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATg                                                                                                                                                >  2:238211/1‑148 (MQ=255)   aaTAACTGGACGCTCTTCTTCTGGGCATGGTGGGTGGCATGGTCGCCGTTTGTCGGCTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTACGGc                                                                                                                                              >  7:1901/1‑148 (MQ=255)   aaTAACTGGACGCTCTTCTTCTGGGCATGGTGGGTGGCATGGTCGCCGTTTGTCGGCTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTACGGc                                                                                                                                              >  7:7730/1‑148 (MQ=255)                 cttcttctGGGCATGGTGGGTGGCATGGTCGCCGTTTGTCGGCTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCgg                                                                                                                                 >  3:340736/1‑147 (MQ=255)                       ctGCGCATGGTTGGTGGCATGGTCGCCGTTTGTCGGCTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAAt                                                                                                                             <  2:69337/145‑1 (MQ=255)                                   ggtggCATGGTCGCCGTTTGTCGGCTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGaa                                                                                                              <  1:238211/148‑1 (MQ=255)                                   ggtggCATGGTCGCCGTTTGTCGGCTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGa                                                                                                               >  5:353047/1‑147 (MQ=255)                                   gggggCATGGTCGCCGTGTGTCGGCTTGTGCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCGGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGaa                                                                                                              <  1:221223/145‑1 (MQ=255)                                           ggTCGCCGTTTGTCGGCTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCa                                                                                                      >  4:69997/1‑148 (MQ=255)                                                ccGTTTGTGGGCATGTTCCTGGCGCGTATCTCGCGTGGGCGGACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCAcggcg                                                                                                 <  6:20270/148‑1 (MQ=255)                                                ccGTTTGTCGGCTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTATCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCAcggcg                                                                                                 >  6:332539/1‑148 (MQ=255)                                                 cGTTTGTCGGCTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCAcggcgg                                                                                                >  3:340896/1‑148 (MQ=255)                                                    ttGTCGGCTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGTCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGgcggcgc                                                                                             >  8:346312/1‑148 (MQ=255)                                                         ggCTTGTTCCTGGCGCGCATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCAt                                                                                        <  2:335502/148‑1 (MQ=255)                                                           cTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCAttt                                                                                      >  1:89534/1‑148 (MQ=255)                                                           cTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCAACGGCGGCGCGGCAttt                                                                                      >  3:368290/1‑148 (MQ=255)                                                             tGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGc                                                                                    >  8:269579/1‑148 (MQ=255)                                                              gTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCgg                                                                                                                                        <  7:162776/95‑1 (MQ=255)                                                              gTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCgg                                                                                                                                        >  8:162776/1‑95 (MQ=255)                                                              gTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGc                                                                                    <  3:316074/147‑1 (MQ=255)                                                                   tGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCg                                                                                  >  7:281895/1‑144 (MQ=255)                                                                   tGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCg                                                                                  <  8:281895/144‑1 (MQ=255)                                                                       gcgTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTTTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCg                                                                          <  2:285838/148‑1 (MQ=255)                                                                       gcgTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCg                                                                          <  1:373796/148‑1 (MQ=255)                                                                       gcgTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCg                                                                          <  6:353047/148‑1 (MQ=255)                                                                        cgTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGa                                                                         <  8:280903/148‑1 (MQ=255)                                                                         gTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGAt                                                                        <  3:69997/148‑1 (MQ=255)                                                                           aTCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGAt                                                                        <  8:331640/146‑1 (MQ=255)                                                                           aTCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGAt                                                                        >  7:331640/1‑146 (MQ=255)                                                                            tctcGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATgg                                                                      <  2:89534/147‑1 (MQ=255)                                                                             ctcGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTAGGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTc                                                                    >  8:11978/1‑148 (MQ=255)                                                                               cgcgTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGAt                                                                        >  5:300394/1‑142 (MQ=255)                                                                               cgcgTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGAt                                                                        <  6:300394/142‑1 (MQ=255)                                                                                 cgTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATc                                                                <  1:100685/148‑1 (MQ=255)                                                                                      gCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGAt                                                                        <  1:74564/135‑1 (MQ=255)                                                                                      gCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGAt                                                                        >  2:74564/1‑135 (MQ=255)                                                                                         tACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAgcgc                                                        >  3:361713/1‑148 (MQ=255)                                                                                          aCCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATcc                                                                                                       >  2:191518/1‑100 (MQ=255)                                                                                          aCCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATcc                                                                                                       <  1:191518/100‑1 (MQ=255)                                                                                          aCCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGcgcg                                                       >  2:13785/1‑148 (MQ=255)                                                                                              ttCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCtt                                                   >  1:7215/1‑148 (MQ=255)                                                                                                cGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTc                                                  >  8:96499/1‑147 (MQ=255)                                                                                                cGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTc                                                  >  2:374574/1‑147 (MQ=255)                                                                                                 gCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTAAGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCt                                                 <  7:11978/147‑1 (MQ=255)                                                                                                   cAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCTa                                                >  5:47245/1‑146 (MQ=255)                                                                                                      ttCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGcc                                           >  3:46956/1‑148 (MQ=255)                                                                                                      ttCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGcc                                           >  5:359935/1‑148 (MQ=255)                                                                                                      ttCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGcc                                           >  5:312435/1‑148 (MQ=255)                                                                                                      ttCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGcc                                           >  7:250325/1‑148 (MQ=255)                                                                                                         gTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGCctg                                         <  6:359935/147‑1 (MQ=255)                                                                                                               ggCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGCCTGCTGgcgc                                  <  4:340736/148‑1 (MQ=255)                                                                                                                   cGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGCCTGCTGGCGCAGTa                              >  1:40838/1‑148 (MQ=255)                                                                                                                   cGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGCCTGCTGGCGCAGTa                              <  4:46956/148‑1 (MQ=255)                                                                                                                     ttGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGAGCGATGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGCAAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGCCTGCTGGCGCAGTATc                            >  5:242307/1‑148 (MQ=255)                                                                                                                       gATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGAt                                                                        >  2:94430/1‑102 (MQ=255)                                                                                                                       gATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGAt                                                                        <  1:94430/102‑1 (MQ=255)                                                                                                                       gATTATTCCGTTTACCTTCACGCTGTTACGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGCCTGCTGGCGCAGTATcc                           <  8:7730/147‑1 (MQ=255)                                                                                                                       gATTATTCCGTTTACCTTCACGCTGTTACGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGCCTGCTGGCGCAGTATcc                           <  8:1901/147‑1 (MQ=255)                                                                                                                            ttCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGtt                     <  7:269579/148‑1 (MQ=255)                                                                                                                            ttCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGt                      <  3:307470/147‑1 (MQ=255)                                                                                                                                 tttACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGGAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCt                >  7:144625/1‑148 (MQ=255)                                                                                                                                  ttACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCtt               >  5:119021/1‑148 (MQ=255)                                                                                                                                    aCCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCttt              <  1:374574/147‑1 (MQ=255)                                                                                                                                     ccTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTa             <  8:215834/147‑1 (MQ=255)                                                                                                                                     ccTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTa             <  5:60371/147‑1 (MQ=255)                                                                                                                                        tCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAgcgc         >  5:151986/1‑148 (MQ=255)                                                                                                                                        tCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAgcgc         >  5:44768/1‑148 (MQ=255)                                                                                                                                          aCGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCt       <  7:346312/148‑1 (MQ=255)                                                                                                                                           cGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTc      <  2:7215/148‑1 (MQ=255)                                                                                                                                               gTTAAGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGCTGGTCCATCCGGAGCGCGGCTTCTACAGCCTGCTCGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTc  >  4:229493/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                TGAATAACTGGACGCTCTTCTTCTGGGCATGGTGGGTGGCATGGTCGCCGTTTGTCGGCTTGTTCCTGGCGCGTATCTCGCGTGGGCGTACCATTCGCCAGTTCGTGCTGGGCACGTTGATTATTCCGTTTACCTTCACGCTGTTATGGCTCTCGGTGTTCGGCAATAGCGCGCTGTATGAAATCATCCACGGCGGCGCGGCATTTGCCGAGGAAGCGATGGTCCATCCGGAGCGCGGCTTCTACAGCCTGCTGGCGCAGTATCCGGCGTTTACCTTTAGCGCCTCCGTC  >  NC_000913/330409‑330698

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑