Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,677,420	T→G	5.6%	L408R (CTG→CGG) ‡	thrC →	threonine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,677,420	0	T	G	5.6%	284.5 / 7.1	90	L408R (CTG→CGG) ‡	thrC	threonine synthase
Reads supporting (aligned to +/- strand):  ref base T (46/39);  new base G (3/2);  total (49/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 = 5.01e-01

GACTCGCTGGCCGTCAGCGACGAGCAAACCTGCGAGACCATCGCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGG  >  NC_002947/1677272‑1677567                                                                                                                                                      |                                                                                                                                                    gACTCGCTGGCCGTCAGCGACGAGCAAACCTGCGAGACCATCGCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCt                                                                                                                                                     <  8:205951/149‑1 (MQ=255)    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                                                                                                                                     <  1:200239/135‑3 (MQ=255)                               tGCGAGACCATCGCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCGCGGTACCCAGCAccac                                                                                                                                     >  2:200239/1‑133 (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                                                                                                               <  1:181076/149‑1 (MQ=255)                                       cATCGCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTgg                                                                                                               <  3:288467/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                                                                                               >  8:280170/1‑149 (MQ=255)                                                       gccgccAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTggg                                                                                               <  7:280170/149‑1 (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                                                                                                                    >  8:133258/1‑123 (MQ=255)                                                            cAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGc                                                                                                                    <  7:133258/123‑1 (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                                                            <  8:31544/147‑1 (MQ=255)                                                                                             tCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGAcc                                                            >  7:31544/1‑147 (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)                                                                                                       ggCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCt                                                                                                                                                     <  4:240152/48‑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                                     <  6:367163/127‑1 (MQ=255)                                                                                                                                        gCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGc                                     >  5:367163/1‑127 (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)                                                                                                                                                      |                                                                                                                                                    GACTCGCTGGCCGTCAGCGACGAGCAAACCTGCGAGACCATCGCAGAAGTGTTCGCCGCCAC‑CGGTGAAGTGCTCGACCCGCATACCGCAATCGGCGTCAAGGCTGCCCGCGAGTGCCGCCGCAGTCTGGACACGCCAATGGTGGTGCTGGGTACCGCGCACCCGGTCAAGTTCCCGGAAGCGGTGGAGAAGGCCGGTGTGGGCAAGGCGCTGGAGCTACCCGCGCACCTCAGCGACCTGTTCAGCCGTGAAGAGCGTTGCACGGTCCTGGCCAACGACCTGAAGGCTGTACAGGG  >  NC_002947/1677272‑1677567

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