Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,677,429	C→A	5.5%	A411E (GCG→GAG) ‡	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,677,429	0	C	A	5.5%	291.3 / 4.4	90	A411E (GCG→GAG) ‡	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base C (44/41);  new base A (3/2);  total (47/43)
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.28e-01

CGTCAGCGACGAGCAAACCTGCGAGACCATCGCAGAAGTGTTCGCCGCCACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGC  >  NC_002947/1677283‑1677577                                                                                                                                                   |                                                                                                                                                     cGTCAGCGACGAGCAAACCTGCGAGACCATCGCAGAAGTGTTCGCCGCCACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACcgcg                                                                                                                                                     >  3:347131/1‑148 (MQ=255)     aGCGACGAGCAAACCTGCGAGACCATCGCAGAAGTGTTCGCCGCCACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCAccc                                                                                                                                                <  4:347131/149‑1 (MQ=255)           gaGCAAACCTGCGAGACCATCGCAGAAGTGTTCGCCGCCACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCAGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTc                                                                                                                                            >  3:158553/1‑147 (MQ=255)               aaaCCTGCGAGACCATCGCAGAAGTGTTCGCCGCCACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTc                                                                                                                                      <  3:28456/149‑1 (MQ=255)                 aCCTGCGAGACCATCGCAGAAGTGTTCGCCGCCACCTGTGAAGTGCTCGACCCGCATACCGGAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTccc                                                                                                                                    <  2:242532/149‑1 (MQ=255)                 aCCTGCGAGACCATCGCAGAAGTGTTCGCCGCCACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTccc                                                                                                                                    >  4:99929/1‑149 (MQ=255)                 aCCTGCGAGACCATCGCAGAAGTGTTCGCCGCCACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTccc                                                                                                                                    <  8:193467/149‑1 (MQ=255)                    tGCGAGACCATCGCAGAAGTGTTCGCCGCCACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCGCGGTACCCAGCAccac                                                                                                                                               <  1:200239/135‑3 (MQ=255)                    tGCGAGACCATCGCAGAAGTGTTCGCCGCCACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCGCGGTACCCAGCAccac                                                                                                                                               >  2:200239/1‑133 (MQ=255)                     gCGAGACCATCGCAGAAGTGTTCGCCGCCACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCgg                                                                                                                                  >  2:317645/1‑147 (MQ=255)                      cGAGACCATCGCAGAAGTGTTCGCCGCCACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGaa                                                                                                                                <  3:155004/148‑1 (MQ=255)                         gaCCATCGCAGAAGTGTTCGCCGCCACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCgg                                                                                                                            >  5:148645/1‑149 (MQ=255)                           ccATCGCAGAAGTGTTCGCCGCCACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTg                                                                                                                          <  1:203680/149‑1 (MQ=255)                            cATCGCAGAAGTGTTCGCCGCCACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTgg                                                                                                                         <  1:181076/149‑1 (MQ=255)                            cATCGCAGAAGTGTTCGCCGCCACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTgg                                                                                                                         <  3:288467/149‑1 (MQ=255)                              tCGCAGAAGTGTTCGCCGCCACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGgag                                                                                                                       <  1:221387/149‑1 (MQ=255)                                gCAGAAGTGTTCGCCGCCACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGaa                                                                                                                     <  5:142588/149‑1 (MQ=255)                                    aaGTGTTCGCCGCCACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGcc                                                                                                                 >  2:236468/1‑149 (MQ=255)                                    aaGTGTTCGCCGCCACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCAGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCCGTCAAGTTCCCGGAAGCGGTGGAGAAGGcc                                                                                                                 >  5:128285/1‑149 (MQ=255)                                       tgtTCGCCGCCACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGgt                                                                                                              <  3:99929/149‑1 (MQ=255)                                         ttagccgccACCGGTGAAGTGCTCGACCCCCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTCCCGCCGCAGTCTGGACACGCCAAAGGTGGTGCTGGGTACCCCGCACCCGGTCAAGTTCCCCGAAGCGGTCGAGAAGGCCgcagt                                                                                                            >  6:242794/4‑145 (MQ=255)                                            gccgccACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTggg                                                                                                         <  7:280170/149‑1 (MQ=255)                                            gccgccACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTggg                                                                                                         >  8:280170/1‑149 (MQ=255)                                             ccgccACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGc                                                                                                        <  7:384301/149‑1 (MQ=255)                                               gccACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGc                                                                                                        >  2:89396/1‑147 (MQ=255)                                               gccACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGc                                                                                                        >  8:233302/1‑147 (MQ=255)                                               gccACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCaa                                                                                                      <  2:279899/149‑1 (MQ=255)                                               gccACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCaa                                                                                                      <  7:150107/149‑1 (MQ=255)                                                 cACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGc                                                                                                                              <  7:133258/123‑1 (MQ=255)                                                 cACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGc                                                                                                                              >  8:133258/1‑123 (MQ=255)                                                      gTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTg                                                                                               <  2:116398/149‑1 (MQ=255)                                                      gTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCAGGTGTGGGCAAGGCGCTg                                                                                               >  1:372388/1‑149 (MQ=255)                                                           gTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCt                                                                                          <  2:232879/149‑1 (MQ=255)                                                            tGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTa                                                                                         >  2:186976/1‑149 (MQ=255)                                                             gCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTAc                                                                                        >  5:320737/1‑149 (MQ=255)                                                              cTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCt                                                                                                >  4:299245/1‑140 (MQ=255)                                                              cTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCt                                                                                                <  3:299245/140‑1 (MQ=255)                                                              cTCGACCCGCATACCGCAATAGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACAGGCCAATGGTGGTGCTGGGTACAGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCACGGCGCTGGAGCTAcc                                                                                       >  4:162163/1‑149 (MQ=255)                                                               tCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGa                                                                                             <  6:76126/142‑1 (MQ=255)                                                               tCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGa                                                                                             >  5:76126/1‑142 (MQ=255)                                                                 gACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCcgc                                                                                    >  8:130585/1‑149 (MQ=255)                                                                  aCCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCcgcg                                                                                   >  6:190420/1‑149 (MQ=255)                                                                           ccGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGc                                                                          <  5:190420/149‑1 (MQ=255)                                                                            cGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCGCGGTACCCAGCAc                                                                                                                                                  <  4:240961/76‑1 (MQ=255)                                                                            cGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCGCGGTACCCAGCAc                                                                                                                                                  >  3:240961/1‑76 (MQ=255)                                                                            cGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCAGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCACGGAAGCGGTGGTGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGc                                                                          >  8:377309/1‑148 (MQ=255)                                                                             gCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGa                                                                        >  8:101399/1‑149 (MQ=255)                                                                              cAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGAc                                                                       <  7:130585/149‑1 (MQ=255)                                                                                 tCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGAcc                                                                      >  7:31544/1‑147 (MQ=255)                                                                                 tCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGAcc                                                                      <  8:31544/147‑1 (MQ=255)                                                                                     cGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCa                                                                <  1:236468/149‑1 (MQ=255)                                                                                      gTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGAcc                                                                      <  3:284738/142‑1 (MQ=255)                                                                                      gTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGAcc                                                                      >  4:284738/1‑142 (MQ=255)                                                                                        cAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGcc                                                             >  2:323222/1‑149 (MQ=255)                                                                                         aaGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGGCATTGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGTAGCGGTGGAGAAGGCCGGGGGGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGcc                                                             <  5:242794/148‑1 (MQ=255)                                                                                         aaGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAg                                                               >  2:232829/1‑146 (MQ=255)                                                                                         aaGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAg                                                               <  1:232829/146‑1 (MQ=255)                                                                                            gCTGCCCGCGTGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTg                                                          >  4:398482/1‑148 (MQ=255)                                                                                            gCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGtt                                                                  <  5:3639/140‑1 (MQ=255)                                                                                            gCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGtt                                                                  >  6:3639/1‑140 (MQ=255)                                                                                                cccGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAaga                                                      <  7:101399/148‑1 (MQ=255)                                                                                                 ccGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGc                                                    >  2:225479/1‑149 (MQ=255)                                                                                                 ccGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGc                                                    >  2:391360/1‑149 (MQ=255)                                                                                                    cgAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGtt                                                 <  1:317645/149‑1 (MQ=255)                                                                                                     gAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTg                                                >  6:172672/1‑149 (MQ=255)                                                                                                        tGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCAc                                             >  4:227199/1‑149 (MQ=255)                                                                                                          ccgccgCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACg                                            <  1:323222/148‑1 (MQ=255)                                                                                                                cAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACGCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTgg                                     <  4:158553/149‑1 (MQ=255)                                                                                                                cAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTgg                                     >  3:17575/1‑149 (MQ=255)                                                                                                                       gACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGa                              >  8:210793/1‑149 (MQ=255)                                                                                                                        acacGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCTGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGAc                             <  4:17575/149‑1 (MQ=255)                                                                                                                            gCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGc                                               >  5:367163/1‑127 (MQ=255)                                                                                                                            gCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGc                                               <  6:367163/127‑1 (MQ=255)                                                                                                                            gCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCAGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTg                          >  1:116904/1‑148 (MQ=255)                                                                                                                             ccAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGTGAATGCCGGTGTGGGGAAGGCGGTTGTGCTACCCGGGCACCTAAGCGTCCTTTTCAGCCGTGAAGATCGTTGCACGGTCCTGGCCAACGACCTGaa                        <  6:148645/149‑1 (MQ=255)                                                                                                                             ccAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTg                          >  4:60962/1‑147 (MQ=255)                                                                                                                              ttgaccgggtgCGCGGTACCCAGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAg                       >  4:262986/8‑149 (MQ=255)                                                                                                                              cAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCaa                                 <  4:14152/139‑1 (MQ=255)                                                                                                                              cAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCaa                                 >  3:14152/1‑139 (MQ=255)                                                                                                                              cAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAg                       >  7:144618/1‑149 (MQ=255)                                                                                                                               aaTGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAgg                      <  5:172672/149‑1 (MQ=255)                                                                                                                                 tggtggTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCt                    >  3:29747/1‑149 (MQ=255)                                                                                                                                   gtggtgCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTg                   <  3:398482/148‑1 (MQ=255)                                                                                                                                       tgCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAg              <  7:323962/149‑1 (MQ=255)                                                                                                                                       tgCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAg              >  6:177985/1‑149 (MQ=255)                                                                                                                                         cTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAggg            >  2:85685/1‑149 (MQ=255)                                                                                                                                           gggTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCt          >  8:33440/1‑149 (MQ=255)                                                                                                                                              tACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTg       <  7:136116/149‑1 (MQ=255)                                                                                                                                              tACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTg       >  6:89164/1‑149 (MQ=255)                                                                                                                                               aCCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGa                                                         <  6:384421/98‑1 (MQ=255)                                                                                                                                               aCCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGa                                                         >  5:384421/1‑98 (MQ=255)                                                                                                                                                 cgcgCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTca    <  2:372388/149‑1 (MQ=255)                                                                                                                                                   cgcACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTcagc  <  6:320737/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                     CGTCAGCGACGAGCAAACCTGCGAGACCATCGCAGAAGTGTTCGCCGCCACCGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCTTTGTCAGC  >  NC_002947/1677283‑1677577

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