Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,691,772	T→C	6.3%	intergenic (+39/‑132)	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,772	0	T	C	6.3%	187.9 / 2.5	63	intergenic (+39/‑132)	PP_1487/PP_1488	hypothetical protein/methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base T (28/31);  new base C (2/2);  total (30/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.16e-01

AGGAAAACTTCATGGAGCGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCT  >  NC_002947/1691624‑1691919                                                                                                                                                     |                                                                                                                                                     aGGAAAACTTCATGGAGCGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGAt                                                                                                                                                      <  2:4458/149‑1 (MQ=255)   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                                                                                            >  8:117637/1‑119 (MQ=255)                                                                                        aGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATca                                                                                            <  7:117637/119‑1 (MQ=255)                                                                                             cGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCgg                                                         <  2:326986/149‑1 (MQ=255)                                                                                             cGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCgg                                                         >  6:229161/1‑149 (MQ=255)                                                                                             cGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCgg                                                         <  2:6711/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                                                                                                                   >  6:326170/1‑64 (MQ=255)                                                                                                                        ttGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAg‑g                                                                                                                   <  5:326170/64‑1 (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                               >  4:267718/1‑137 (MQ=255)                                                                                                                                   tATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTc                               <  3:267718/137‑1 (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                               >  4:131311/1‑127 (MQ=255)                                                                                                                                             cTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTc                               <  3:131311/127‑1 (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)                                                                                                                                                     |                                                                                                                                                     AGGAAAACTTCATGGAGCGTGCCTACCTGGACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCT  >  NC_002947/1691624‑1691919

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