Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,471,398	A→C	5.2%	S477A (TCG→GCG)	alg8 ←	glycosyltransferase Alg8

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,471,398	0	A	C	5.2%	311.4 / 7.0	97	S477A (TCG→GCG)	alg8	glycosyltransferase Alg8
Reads supporting (aligned to +/- strand):  ref base A (43/49);  new base C (2/3);  total (45/52)
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.72e-01

CTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCA  >  NC_002947/1471256‑1471541                                                                                                                                                |                                                                                                                                                 cTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCa                                                                                                                                            <  4:167568/149‑1 (MQ=255)  ttCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAg                                                                                                                                           >  4:121893/1‑149 (MQ=255)  ttCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAg                                                                                                                                           >  1:59928/1‑149 (MQ=255)     ggACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCTCGCGGACCatga                                                                                                                                         <  1:137044/148‑4 (MQ=255)     ggACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGgtg                                                                                                                                        <  5:138553/149‑1 (MQ=255)     ggACTCATGCACGACATTGACGTTCACGGCGGTATTCATGCTGG‑TTTGTTCCTGTTAGATCCGCTCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCTCGCGGACCatga                                                                                                                                         <  5:361016/148‑4 (MQ=255)       aCTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGtt                                                                                                                                      >  5:248942/1‑149 (MQ=255)        cTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTg                                                                                                                                     <  5:44989/149‑1 (MQ=255)         tCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTg                                                                                                                                     >  1:350020/1‑148 (MQ=255)           aTGCACGACATTGACGTTCACGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGa                                                                                                                                    <  5:153930/147‑1 (MQ=255)           aTGCACGACATTGACGTTCACGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAAc                                                                                                                                  >  3:129746/1‑149 (MQ=255)           aTGCACGACAATGACGTTCACGGCGGAATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCAGGTTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGAAGCCGAGAAGGTCATGGTCCGAGAGGACCAGGTGTTCaa                                                                                                                                   >  4:112254/1‑148 (MQ=255)            tGCACGACATTGACGTTCACGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAAc                                                                                                                                  >  5:142307/1‑148 (MQ=255)              cACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTTGTTCCTGTTAGATCTGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCa                                                                                                                                >  1:60647/1‑149 (MQ=255)               aCGACATTGACGTTCACGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGGCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCAtt                                                                                                                              >  5:357501/1‑149 (MQ=255)               aCGACATTGACGTTCACGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCAtt                                                                                                                              <  7:198942/149‑1 (MQ=255)               aCGACATTGACGTTCACGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCAtt                                                                                                                              >  4:47934/1‑149 (MQ=255)                  aCATTGACGTTCACGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCAttgtt                                                                                                                           <  8:117344/149‑1 (MQ=255)                    aTTGACGTTCACGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCAttgttg                                                                                                                          >  1:213462/1‑148 (MQ=255)                     ttGACGTTCACGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGaa                                                                                                                        <  3:98602/149‑1 (MQ=255)                       gACGTTCACGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGaaa                                                                                                                       <  4:129746/148‑1 (MQ=255)                        aCGTTCACGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCGTGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGc                                                                                                                     >  8:302291/1‑149 (MQ=255)                         cGTTCACGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCTGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCt                                                                                                                    <  2:60647/149‑1 (MQ=255)                          gTTCACGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTg                                                                                                                   >  2:21865/1‑149 (MQ=255)                           ttCACGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTgg                                                                                                                  <  6:142307/149‑1 (MQ=255)                              aCGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGc                                                                                                                 >  1:310343/1‑147 (MQ=255)                              aCGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGa                                                                                                               <  3:401616/149‑1 (MQ=255)                               cggcggTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAg                                                                                                              <  3:47934/149‑1 (MQ=255)                                 gcggTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGt                                                                                                            >  8:99636/1‑149 (MQ=255)                                          aTGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGa                                                                                                   <  5:48746/149‑1 (MQ=255)                                          aTGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGa                                                                                                   <  7:99636/149‑1 (MQ=255)                                          aTGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGa                                                                                                   <  3:121893/149‑1 (MQ=255)                                          aTGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGa                                                                                                   <  3:337561/149‑1 (MQ=255)                                            ggtgg‑tTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGaaa                                                                                                                       >  2:376059/1‑127 (MQ=255)                                            ggtgg‑tTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGaaa                                                                                                                       <  1:376059/127‑1 (MQ=255)                                                    tGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCGTGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTg                                                                                          <  7:302291/148‑1 (MQ=255)                                                            ttAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAg                                                                                                  <  6:212604/133‑1 (MQ=255)                                                            ttAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAg                                                                                                  >  5:212604/1‑133 (MQ=255)                                                            ttAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGAc                                                                                   <  2:59928/147‑1 (MQ=255)                                                               gatccgatccGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAg                                                                                                              <  4:349018/118‑1 (MQ=255)                                                               gatccgatccGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAg                                                                                                              >  3:349018/1‑118 (MQ=255)                                                               gatccgatccGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTc                                                                              >  2:274165/1‑149 (MQ=255)                                                                 tccgatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCa                                                                            >  8:144698/1‑149 (MQ=255)                                                                  ccgatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGa                                                                                                   <  1:357083/126‑1 (MQ=255)                                                                  ccgatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGa                                                                                                   >  2:357083/1‑126 (MQ=255)                                                                   cgatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCAc                                                                           <  6:357501/148‑1 (MQ=255)                                                                   cgatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCAc                                                                           <  2:213462/148‑1 (MQ=255)                                                                    gatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGTG‑AAGTGGGCTGACGCGTCCACGa                                                                         >  6:119371/1‑149 (MQ=255)                                                                         gggTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCt                                                                    <  1:21865/149‑1 (MQ=255)                                                                          ggTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCtt                                                                   <  7:144698/149‑1 (MQ=255)                                                                                  tCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAg                                                                                                  <  5:84819/111‑1 (MQ=255)                                                                                  tCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAg                                                                                                  >  6:84819/1‑111 (MQ=255)                                                                                  tCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGc                                                           >  6:387531/1‑149 (MQ=255)                                                                                      acGACCATGAACAGCACAGCAACTAAGATGATGGCAGCCTAGAAGGTCATGTTCCGCTAGGACCAGGTGTTGAACCATTGTTGAAAGCTTGCGAGGTCACGCTTTAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGgaa                                                       <  3:137285/149‑1 (MQ=255)                                                                                      acGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAg                                                                                                  >  7:354902/1‑107 (MQ=255)                                                                                      acGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAg                                                                                                  <  8:354902/107‑1 (MQ=255)                                                                                      acGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGc                                                           >  2:253273/1‑145 (MQ=255)                                                                                      acGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGc                                                           <  1:253273/145‑1 (MQ=255)                                                                                        gACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAgg                                                            <  4:322175/142‑1 (MQ=255)                                                                                        gACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAgg                                                            >  3:322175/1‑142 (MQ=255)                                                                                        gACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGgaaga                                                     <  2:350020/149‑1 (MQ=255)                                                                                        gACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGgaag                                                      >  3:44200/1‑148 (MQ=255)                                                                                        gACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGgaag                                                      <  4:44200/148‑1 (MQ=255)                                                                                             tGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACg                                                >  4:291854/1‑149 (MQ=255)                                                                                                 cagcacagcaACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACg                                                >  5:52754/1‑145 (MQ=255)                                                                                                   gcacagcaACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATc                                          >  1:189145/1‑149 (MQ=255)                                                                                                    cacagcaACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCt                                         <  3:291854/149‑1 (MQ=255)                                                                                                      cagcaACGAAGATGCTGGCAGCCTAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGTCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTtc                                       <  5:119371/149‑1 (MQ=255)                                                                                                       agcaACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCtt                                                                   >  4:361809/1‑120 (MQ=255)                                                                                                       agcaACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCtt                                                                   <  3:361809/120‑1 (MQ=255)                                                                                                                   gCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCga                          >  1:97605/1‑149 (MQ=255)                                                                                                                       gCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑Ca                                                                                             >  7:146120/1‑78 (MQ=255)                                                                                                                       gCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑Ca                                                                                             <  8:146120/78‑1 (MQ=255)                                                                                                                        cAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATc                     <  2:310343/149‑1 (MQ=255)                                                                                                                           ccGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTgg                  >  2:354197/1‑149 (MQ=255)                                                                                                                            cGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTtc                                       >  3:108699/1‑127 (MQ=255)                                                                                                                            cGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTtc                                       <  4:108699/127‑1 (MQ=255)                                                                                                                            cGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGt                 >  1:205037/1‑149 (MQ=255)                                                                                                                             gagaAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGtt                <  6:52754/149‑1 (MQ=255)                                                                                                                                    acctGGTCCTCGCGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGgaa                                                       >  1:171020/4‑103 (MQ=255)                                                                                                                                    acctGGTCCTCGCGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGgaa                                                       <  2:171020/100‑1 (MQ=255)                                                                                                                                     cctGGTCCTCGCGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGAGAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATAAGCGCGCCGATGATCTGGTTGTaa            >  1:328049/3‑146 (MQ=255)                                                                                                                                      aTGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTgg                  >  3:295053/1‑138 (MQ=255)                                                                                                                                      aTGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTgg                  <  4:295053/138‑1 (MQ=255)                                                                                                                                      aTGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATaga        <  2:205037/148‑1 (MQ=255)                                                                                                                                         gTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCTAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAAt    <  2:189145/149‑1 (MQ=255)                                                                                                                                         gTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAAt    <  2:97605/149‑1 (MQ=255)                                                                                                                                          tCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAAt    >  4:185489/1‑148 (MQ=255)                                                                                                                                          tCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAAt    >  2:309348/1‑148 (MQ=255)                                                                                                                                            cgcgAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAacac                                                 >  4:158459/1‑101 (MQ=255)                                                                                                                                            cgcgAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAacac                                                 <  3:158459/101‑1 (MQ=255)                                                                                                                                            cgcgAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCa  <  5:387531/148‑1 (MQ=255)                                                                                                                                            cgcgAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCa  <  6:248942/148‑1 (MQ=255)                                                                                                                                               gAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGgaa                                                       <  1:332032/92‑1 (MQ=255)                                                                                                                                               gAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGgaa                                                       >  2:332032/1‑92 (MQ=255)                                                                                                                                               gAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTtc                                       <  1:349259/108‑1 (MQ=255)                                                                                                                                               gAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCGTGTCGAGGCGGAAGAACACGTAGATCTtc                                       >  2:349259/1‑108 (MQ=255)                                                                                                                                                |                                                                                                                                                 CTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGG‑TTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTCCGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGG‑CAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCA  >  NC_002947/1471256‑1471541

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