Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,677,421	G→C	5.7%	L408L (CTG→CTC) ‡	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,677,421	0	G	C	5.7%	292.3 / 5.9	88	L408L (CTG→CTC) ‡	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base G (44/39);  new base C (3/2);  total (47/41)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.96e-01

TCGCTGGCCGTCAGCGACGAGCAAACCTGCGAGACCATCGCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCT  >  NC_002947/1677275‑1677569                                                                                                                                                    |                                                                                                                                                     tCGCTGGCCGTCAGCGACGAGCAAACCTGCGAGACCATCGCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTggg                                                                                                                                                    >  8:384301/1‑149 (MQ=255)  cGCTGGCCGTCAGCGACGAGCAAACCTGCGAGACCATCGCAGAAGTGTTCGCCGCCACTAGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTggg                                                                                                                                                    >  1:116398/1‑149 (MQ=255)     tGGCCGTCAGCGACGAGCAAACCTGCGAGACCATCGCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTAcc                                                                                                                                                >  2:221387/1‑149 (MQ=255)     tGGCCGTCAGCGACGAGCAAACCTGCGAGACCATCGCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGCGTAcc                                                                                                                                                >  7:193467/1‑149 (MQ=255)         cGTCAGCGACGAGCAAACCTGCGAGACCATCGCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACcgcg                                                                                                                                             >  3:347131/1‑148 (MQ=255)             aGCGACGAGCAAACCTGCGAGACCATCGCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCAccc                                                                                                                                        <  4:347131/149‑1 (MQ=255)                   gaGCAAACCTGCGAGACCATCGCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCAGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTc                                                                                                                                    >  3:158553/1‑147 (MQ=255)                       aaaCCTGCGAGACCATCGCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTc                                                                                                                              <  3:28456/149‑1 (MQ=255)                         aCCTGCGAGACCATCGCAGAAGTGTTCGCCGCCAC‑CTGTGAAGTGCTCGACCCGCATACCGGAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTccc                                                                                                                            <  2:242532/149‑1 (MQ=255)                         aCCTGCGAGACCATCGCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTccc                                                                                                                            <  8:193467/149‑1 (MQ=255)                         aCCTGCGAGACCATCGCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTccc                                                                                                                            >  4:99929/1‑149 (MQ=255)                            tGCGAGACCATCGCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCGCGGTACCCAGCAccac                                                                                                                                       >  2:200239/1‑133 (MQ=255)                            tGCGAGACCATCGCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCGCGGTACCCAGCAccac                                                                                                                                       <  1:200239/135‑3 (MQ=255)                             gCGAGACCATCGCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCgg                                                                                                                          >  2:317645/1‑147 (MQ=255)                              cGAGACCATCGCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGaa                                                                                                                        <  3:155004/148‑1 (MQ=255)                                 gaCCATCGCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCgg                                                                                                                    >  5:148645/1‑149 (MQ=255)                                   ccATCGCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTg                                                                                                                  <  1:203680/149‑1 (MQ=255)                                    cATCGCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTgg                                                                                                                 <  3:288467/149‑1 (MQ=255)                                    cATCGCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTgg                                                                                                                 <  1:181076/149‑1 (MQ=255)                                      tCGCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGgag                                                                                                               <  1:221387/149‑1 (MQ=255)                                        gCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGaa                                                                                                             <  5:142588/149‑1 (MQ=255)                                            aaGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGcc                                                                                                         >  2:236468/1‑149 (MQ=255)                                            aaGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCAGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCCGTCAAGTTCCCGGAAGCGGTGGAGAAGGcc                                                                                                         >  5:128285/1‑149 (MQ=255)                                               tgtTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGgt                                                                                                      <  3:99929/149‑1 (MQ=255)                                                 ttagccgccAC‑CGGTGAAGTGCTCGACCCCCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTCCCGCCGCAGTCTGGACACGCCAAAGGTGGTGCTGGGTACCCCGCACCCGGTCAAGTTCCCCGAAGCGGTCGAGAAGGCCgcagt                                                                                                    >  6:242794/4‑145 (MQ=255)                                                    gccgccAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTggg                                                                                                 <  7:280170/149‑1 (MQ=255)                                                    gccgccAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTggg                                                                                                 >  8:280170/1‑149 (MQ=255)                                                     ccgccAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGc                                                                                                <  7:384301/149‑1 (MQ=255)                                                       gccAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGc                                                                                                >  8:233302/1‑147 (MQ=255)                                                       gccAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGc                                                                                                >  2:89396/1‑147 (MQ=255)                                                       gccAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCaa                                                                                              <  7:150107/149‑1 (MQ=255)                                                       gccAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCaa                                                                                              <  2:279899/149‑1 (MQ=255)                                                         cAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGc                                                                                                                      <  7:133258/123‑1 (MQ=255)                                                         cAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGc                                                                                                                      >  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                                                                                        <  3:299245/140‑1 (MQ=255)                                                                       cTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCt                                                                                        >  4:299245/1‑140 (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                                                                                                                                          >  3:240961/1‑76 (MQ=255)                                                                                     cGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCGCGGTACCCAGCAc                                                                                                                                          <  4:240961/76‑1 (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                                                              >  4:284738/1‑142 (MQ=255)                                                                                               gTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGAcc                                                              <  3:284738/142‑1 (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:391360/1‑149 (MQ=255)                                                                                                          ccGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGc                                            >  2:225479/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)                                                                                                                                                    |                                                                                                                                                     TCGCTGGCCGTCAGCGACGAGCAAACCTGCGAGACCATCGCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGGCT  >  NC_002947/1677275‑1677569

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