Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,783,115	T→C	10.1%	E1400G (GAG→GGG)	pvdJ ←	non‑ribosomal peptide synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,783,115	0	T	C	10.1%	181.8 / 9.1	69	E1400G (GAG→GGG)	pvdJ	non‑ribosomal peptide synthetase
Reads supporting (aligned to +/- strand):  ref base T (29/33);  new base C (3/4);  total (32/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.71e-01

CCTGCTCGCTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCAGGTCGAGGTCGTCGAACAGGTCCTCGCGGCCGTGG  >  NC_002947/4782975‑4783252                                                                                                                                             |                                                                                                                                          ccTGCTCGCTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGAGGTCCCGGCCa                                                                                                                                     >  1:142367/1‑147 (MQ=255) ccTGCTCGCTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGAGGTCCCGGCCa                                                                                                                                     >  7:241428/1‑147 (MQ=255)  cTGCTCGCTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGAGGTCCCGGCCaaag                                                                                                                                  <  4:80979/149‑3 (MQ=255)  cTGCTCGCTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGAGGTCCCGGCCaaag                                                                                                                                  <  8:241428/149‑3 (MQ=255)  cTGCTCGCTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGAGGTCCCGGCCaaag                                                                                                                                  <  1:118524/149‑3 (MQ=255)   tGCTCGCTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAActg                                                                                                                                 >  1:76230/1‑149 (MQ=255)     cTCGCTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAActgct                                                                                                                               <  1:194087/149‑1 (MQ=255)      tCGCTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGAGGTCCCGGCCaaa                                                                                                                                   >  6:140702/1‑143 (MQ=255)      tCGCTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGAGGTCCCGGCCaaa                                                                                                                                   <  5:140702/144‑2 (MQ=255)         cTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCgg                                                                                                                                        <  2:266956/136‑1 (MQ=255)         cTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCtt                                                                                                                           <  6:41248/149‑1 (MQ=255)         cTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCAGGGACCTCgg                                                                                                                                        >  1:266956/1‑136 (MQ=255)            aaGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGAt                                                                                                                        <  1:212354/149‑1 (MQ=255)             aGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCg                                                                                                                                         <  5:103037/131‑1 (MQ=255)             aGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCg                                                                                                                                         >  6:103037/1‑131 (MQ=255)                 gCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACtgctgc                                                                                                                             >  4:345982/1‑139 (MQ=255)                 gCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACtgctgc                                                                                                                             <  3:345982/139‑1 (MQ=255)                    aTCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCgg                                                                                                                                        >  7:208948/1‑125 (MQ=255)                    aTCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCgg                                                                                                                                        <  8:208948/125‑1 (MQ=255)                        aaCTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCt                                                                                                            >  2:345981/1‑149 (MQ=255)                          cTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTgg                                                                                                          >  5:145964/1‑149 (MQ=255)                          cTGGCCGAGGTACGTGAAGGGCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTACTGCTTGATGTGTTTGATGCTgg                                                                                                          <  7:140861/149‑1 (MQ=255)                             gCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGgccg                                                                                                       >  2:65979/1‑149 (MQ=255)                                         gAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTgg                                                                                             >  5:37179/1‑147 (MQ=255)                                         gAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTgg                                                                                             >  7:210147/1‑147 (MQ=255)                                         gAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTgg                                                                                             >  8:105518/1‑147 (MQ=255)                                          aaGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAg                                                                                          >  6:187203/1‑149 (MQ=255)                                           aGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGc                                                                                         <  6:29809/149‑1 (MQ=255)                                              tCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGgt                                                                                      <  5:341558/149‑1 (MQ=255)                                              tCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGgt                                                                                      <  5:262369/149‑1 (MQ=255)                                               cACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGAt                                                                                                                        >  6:192213/1‑114 (MQ=255)                                               cACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGAt                                                                                                                        <  5:192213/114‑1 (MQ=255)                                                aCCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGgtgt                                                                                    <  6:37179/149‑1 (MQ=255)                                                  ccGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGc                                                                                                 >  8:305359/1‑134 (MQ=255)                                                  ccGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGc                                                                                                 <  7:305359/134‑1 (MQ=255)                                                     cgGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGtt                                                                               <  2:76230/149‑1 (MQ=255)                                                        cTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCa                                                                            >  5:77301/1‑149 (MQ=255)                                                            ggCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCgg                                                                        <  5:37213/149‑1 (MQ=255)                                                                  gCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCaa                                                                                                                                    >  5:178743/1‑83 (MQ=255)                                                                  gCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCaa                                                                                                                                    <  6:178743/83‑1 (MQ=255)                                                                  gCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGaa                                                                     >  4:194323/1‑146 (MQ=255)                                                                      ccAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCAcc                                                                          <  1:276291/137‑1 (MQ=255)                                                                      ccAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCAcc                                                                          >  2:276291/1‑137 (MQ=255)                                                                      ccAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTg                                                              >  5:42899/1‑149 (MQ=255)                                                                      ccAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTg                                                              >  3:97154/1‑149 (MQ=255)                                                                      ccAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTg                                                              <  3:194323/149‑1 (MQ=255)                                                                       cAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTgg                                                             <  1:217424/149‑1 (MQ=255)                                                                                  cGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACtgctgc                                                                                                                             >  7:79984/1‑74 (MQ=255)                                                                                  cGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACtgctgc                                                                                                                             <  8:79984/74‑1 (MQ=255)                                                                                      cgcgGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACgg                                              <  8:210147/149‑1 (MQ=255)                                                                                       gcgGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGAt                                                                                                                        >  2:355164/1‑74 (MQ=255)                                                                                       gcgGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGAt                                                                                                                        <  1:355164/74‑1 (MQ=255)                                                                                         ggCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGc                                           <  1:345981/149‑1 (MQ=255)                                                                                               gCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCt                                       >  5:214918/1‑147 (MQ=255)                                                                                               gCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCa                                     <  6:77301/149‑1 (MQ=255)                                                                                               gCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCa                                     <  1:65979/149‑1 (MQ=255)                                                                                                    ggTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCCCCGGGAACAGAGTGGTGAACCAGCCAACGGTGCGGCGCgggtcg                                >  4:330408/1‑149 (MQ=255)                                                                                                        aCGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCaggtcgagg                             <  5:187203/148‑1 (MQ=255)                                                                                                            agcaCGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCAGGTCGAGgtcgtc                        <  6:214918/149‑1 (MQ=255)                                                                                                               aCGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCAGGTCGAGGTCGTCGaa                     <  1:338595/149‑1 (MQ=255)                                                                                                                 gCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCAGGTCGAGGTCGTCGaa                     >  8:5241/1‑147 (MQ=255)                                                                                                                 gCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCAGGTCGAGGTCGTCGAACa                   <  4:97154/149‑1 (MQ=255)                                                                                                                           aTACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCAGGTCGAGGTCGTCGAACAGGTCCTcgcg         <  6:42899/149‑1 (MQ=255)                                                                                                                             aCCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCAGGTCGAGGTCGTCGAACAGGTCCTCGCGGc       >  3:86826/1‑149 (MQ=255)                                                                                                                               cTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCAACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCAGGTCGAGGTCGTCGAACAGGTCCTCGCGGCCg     >  4:129491/1‑149 (MQ=255)                                                                                                                                tttGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCaggtcgag                              <  4:285736/123‑1 (MQ=255)                                                                                                                                  tGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCAGGTCGAGGTCGTCGAACAGGTCCTCGCGGCCGTgg  >  6:210715/1‑149 (MQ=255)                                                                                                                                  tGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCAGGTCGAGGTCGTCGAACAGGTCCTCGCGGCCGTgg  >  6:215840/1‑149 (MQ=255)                                                                                                                                   ggCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCAGGTCGAGGTCGTCGAACAGGTCCTCGCGGCCGTgg  <  7:105518/148‑1 (MQ=255)                                                                                                                                             |                                                                                                                                          CCTGCTCGCTGAAGCTGCCATCGAACTGGCCGAGGTAGTTGAAGGTCACCCGCGGCTGCGGCAACGCAGCCAGGCGCTGGGCGAACGCGGCTTCGCCCAGGTAACGCAGCACGCCGTAACCGATACCTTTGGCCGGGACCTCGGCCAACTGCTGCTTGATGTGTTTGATGCTGGCCGCCACGCTTGGCAGCGGTGTCAGTTGCACCGGGAACAGACTGGTGAACCAGCCAACGGTGCGGCTCAGGTCGAGGTCGTCGAACAGGTCCTCGCGGCCGTGG  >  NC_002947/4782975‑4783252

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