Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,691,773	C→G	6.2%	intergenic (+40/‑131)	PP_1487 → / → PP_1488	hypothetical protein/methyl‑accepting chemotaxis transducer

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,691,773	0	C	G	6.2%	203.0 / 3.4	64	intergenic (+40/‑131)	PP_1487/PP_1488	hypothetical protein/methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base C (28/32);  new base G (2/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.41e-01

GAAAACTTCATGGAGCGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTG  >  NC_002947/1691626‑1691921                                                                                                                                                    |                                                                                                                                                      gAAAACTTCATGGAGCGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCg                                                                                                                                                      <  4:278151/149‑1 (MQ=255)  aaaaCTTCATGGAGCGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATc                                                                                                                                                       >  1:166126/1‑147 (MQ=255)       ttCATGGAGCGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGcgcc                                                                                                                                                >  8:295601/1‑149 (MQ=255)                 gTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTgcgg                                                                                                                                      <  2:328028/149‑1 (MQ=255)                    ccTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGgcgc                                                                                                                                   <  2:166126/149‑1 (MQ=255)                    ccTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGgcgc                                                                                                                                   >  4:210850/1‑149 (MQ=255)                            gACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGtctc                                                                                                                           <  1:168007/149‑1 (MQ=255)                              cGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTcac                                                                                                                         <  4:320101/149‑1 (MQ=255)                               gAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTca                                                                                                                          >  5:243710/1‑147 (MQ=255)                                 gCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAg‑g                                                                                                                     >  1:326986/1‑149 (MQ=255)                                  cTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑gc                                                                                                                    <  8:9789/149‑1 (MQ=255)                                  cTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑gc                                                                                                                    <  4:60281/149‑1 (MQ=255)                                           cAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGCGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTc                                                                                                           <  1:275058/149‑1 (MQ=255)                                                 cGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTgac                                                                                                      <  3:210850/148‑1 (MQ=255)                                                     tGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACaa                                                                                                 >  8:228200/1‑149 (MQ=255)                                                       gACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATc                                                                                               <  7:217430/149‑1 (MQ=255)                                                            gggCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACAt                                                                                           >  7:293499/1‑148 (MQ=255)                                                                    gcCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGt                                                                                  <  7:96685/149‑1 (MQ=255)                                                                     ccAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAGTGCGCAAAATTCCTGACGACAATCACATCAATACCGt                                                                                  <  1:254503/149‑1 (MQ=255)                                                                      cAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTa                                                                                 >  4:51714/1‑148 (MQ=255)                                                                         gAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGAc                                                                             >  8:305665/1‑149 (MQ=255)                                                                             tGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACctct                                                                         >  4:86367/1‑149 (MQ=255)                                                                              ggAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGGGCCTGTGAGACGGAGCGCCGCACGGGCGGCGCTCGATCTCACAg‑g                                                                                                                     <  1:99574/104‑1 (MQ=21)                                                                              ggAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGAGCCTGTGAGACGGAGCGCCGCACGGGCGGCGCTCGATCTCACAg‑g                                                                                                                     >  2:99574/1‑104 (MQ=21)                                                                               gAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACTACAATcac                                                                                             <  4:168660/127‑1 (MQ=255)                                                                               gAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATcac                                                                                             >  3:168660/1‑127 (MQ=255)                                                                                    ccAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCta                                                                   >  5:345089/1‑148 (MQ=255)                                                                                      aGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTga                                                                                                       <  2:311385/110‑1 (MQ=255)                                                                                      aGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTga                                                                                                       >  1:311385/1‑110 (MQ=255)                                                                                      aGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATca                                                                                              <  7:117637/119‑1 (MQ=255)                                                                                      aGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATca                                                                                              >  8:117637/1‑119 (MQ=255)                                                                                           cGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCgg                                                           <  2:6711/149‑1 (MQ=255)                                                                                           cGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCgg                                                           >  6:229161/1‑149 (MQ=255)                                                                                           cGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCgg                                                           <  2:326986/149‑1 (MQ=255)                                                                                             gTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCTGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTg                                                                        <  6:22253/134‑1 (MQ=255)                                                                                             gTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTg                                                                        >  5:22253/1‑134 (MQ=255)                                                                                                         gTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCt                                             >  1:361696/1‑149 (MQ=255)                                                                                                           aGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGC‑AAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGt                                          >  8:1304/1‑149 (MQ=255)                                                                                                             tGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGtt                                         <  7:228200/149‑1 (MQ=255)                                                                                                               gCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTtgt                                       <  6:345089/149‑1 (MQ=255)                                                                                                               gCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTtg                                        <  5:229161/148‑1 (MQ=255)                                                                                                                     cTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTc                                 <  8:293499/149‑1 (MQ=255)                                                                                                                      ttGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAg‑g                                                                                                                     <  5:326170/64‑1 (MQ=255)                                                                                                                      ttGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAg‑g                                                                                                                     >  6:326170/1‑64 (MQ=255)                                                                                                                      ttGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTcc                                >  6:16121/1‑149 (MQ=255)                                                                                                                          cACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTGGTTGTGAATTCCTGAg                            >  7:234762/1‑149 (MQ=255)                                                                                                                          cACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAg                            <  7:311438/149‑1 (MQ=255)                                                                                                                            cGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAggg                          >  4:276563/1‑149 (MQ=255)                                                                                                                              gCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCt                        >  6:245139/1‑149 (MQ=255)                                                                                                                               ccTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCt                        >  1:363234/1‑148 (MQ=255)                                                                                                                                 tATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTc                                 <  3:267718/137‑1 (MQ=255)                                                                                                                                 tATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTc                                 >  4:267718/1‑137 (MQ=255)                                                                                                                                     gaccgcCTGTCAGATCCAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAAtgt                  <  2:361696/145‑1 (MQ=255)                                                                                                                                     gAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAAtgt                  <  7:305665/148‑1 (MQ=255)                                                                                                                                       gcgcCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAACGTgaa               <  5:16121/149‑1 (MQ=255)                                                                                                                                        cgcATGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTgaa               <  3:51714/148‑1 (MQ=255)                                                                                                                                           cTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTc                                 <  3:131311/127‑1 (MQ=255)                                                                                                                                           cTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTc                                 >  4:131311/1‑127 (MQ=255)                                                                                                                                             gtgAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTg         >  6:367959/1‑149 (MQ=255)                                                                                                                                             gtgAGACGGAGCGCCGCACGGGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAAtgt                  <  8:112876/140‑1 (MQ=255)                                                                                                                                             gtgAGACGGAGCGCCGCACGGGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAAtgt                  >  7:112876/1‑140 (MQ=255)                                                                                                                                              tgAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTg         <  3:86367/148‑1 (MQ=255)                                                                                                                                                  aTCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCt    >  7:48453/1‑149 (MQ=255)                                                                                                                                                   tCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAAtgt                  >  3:312660/1‑134 (MQ=255)                                                                                                                                                   tCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAAtgt                  <  4:312660/134‑1 (MQ=255)                                                                                                                                                    cGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTg  >  7:11555/1‑149 (MQ=255)                                                                                                                                                    |                                                                                                                                                      GAAAACTTCATGGAGCGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTG  >  NC_002947/1691626‑1691921

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