Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,691,797	G→A	6.8%	intergenic (+64/‑107)	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,797	0	G	A	6.8%	216.5 / 4.1	73	intergenic (+64/‑107)	PP_1487/PP_1488	hypothetical protein/methyl‑accepting chemotaxis transducer
Reads supporting (aligned to +/- strand):  ref base G (32/36);  new base A (3/2);  total (35/38)
Fisher's exact test for biased strand distribution p-value = 6.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.18e-01

GACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAG  >  NC_002947/1691653‑1691941                                                                                                                                                 |                                                                                                                                                  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                                                                                                                           >  1:311385/1‑110 (MQ=255)                                                           aGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTga                                                                                                                           <  2:311385/110‑1 (MQ=255)                                                           aGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATca                                                                                                                  >  8:117637/1‑119 (MQ=255)                                                           aGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATca                                                                                                                  <  7:117637/119‑1 (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                                                                                                                                         >  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                                                     <  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)                                                                                                                     agacggAGCGCCGCACGGGCGGCGCTCGATCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGGTGGGTGTGACTTCCTGAGGGCTGCAAGTGAAGAACTGGAc                          >  8:343630/6‑149 (MQ=37)                                                                                                                     agacggAGCGCCGCACGGGCGGCGCTCGATCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGAc                          >  1:347304/6‑149 (MQ=38)                                                                                                                       aTCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCt                        >  7:48453/1‑149 (MQ=255)                                                                                                                        tCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAAtgt                                      <  4:312660/134‑1 (MQ=255)                                                                                                                        tCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAAtgt                                      >  3:312660/1‑134 (MQ=255)                                                                                                                         ggAGCGCCGCACGGGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCg                                                                                <  1:337403/90‑1 (MQ=255)                                                                                                                         ggAGCGCCGCACGGGCGGCGCTCCGTCTCACAG‑GCGCAAAATTACTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCg                                                                                >  2:337403/2‑91 (MQ=38)                                                                                                                         cGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTg                      >  7:11555/1‑149 (MQ=255)                                                                                                                          gAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTgc                     <  7:295601/149‑1 (MQ=255)                                                                                                                           aGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTgcg                    >  3:170362/1‑149 (MQ=255)                                                                                                                               ccgccCGTGCGGCGCTCCGTCTCACAG‑GCGC‑AAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAAc               <  7:1304/149‑1 (MQ=255)                                                                                                                               ccgcaCGGGCGGCGCTCGATCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCaa                <  2:347304/149‑1 (MQ=255)                                                                                                                                cgccCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAAc               <  5:367959/149‑1 (MQ=255)                                                                                                                                  cccGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTga                                    >  3:4205/1‑126 (MQ=255)                                                                                                                                  cccGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTga                                    <  4:4205/126‑1 (MQ=255)                                                                                                                                     gTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATc          <  5:245139/149‑1 (MQ=255)                                                                                                                                             cTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAg  >  5:6977/1‑149 (MQ=255)                                                                                                                                                 |                                                                                                                                                  GACGAGCTGGCCCGCCAGCTGCGGCTGGACCCGGGCTTGCGCCAGGAACTGGAAAGCCAGGTGCGGTTGGCTGCAGCGTAGTGGCGGTGCTTGCCACGGCCTATTGAGCGCCTGTGAGATCGAGCGCCGCCCGTGCGGCGCTCCGTCTCACAG‑GCGCAAAATTCCTGACGACAATCACATCAATACCGTAAGACCTCTGGGCTATACTTCGGCGATTTTTCCCGCTCGTTGTGAATTCCTGAGGGCTGAATGTGAAGAACTGGACCTTGCGCCAACGGATCCTGGCAAG  >  NC_002947/1691653‑1691941

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