Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,874,820	A→C	5.6%	S451A (TCG→GCG) ‡	ilvA‑II ←	threonine deaminase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,874,820	0	A	C	5.6%	225.3 / 4.1	71	S451A (TCG→GCG) ‡	ilvA‑II	threonine deaminase
Reads supporting (aligned to +/- strand):  ref base A (32/35);  new base C (2/2);  total (34/37)
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.98e-01

GTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCAATTCGTTG  >  NC_002947/5874674‑5874968                                                                                                                                                   |                                                                                                                                                         gTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCgag                                                                                                                                                         <  4:335433/148‑1 (MQ=255) gTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCgag                                                                                                                                                         >  3:335433/1‑148 (MQ=255)    ggCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATg                                                                                                                                                      <  3:209405/148‑1 (MQ=255)     gCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGtt                                                                                                                                                    >  2:121471/1‑149 (MQ=255)      ccGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTc                                                                                                                                                   >  1:2841/1‑149 (MQ=255)           ttttCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCTCGATGTTCCacta                                                                                                                                              <  6:297601/149‑4 (MQ=255)               cGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑Gcgcc‑cg                                                                                                                                        >  4:222863/1‑149 (MQ=255)                 gTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑Gcgcc‑cgcc                                                                                                                                      >  1:155004/1‑149 (MQ=255)                  tctcGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑Gcgcc‑cgcc                                                                                                                                      >  6:123862/1‑148 (MQ=255)                    tcGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAg                                                                                                                                   >  3:170491/1‑149 (MQ=255)                        cccAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCCCGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTACA‑GCGCC‑CGCCCCGCttg                                                                                                                               >  5:55010/1‑149 (MQ=255)                          cAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCttgt‑‑t                                                                                                                           <  3:222863/149‑1 (MQ=255)                           aGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGTAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCttgt‑‑t                                                                                                                           >  7:331978/1‑148 (MQ=255)                              aCGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAg                                                                                                                        <  5:89443/148‑1 (MQ=255)                               cGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGa                                                                                                                      <  2:134490/149‑1 (MQ=255)                                gggTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAgg                                                                                                                       >  6:345392/1‑147 (MQ=255)                                gggTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGaa                                                                                                                     <  4:348789/149‑1 (MQ=255)                                   tAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGTAAGtt                                                                                                                  <  2:132387/149‑1 (MQ=255)                                     aCCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGa                                                                                                                <  5:345392/149‑1 (MQ=255)                                      ccGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCATCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGaa                                                                                                               <  8:132794/149‑1 (MQ=255)                                       cGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGaa                                                                                                               >  1:336847/1‑148 (MQ=255)                                       cGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGaa                                                                                                               >  1:331665/1‑148 (MQ=255)                                          tCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCCCCTGCAGCCCCGCCACCACGCGCCCGTAAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGCACAGt                                                                                                           >  4:268894/1‑149 (MQ=255)                                          tCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGt                                                                                                           >  6:363973/1‑149 (MQ=255)                                              ggCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGTTGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACa                                                                                                             <  1:188799/143‑1 (MQ=255)                                              ggCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAAc                                                                                                              <  6:197752/142‑1 (MQ=255)                                              ggCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACa                                                                                                             <  5:289694/143‑1 (MQ=255)                                               gCCAGCGCTGCCGGCACCAGGTGACGTGCGTCCTCCGGCACCTTCATCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGcccc                                                                                                      <  6:55010/149‑1 (MQ=255)                                                  aGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTgg                                                                                                   <  7:365834/149‑1 (MQ=255)                                                    cgcTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGcc                                                                                                 <  4:203417/149‑1 (MQ=255)                                                           ggCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCggg                                                                                          >  7:131343/1‑149 (MQ=255)                                                           ggCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCggg                                                                                          <  4:129789/149‑1 (MQ=255)                                                           ggCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCggg                                                                                          <  8:196735/149‑1 (MQ=255)                                                               ccAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAAt                                                                                      <  4:375888/149‑1 (MQ=255)                                                                cAGGTGACGTTCGTCCTCCGGCACCTGCAGTCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAAtt                                                                                     >  6:178857/1‑149 (MQ=255)                                                                 aGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTc                                                                                    >  8:315282/1‑149 (MQ=255)                                                                   gTGACGTTCTTCCTACGGCACCTGTAGCCCCGCCACCACGCGAGCGGCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCg                                                                                   <  8:331978/148‑1 (MQ=255)                                                                    tGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGaa                                                                                 <  6:66566/149‑1 (MQ=255)                                                                     gACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAg                                                                                >  5:133223/1‑149 (MQ=255)                                                                      aCGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCGATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAgc                                                                               <  6:237851/149‑1 (MQ=255)                                                                            gTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCgcagc                                                                          <  5:123862/148‑1 (MQ=255)                                                                               ctccGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCAcc                                                                       <  4:156651/148‑1 (MQ=255)                                                                               ctccGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCa                                                                      >  2:123282/1‑149 (MQ=255)                                                                                tccGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCAt                                                                     >  3:312735/1‑149 (MQ=255)                                                                                 ccGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCAtt                                                                    >  3:265116/1‑149 (MQ=255)                                                                                    gCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTtcg                                                                 >  6:284985/1‑149 (MQ=255)                                                                                     cACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTtcgt                                                                >  4:126072/1‑149 (MQ=255)                                                                                        cTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGc                                                             >  7:37504/1‑149 (MQ=255)                                                                                          gCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTg                                                           >  2:48238/1‑149 (MQ=255)                                                                                          gCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTg                                                           >  2:336664/1‑149 (MQ=255)                                                                                          gCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTg                                                           >  4:360549/1‑149 (MQ=255)                                                                                           cAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTg                                                           <  2:276182/148‑1 (MQ=255)                                                                                                ccGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCg                                                     >  5:45268/1‑149 (MQ=255)                                                                                                ccGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCg                                                     >  1:230906/1‑149 (MQ=255)                                                                                                         cgcgCCCGTGATCCGCGCGATGGTTTTGGTATTGGTAAATAGATATGTTACATGGGCCTCGCCCAGCTTTTTCTCTAGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGa                                                <  1:121471/149‑1 (MQ=255)                                                                                                          gcgcCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGgc                                           <  3:126072/149‑1 (MQ=255)                                                                                                           cgcCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGGCGCTGGCACCGGCCGAGTGgcc                                          <  4:312735/149‑1 (MQ=255)                                                                                                                 tCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGAc                                    <  4:265116/149‑1 (MQ=255)                                                                                                                   aGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGAcc                                   <  1:48238/148‑1 (MQ=255)                                                                                                                   aGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCa                                  <  5:284985/149‑1 (MQ=255)                                                                                                                     ccGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCata                                >  6:116016/1‑149 (MQ=255)                                                                                                                     ccGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCata                                >  2:326508/1‑149 (MQ=255)                                                                                                                                ttGCGGTTTTTGAACATCGATATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTTAACAGTGCCCCTTGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATg                     <  7:58727/149‑1 (MQ=255)                                                                                                                                ttGCGGTAGTGGAACATCGATATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATg                     <  5:281287/149‑1 (MQ=255)                                                                                                                                 cagattcgTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGc                    <  3:304902/142‑1 (MQ=255)                                                                                                                                  gCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGc                    <  1:326508/148‑1 (MQ=255)                                                                                                                                   cGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAg                  <  5:363973/149‑1 (MQ=255)                                                                                                                                     gTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCt                >  1:13631/1‑149 (MQ=255)                                                                                                                                      tAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGACCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTAGTCGCTGGAACCGGCAGAGTGgccgcc                                       >  5:202487/1‑125 (MQ=255)                                                                                                                                      tAGTGGAACATCGAGATGGTCCA‑TCGCC‑CGCCCATCTTGT‑‑TGAGGAAGTTTAACAGTTCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGTCCGAGTGtccgcc                                       <  6:202487/125‑1 (MQ=255)                                                                                                                                       gctGGAACATCTCGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTg              >  7:60501/3‑149 (MQ=255)                                                                                                                                         tGGAACATCTCGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATg                     <  5:250454/140‑1 (MQ=255)                                                                                                                                         tGGAACATCTCGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATg                     >  6:250454/1‑140 (MQ=255)                                                                                                                                                  gagaTGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCAATTCGtt   <  7:81468/149‑1 (MQ=255)                                                                                                                                                  gagaTGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTCGCCAATTCGt    <  2:230906/148‑1 (MQ=255)                                                                                                                                                  gagaTGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGGTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCAATTCGt    <  2:2841/148‑1 (MQ=255)                                                                                                                                                   agaTGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCAATTCGTTg  <  2:331665/149‑1 (MQ=255)                                                                                                                                                   agaTGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCAATTCGTTg  <  2:336847/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                         GTAGGCCGGGTTTTCGGTCTCGTCCCAGTACGGGTAACCGATCTTGGCCAGCGCTGCCGGCACCAGGTGACGTTCGTCCTCCGGCACCTGCAGCCCCGCCACCACGCGCCCGTCAGCCGCGCCATGGTTGCGGTAGTGGAACATCGAGATGTTCCA‑GCGCC‑CGCCCAGCTTGT‑‑TGAGGAAGTTGAACAGTGCCCCTGGCCGCTCGGGGAATTCGAAGCGCAGCACCATTTCGTCGCTGGCACCGGCCGAGTGGCCGCCGACCATATGACGAATATGCAGCTTGGCCAATTCGTTG  >  NC_002947/5874674‑5874968

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