Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,359,184	A→G	5.8%	L258L (TTG→CTG)	PP_3834 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,359,184	0	A	G	5.8%	254.2 / 2.9	84	L258L (TTG→CTG)	PP_3834	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base A (41/38);  new base G (3/2);  total (44/40)
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.80e-01

GCCTTGCTGCGCAGGCGTTCGGGTGTTTCGTAAAAACCGAAATGCAACCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAGACCTGCATGCCGGACTCATCGATAAGGCATTCGTGAC  >  NC_002947/4359036‑4359330                                                                                                                                                     |                                                                                                                                                     gCCTTGCTGCGCAGGCGTTCGGGTGTTTCGTAAAAACCGAAATGCAACCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCaa                                                                                                                                                      <  2:294263/149‑1 (MQ=255) gCCTTGCTGCGCAGGCGTTCGGGTGTTTCGTAAAAACCGAAATGCAACCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCaa                                                                                                                                                      <  8:46923/149‑1 (MQ=255) gCCTTGCTGCGCAGGCGTTCGGGTGTTTCGTAAAAACCGAAATGCAACCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCaa                                                                                                                                                      <  3:127719/149‑1 (MQ=255)        tgcGCAGGCGTTCGGGTGTTTCGTAAAAACCGAAATGCAACCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCg                                                                                                                                              >  3:8293/1‑149 (MQ=255)          cgcAGGCGTTCGGGTGTTTCGTAAAAACCGAAATGCAACCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATTCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCTTG‑GCCACCgg                                                                                                                                             <  3:321112/148‑2 (MQ=255)          cgcAGGCGTTCGGGTGTTTCGTAAAAACCGAAATGCAACCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCTTG‑GCCACCgg                                                                                                                                             <  6:119022/148‑2 (MQ=255)            cAGGCGTTCGGGTGTTTCGTAAAAACCGAAATGCAACCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCaa                                                                                                                                          <  1:153901/149‑1 (MQ=255)            cAGGCGTTCGGGTGTTTCGTAAAAACCGAAATGCAACCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCaa                                                                                                                                          <  1:60265/149‑1 (MQ=255)               gCGTTCGGGTGTTTCGTAAAAACCGAAATGCAACCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCa                                                                                                                                       >  6:135685/1‑149 (MQ=255)                  ttCGGGTGTTTCGTAAAAACCGAAATGCAACCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGTCCCGGTGGCCAA‑GCCACCGCCAAGCAttt                                                                                                                                    <  8:1382/149‑1 (MQ=255)                  ttCGGGTGTTTCGTAAAAACCGAAATGCAACCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCAtt                                                                                                                                     >  7:172330/1‑148 (MQ=255)                    cGGGTGTTTCGTAAAAACCGAAATGCAACCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCa                                                                                                                                  <  4:174338/149‑1 (MQ=255)                     gggTGTTTCGTAAAAACCGAAATGCAACCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCaa                                                                                                                                                      <  1:223295/129‑1 (MQ=255)                     gggTGTTTCGTAAAAACCGAAATGCAACCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCaa                                                                                                                                                      >  2:223295/1‑129 (MQ=255)                          tttCGTAAAAACCGAAATGCAACCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAg                                                                                                                            <  3:274818/149‑1 (MQ=255)                          tttCGTAAAAACCGAAATGCAACCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAg                                                                                                                            <  5:27022/149‑1 (MQ=255)                            tCGTAAAAACCGAAATGCAACCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAgg                                                                                                                           <  3:99836/148‑1 (MQ=255)                             cGTAAAAACCGAAATGCAACCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGtt                                                                                                                         >  1:250364/1‑149 (MQ=255)                               tAAAAACCGAAATGCAACCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTgg                                                                                                                       <  5:206576/149‑1 (MQ=255)                                 aaaaCCGAAATGCAACCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGt                                                                                                                      <  6:279240/148‑1 (MQ=255)                                  aaaCCGAAATGCAACCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTAc                                                                                                                    >  1:148598/1‑149 (MQ=255)                                    aCCGAAATGCAACCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTAcgc                                                                                                                  >  2:91952/1‑149 (MQ=255)                                       gAAATGCAACCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTAcgccgc                                                                                                               >  2:225362/1‑149 (MQ=255)                                              agccccGGTGCCTGCGGGGTGCTGGCACGCAGAAGGGCCAGGCCCAATGATGTCTTGGCGATCCGCGAAGGCCCCAGCACCAGGCTGACTGACGCGGTGGCCAG‑GCCACCGGCAGGCCTGGGATCGAAGGGTGGTAGGCCGctgctgct                                                                                                        >  8:134294/3‑147 (MQ=255)                                              aaCCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGa                                                                                                         >  8:134814/1‑148 (MQ=255)                                              aaCCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGAt                                                                                                        <  4:8293/149‑1 (MQ=255)                                               aCCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATa                                                                                                       >  3:234322/1‑149 (MQ=255)                                                  ccGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATAcgc                                                                                                    <  5:135685/149‑1 (MQ=255)                                                    ggTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACgcgc                                                                                                  <  7:183320/149‑1 (MQ=255)                                                    ggTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACgcgc                                                                                                  >  8:249022/1‑149 (MQ=255)                                                        ccTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCaa                                                                                               <  8:172330/148‑1 (MQ=255)                                                        ccTGCGGGGTGCTGGCACCCAGAAAGGCAAGGCCCAATGATGTCTTGCCAATACACGAAGGCCCCATAAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCAGCGAAGGTTGGTACGCCGCAGGTGATAACCGCCAAg                                                                                              >  7:234988/1‑149 (MQ=255)                                                          tGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGt                                                                                            >  5:27743/1‑149 (MQ=255)                                                               ggTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCAc                                                                                       >  3:284526/1‑149 (MQ=255)                                                                   cTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTAcc                                                                                   <  7:134814/149‑1 (MQ=255)                                                                   cTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTAcc                                                                                   >  7:139079/1‑149 (MQ=255)                                                                   cTGGCACCCAGAAAGGCCAGGCCCAATGATGGCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTAcc                                                                                   <  4:234322/149‑1 (MQ=255)                                                                    tGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTAcc                                                                                   >  1:43390/1‑148 (MQ=255)                                                                     ggCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTAccc                                                                                  >  7:56151/1‑148 (MQ=255)                                                                      gCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCaa                                                                                >  4:293464/1‑149 (MQ=255)                                                                          ccAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTgg                                                                            >  2:247693/1‑149 (MQ=255)                                                                          ccAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGATGCCTGACCCGGTGGCCAA‑GCCCCCGCCAAGCAGTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCAAATTgg                                                                            >  7:261737/1‑149 (MQ=255)                                                                            aGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCt                                                                          >  6:120726/1‑149 (MQ=255)                                                                            aGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCt                                                                          >  5:139615/1‑149 (MQ=255)                                                                             gAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCt                                                                          <  6:319226/148‑1 (MQ=255)                                                                             gAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCt                                                                          >  5:319226/1‑148 (MQ=255)                                                                               aaGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTg‑g                                                                       <  6:27743/148‑1 (MQ=255)                                                                                       gccCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTg                                                               >  7:75578/1‑148 (MQ=255)                                                                                        cccAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGaa                                                             <  3:154692/149‑1 (MQ=255)                                                                                          cAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGc                                                                 >  2:199268/1‑143 (MQ=255)                                                                                          cAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGc                                                                 <  1:199268/143‑1 (MQ=255)                                                                                           aaTGATGTCTTGCCAATACCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGACGCTGCTGATACGCGCCAAGGTGCGACTACCCAATTGGCTC‑GGGTGGCTGaacaa                                                          >  4:240543/1‑149 (MQ=255)                                                                                            atgatgTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGaacaac                                                         <  1:225362/149‑1 (MQ=255)                                                                                            atgatgTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGaacaac                                                         <  7:84304/149‑1 (MQ=255)                                                                                            atgatgTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGaacaac                                                         <  6:138883/149‑1 (MQ=255)                                                                                                      gCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGc                                                       <  2:351093/141‑1 (MQ=255)                                                                                                      gCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGc                                               >  5:281965/1‑149 (MQ=255)                                                                                                      gCCAATCCCCAAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCAGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCGGAACAACGcaagatcg                                                >  1:351093/1‑141 (MQ=255)                                                                                                          aTCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAgg                                           <  4:284526/149‑1 (MQ=255)                                                                                                             cccGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATa                                        <  2:43390/149‑1 (MQ=255)                                                                                                             cccGAAGGCCCCATCAACAGGCTGACTGACCCGGGGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATa                                        <  8:234988/149‑1 (MQ=255)                                                                                                              ccGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAg                                       <  8:139079/149‑1 (MQ=255)                                                                                                                 aaGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTGTGGGTGGCTGAACAACGCTTCCAGGCGAGGATAGAc                                     <  1:56237/149‑1 (MQ=255)                                                                                                                      cccATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAGACCTGCAt                               >  3:115433/1‑149 (MQ=255)                                                                                                                            aaCAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAGACCTGCATGCCGGa                         >  6:137955/1‑149 (MQ=255)                                                                                                                               aGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAGGCCTGCATGCCGGACTc                      <  8:75578/149‑1 (MQ=255)                                                                                                                                 gCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAGACCTGCATGCCGGACTCAt                    >  6:256923/1‑149 (MQ=255)                                                                                                                                 gCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAGACCTGCATGCCGGACTCAt                    >  4:150722/1‑149 (MQ=255)                                                                                                                                 gCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAGACCTGCATGCCGGACTCAt                    >  4:293334/1‑149 (MQ=255)                                                                                                                                  ctgactgaCCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAGACCTGCATGCCGGACTCATc                   <  2:250364/149‑1 (MQ=255)                                                                                                                                     actgacCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAGACCTGCATGCCGGACTCATCg                  <  1:248330/147‑1 (MQ=255)                                                                                                                                      ctgacCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAGACCTGCATGCCg                           <  1:131335/137‑1 (MQ=255)                                                                                                                                      ctgacCCGGTGGCAAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAGACCTGCATGCCg                           >  2:131335/1‑137 (MQ=255)                                                                                                                                       tgacCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAGACCTGCATGCCGGACTCATCGAt                >  5:282794/1‑147 (MQ=255)                                                                                                                                        ttggcGGTGGCTTG‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGCTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGCGGATAGACCTGCATGCCGGACTCAGCGATatg             >  8:306460/5‑147 (MQ=255)                                                                                                                                         tggcGGTGGCTTG‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAGACCTGCATGCCGGACTCATCGATAAgg            >  1:169615/4‑149 (MQ=255)                                                                                                                                         acCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCGAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAGACCTGCATGCCGGACTCATctcagaat            >  2:292731/1‑142 (MQ=255)                                                                                                                                         acCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAGACCTGCATGCCGGACTCATCGATAAg             <  7:134294/148‑1 (MQ=255)                                                                                                                                         acCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGGGGCTGAACAACGCTTCCAGGCGAGGATAGACCTGCATGCCGGACTCATCGATaa              <  1:247693/147‑1 (MQ=255)                                                                                                                                           ccGGTGGGCAATGCCACGGCCAAGCATTTCATCGATGGTGGGTACCCGGCTGCTTAAAGGCGCCACGGGCCAACTACCCAATTGGCGG‑GGGTGGCGTACCAACGCTTCCAGGCGAGGATAGACCTGCATGCCGGACTCATCGATAAGGc           <  1:91952/149‑1 (MQ=255)                                                                                                                                           ccGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAGACCTGCATGCCGGACTCATCGATAAg             >  8:345380/1‑146 (MQ=255)                                                                                                                                           ccGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAGACCTGCATGCCGGACTCATCGATAAg             <  7:345380/146‑1 (MQ=255)                                                                                                                                           ccGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAGACCTGCATGCCGGACTCATCGATAAGGCa          >  3:130004/1‑149 (MQ=255)                                                                                                                                                 gCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAGACCTGCATGCCGGACTCATCGATAAGGCATTCGt     >  1:42799/1‑148 (MQ=255)                                                                                                                                                 gCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAGACCTGCATGCCGGACTCATCGATAAGGCATTCGt     >  2:33266/1‑148 (MQ=255)                                                                                                                                                 gCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAGACCTGCATGCCGGACTCATCGATAAGGCATTCGt     >  6:229049/1‑148 (MQ=255)                                                                                                                                                 gCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAGACCTGCATGCCGGACTCATCGATAAGGCATTCGTg    <  3:293334/149‑1 (MQ=255)                                                                                                                                                  ccAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAGACCTGCATGCCGGACTCATCGATAAGGCATTCGTGa   <  2:42799/149‑1 (MQ=255)                                                                                                                                                   cAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAGACCTGCATGCCGGACTCATCGATAAGGCATTCGTGAc  >  1:221213/1‑149 (MQ=255)                                                                                                                                                   aaa‑GCCACCGCCAAGCATTTCATCGATGGTTGGTACGCCGCTGCTGATACGCGCCAACGTGCCACTACCCAATTGGCTG‑GGGTGGGTGAACAACGCTTCCAGGCGAGGATAGACCTGCATGCCTGACTCATCGATAAGGCATTCGTGAc  <  3:240543/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GCCTTGCTGCGCAGGCGTTCGGGTGTTTCGTAAAAACCGAAATGCAACCCCGGTGCCTGCGGGGTGCTGGCACCCAGAAAGGCCAGGCCCAATGATGTCTTGCCAATCCCCGAAGGCCCCATCAACAGGCTGACTGACCCGGTGGCCAA‑GCCACCGCCAAGCATTTCATCGAAGGTTGGTACGCCGCTGCTGATACGCGCCAAGGTGCCACTACCCAATTGGCTG‑GGGTGGCTGAACAACGCTTCCAGGCGAGGATAGACCTGCATGCCGGACTCATCGATAAGGCATTCGTGAC  >  NC_002947/4359036‑4359330

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