Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,531,234	G→A	8.2%	intergenic (+36/+12)	lpxC → / ← PP_1344	UDP‑3‑O‑[3‑hydroxymyristoyl] N‑acetylglucosamine deacetylase/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,531,234	0	G	A	8.2%	169.0 / 4.7	61	intergenic (+36/+12)	lpxC/PP_1344	UDP‑3‑O‑[3‑hydroxymyristoyl] N‑acetylglucosamine deacetylase/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (30/26);  new base A (3/2);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.73e-01

CAACCAGCTGTTGCGCAAGCTGATTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTCGGCCGCATGCCCGCCACGCTT  >  NC_002947/1531087‑1531380                                                                                                                                                     |                                                                                                                                                    caacCAGCTGTTGCGCAAGCTGATTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCgg                                                                                                                                                    >  1:152296/1‑149 (MQ=255)  aacCAGCTGTTGCGCAAGCTGATTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCggg                                                                                                                                                   >  2:178811/1‑149 (MQ=255)  aacCAGCTGTTGCGCAAGCTGATTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑ATGTTGCGGCAGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCggg                                                                                                                                                   >  2:21514/1‑149 (MQ=255)    ccAGCTGTTGCGCAAGCTGATTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCggg                                                                                                                                                   <  2:23262/147‑1 (MQ=255)    ccAGCTGTTGCGCAAGCTGATTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTg                                                                                                                                                 >  2:15791/1‑149 (MQ=255)    ccAGCTGTTGCGCAAGCTGATTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATCAAGGCCACCTTCGGGTg                                                                                                                                                 >  4:28390/1‑149 (MQ=255)    ccAGCTGTTGCGCAAGCTGATTGCTGAAAACGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCggg                                                                                                                                                   >  1:23262/1‑147 (MQ=255)             gcgcAAGCTGATTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCttttt                                                                                                                                         <  2:152296/148‑1 (MQ=255)                       aTTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATgg                                                                                                                                    >  6:6705/1‑143 (MQ=255)                       aTTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATgg                                                                                                                                    <  5:6705/143‑1 (MQ=255)                       aTTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCCGAAGGTGGCCttaatgaa                                                                                                                                      >  6:43011/1‑135 (MQ=255)                       aTTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCCGAAGGTGGCCttaatgaa                                                                                                                                      <  5:43011/141‑7 (MQ=255)                           cTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTgg                                                                                                                          <  5:108030/149‑1 (MQ=255)                            tGAAACCGATGCCT‑GGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGc                                                                                                                        >  2:186760/1‑149 (MQ=255)                              aaaCCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGc                                                                                                                        <  2:173690/148‑1 (MQ=255)                                    aTGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGgcc                                                                                                                 >  7:233289/1‑149 (MQ=255)                                        cTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGgccgcc                                                                                                              >  4:137033/1‑148 (MQ=255)                                         tGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCaa                                                                                                            >  1:269974/1‑149 (MQ=255)                                                 ggTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTc                                                                                                      >  5:20697/1‑147 (MQ=255)                                                  gTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCCACTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTc                                                                                                      >  1:70553/1‑147 (MQ=255)                                                  gTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGgcc                                                                                                                 <  1:124207/135‑1 (MQ=255)                                                  gTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCa                                                                                                    <  3:137033/148‑1 (MQ=255)                                                  gTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAg                                                                                                   >  4:178910/1‑149 (MQ=255)                                                      cTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCg                                                                                               <  2:70553/149‑1 (MQ=255)                                                         tCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCa                                                                                            >  8:274351/1‑149 (MQ=255)                                                            aaGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCCCCTTCGGGTGGCCGTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACg                                                                                          >  1:235451/1‑148 (MQ=255)                                                                            aCCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCg                                                                             <  4:287358/145‑1 (MQ=255)                                                                            aCCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCg                                                                             >  3:287358/1‑145 (MQ=255)                                                                             ccGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAAt                                                                        <  3:178910/149‑1 (MQ=255)                                                                                     gTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCAc                                                                                           >  3:102507/1‑122 (MQ=255)                                                                                     gTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCAc                                                                                           <  4:102507/122‑1 (MQ=255)                                                                                         aTGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCt                                                                                                                            <  1:146317/85‑1 (MQ=255)                                                                                         aTGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCt                                                                                                                            >  2:146317/1‑85 (MQ=255)                                                                                                 cTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTcgcg                                                      >  5:157170/1‑148 (MQ=255)                                                                                                  tGTTGCGGCCGTGTAAGTTAGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGACGCAAACCGTTCCAGTGCGGCAAGCAGTTTGGGGTCGGTAATGACCTcgg                                                                >  3:206249/1‑137 (MQ=255)                                                                                                   gTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAg                                                   >  5:150854/1‑149 (MQ=255)                                                                                                      gCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGct                                                >  5:243170/1‑149 (MQ=255)                                                                                                         gCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCTGTAATGCCCTCGGCGGTGTCGCGCAGGCtctc                                             <  6:243170/149‑1 (MQ=255)                                                                                                           cGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCgg                                           <  5:210420/149‑1 (MQ=255)                                                                                                                  gTTCGAACACTCTCTAGTTCATCAAGGCCACCTTCGGGTTGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCt                                    <  3:28390/149‑1 (MQ=255)                                                                                                                      gAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGtt                                 <  6:157170/148‑1 (MQ=255)                                                                                                                        acacTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGGCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGa                              <  3:278432/149‑1 (MQ=255)                                                                                                                            tctctAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACAGTACCAGTGCGGCAGGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCCAA‑AGc                         >  4:219651/1‑149 (MQ=255)                                                                                                                              tctAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTc                       >  4:110118/1‑149 (MQ=255)                                                                                                                                tAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCTCTGTTCGAA‑AGCTc                       <  1:21514/147‑1 (MQ=255)                                                                                                                                  gTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCcgcg                                      >  7:235791/1‑131 (MQ=255)                                                                                                                                  gTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCcgcg                                      <  8:235791/131‑1 (MQ=255)                                                                                                                                   ttCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGtt                                 <  8:130038/135‑1 (MQ=255)                                                                                                                                   ttCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGtt                                 >  7:130038/1‑135 (MQ=255)                                                                                                                                   ttCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTCGGCCg                  >  4:54436/1‑149 (MQ=255)                                                                                                                                    taaaaaaGGCCACCCGAAGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAg                                                   <  1:79693/111‑1 (MQ=255)                                                                                                                                    taaaaaaGGCCACCCGAAGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAg                                                   >  2:79693/6‑116 (MQ=255)                                                                                                                                     cATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTCGGCCGc                 <  7:274351/148‑1 (MQ=255)                                                                                                                                     cATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTCGGCCGCa                >  1:121823/1‑149 (MQ=255)                                                                                                                                      aTTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTCGGCCGCa                <  1:186760/148‑1 (MQ=255)                                                                                                                                      aTTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTCGGCCGCAt               <  2:269974/149‑1 (MQ=255)                                                                                                                                      aTTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGAGCT‑TTCGAACAGCGCGGCCgaga               <  8:233289/149‑4 (MQ=255)                                                                                                                                       ttAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCGAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTCGGCCGCGTg              >  8:300053/1‑149 (MQ=255)                                                                                                                                       ttAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTc                                >  4:140726/1‑132 (MQ=255)                                                                                                                                       ttAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTc                                <  3:140726/132‑1 (MQ=255)                                                                                                                                         aaGGCCACCCGAAGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTCGGCCGCATGcc            >  6:67887/1‑149 (MQ=255)                                                                                                                                          aGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTCGGCCGCATGcc            <  6:20697/148‑1 (MQ=255)                                                                                                                                           ggCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTCGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGATTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCAGTTCGAA‑AGCTCGGCCGCATGCCCg          >  6:158468/1‑149 (MQ=255)                                                                                                                                                 cTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTCGGCCGCATGCCCGCCACGc    >  6:82707/1‑149 (MQ=255)                                                                                                                                                 cTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGGCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTCGGCCGCATGCCCGCCACGc    >  6:262376/1‑149 (MQ=255)                                                                                                                                                    cGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTCGGCCGCATGCCCGCCACGCtt  <  3:110118/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CAACCAGCTGTTGCGCAAGCTGATTGCTGAAACCGATGCCTGGGAAGTGGTCACTTTCGAAGATGCCAGCACGGCACCGATCTCGTACATGCGCC‑CTGTTGCGGCCGTGTAAGTTCGAACACTCTCTAGTTCATTAAGGCCACCTTCGGGTGGCCTTTTTTATGGCTTGCCCTGGGCGTGGGCCGCCAACCGTTCCAGTGCGGCACGCAGTTTGGGGTCGGTAATGCCCTCGGCGGTGTCGCGCAGGCTCTCGGCCGCGCTGTTCGAA‑AGCTCGGCCGCATGCCCGCCACGCTT  >  NC_002947/1531087‑1531380

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