Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,471,405	A→C	7.6%	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,405	0	A	C	7.6%	179.4 / 5.9	68	S477A (TCG→GCG)	alg8	glycosyltransferase Alg8
Reads supporting (aligned to +/- strand):  ref base A (29/33);  new base C (3/2);  total (33/35)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.42e-01

GCTGGGCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCG  >  NC_002947/1471257‑1471551                                                                                                                                                      |                                                                                                                                                   gcTGGGCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAg                                                                                                                                                   <  1:1216512/150‑1 (MQ=255)   tGGGCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGa                                                                                                                                                 <  1:1074522/150‑1 (MQ=255)     ggCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGAcc                                                                                                                                               <  2:630453/150‑1 (MQ=255)      gCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGa                                                                                                                                                    <  1:1054544/144‑1 (MQ=255)      gCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGa                                                                                                                                                    >  2:1054544/1‑144 (MQ=255)                    cACGACATTGACGTTCACGGCGGTATTNATGGTGGTTTGTNCCTGTTNGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCAtt                                                                                                                                >  1:1171170/1‑150 (MQ=255)                    cACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGCGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCAtt                                                                                                                                <  2:742428/150‑1 (MQ=255)                            tGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAg                                                                                                                        >  2:37540/1‑150 (MQ=255)                                   cACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGa                                                                                                                 <  2:1171170/150‑1 (MQ=255)                                   cACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGa                                                                                                                 >  2:84540/1‑150 (MQ=255)                                      ggcggTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGt                                                                                                              <  1:37540/150‑1 (MQ=255)                                          gTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACg                                                                                                          >  1:43354/1‑150 (MQ=255)                                             ttCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCtt                                                                                                       <  1:1111761/150‑1 (MQ=255)                                               cATGGTGGTTTGTTCCTGTTAGATACGTTCCGGGTTGGGCTCACACGAACATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGa                                                                                                     <  2:707547/150‑1 (MQ=255)                                                aTGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAg                                                                                                                                             <  2:1026211/109‑1 (MQ=255)                                                aTGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAg                                                                                                                                             >  1:1026211/1‑109 (MQ=255)                                                 tggtggTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCAttgt                                                                                                                              <  1:319029/123‑1 (MQ=255)                                                 tggtggTTTGTTCCTGTTAGATCCGATCCGGGTTGAGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAAAAGCACACGTC‑TGAACTCCAGTCAACCAGTGTTTGAACCATTGTTGAAAGCTGGCGAGGGCACGCTTGAgg                                                                                                   >  1:1115209/1‑150 (MQ=255)                                                 cggtggtTTGGTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCAttgt                                                                                                                              >  2:319029/2‑123 (MQ=255)                                                  ggtggtTTGTTCCTGTTAGATCCGATCCGGGTTGGGCCCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAggg                                                                                                  <  1:1035791/150‑1 (MQ=255)                                                    tggtTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGc                                                                                                                   <  1:142882/131‑1 (MQ=255)                                                    tggtTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGc                                                                                                                   >  2:142882/1‑131 (MQ=255)                                                    tggtTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCa                                                                                                >  1:930356/1‑150 (MQ=255)                                                       tttGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTg                                                                                             >  1:141822/1‑150 (MQ=255)                                                         tGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTggg                                                                                           <  1:168616/150‑1 (MQ=255)                                                           ttCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCt                                                                                         <  2:930356/150‑1 (MQ=255)                                                              cTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGAc                                                                                      >  2:829063/1‑150 (MQ=255)                                                               tGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGAcg                                                                                     >  1:528066/1‑150 (MQ=255)                                                                gTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGAcgc                                                                                    <  1:829063/150‑1 (MQ=255)                                                                   aGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTc                                                                                 <  2:528066/150‑1 (MQ=255)                                                                     atccgatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGgtg                                                                                                                                          <  2:560807/91‑1 (MQ=255)                                                                     atccgatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGgtg                                                                                                                                          >  1:560807/1‑91 (MQ=255)                                                                      tccgatacgGGTTGAGGATAGCCGACAATGAACAGCACAGCTAGGAAGAGGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCAc                                                                              >  2:338467/1‑150 (MQ=255)                                                                       ccgatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACg                                                                             >  2:1178757/1‑150 (MQ=255)                                                                        cgatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGa                                                                            <  1:84540/150‑1 (MQ=255)                                                                         gatccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGAc                                                                           >  2:776333/1‑150 (MQ=255)                                                                          atccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACt                                                                          >  2:243794/1‑150 (MQ=255)                                                                            ccgGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGc                                                                        >  1:961592/1‑150 (MQ=255)                                                                              gggTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAg                                                                                               <  1:560510/125‑1 (MQ=255)                                                                              gggTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAg                                                                                               >  2:560510/1‑125 (MQ=255)                                                                              gggTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCtt                                                                      <  2:141822/150‑1 (MQ=255)                                                                                          cacGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGCCGCGTCCACGACTGCTTGTCGAGGCGgaa                                                          <  2:716976/150‑1 (MQ=255)                                                                                          cacGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGgaa                                                          >  1:637115/1‑150 (MQ=255)                                                                                             gACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCAttgtt                                                                                                                             >  1:869347/1‑80 (MQ=255)                                                                                             gACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCAttgtt                                                                                                                             <  2:869347/80‑1 (MQ=255)                                                                                              aCCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAac                                                      <  2:1115209/150‑1 (MQ=255)                                                                                                 aTGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACg                                                   >  1:775779/1‑150 (MQ=255)                                                                                                    aaCAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCtt                                                                      >  1:1168708/1‑128 (MQ=255)                                                                                                    aaCAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCtt                                                                      <  2:1168708/128‑1 (MQ=255)                                                                                                               aaCGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAg                                     >  1:494698/1‑150 (MQ=255)                                                                                                                    aGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAgcgcgc                                <  2:775779/150‑1 (MQ=255)                                                                                                                     gATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGcc                               <  1:1178757/150‑1 (MQ=255)                                                                                                                       tGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAggg                                                                                                  <  1:530319/81‑1 (MQ=255)                                                                                                                       tGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAggg                                                                                                  >  2:530319/1‑81 (MQ=255)                                                                                                                        gCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCgat                            <  2:494698/150‑1 (MQ=255)                                                                                                                           ggCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCgatgat                         <  1:338467/150‑1 (MQ=255)                                                                                                                               gCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTgg                     <  1:243794/150‑1 (MQ=255)                                                                                                                                 cGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGtt                   >  1:1264056/1‑150 (MQ=255)                                                                                                                                  gagaAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTAAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTg                  >  2:930222/1‑150 (MQ=255)                                                                                                                                      aacacctGGTC‑CTCGCGGACCATGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTtcat                                        <  2:354391/118‑1 (MQ=255)                                                                                                                                      aacacctGGTC‑CTCGCGGACCATGTGGTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTtcat                                        >  1:354391/7‑124 (MQ=255)                                                                                                                                        gTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGTAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTACATCTTCATCAGCGCGCCGATGATCTGGTTGTAATag            >  2:743389/1‑150 (MQ=255)                                                                                                                                         tCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATaga           <  2:43354/150‑1 (MQ=255)                                                                                                                                           ctGGTC‑CTCGCGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGAc                                                                           <  1:709839/83‑1 (MQ=255)                                                                                                                                           ctGGTC‑CTCGCGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGAc                                                                           >  2:709839/2‑84 (MQ=255)                                                                                                                                           aTGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAgaga         <  1:776333/150‑1 (MQ=255)                                                                                                                                            tGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGaa        >  2:1248317/1‑150 (MQ=255)                                                                                                                                             ggTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAAt       <  1:1248317/150‑1 (MQ=255)                                                                                                                                              gTC‑CTCGCGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGATCACGCATGAGGACAGTGGGCTGACACGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATATTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATc      >  2:68335/1‑150 (MQ=255)                                                                                                                                                   gcgAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCg  >  1:628609/1‑150 (MQ=255)                                                                                                                                                      |                                                                                                                                                   GCTGGGCTTCGGACTCATGCACGACATTGACGTTCACGGCGGTATTCATGGTGGTTTGTTCCTGTTAGATCCGATCCGGGTTGGGCTCACACGACCATGAACAGCACAGCAACGAAGATGCTGGCAGCCGAGAAGGTCATGGTC‑CGCGAGGACCAGGTGTTGAACCATTGTTGAAAGCTGGCGAGGTCACGCTTGAGGGCAGTGGGCTGACGCGTCCACGACTGCTTGTCGAGGCGGAAGAACACGTAGATCTTCATCAGCGCGCCGATGATCTGGTTGTAATAGAGAATCAGCG  >  NC_002947/1471257‑1471551

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