Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,107,497	C→G	6.7%	A205A (GCG→GCC)	ydiJ ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,107,497	0	C	G	6.7%	192.0 / 3.5	60	A205A (GCG→GCC)	ydiJ	putative oxidoreductase
Reads supporting (aligned to +/- strand):  ref base C (29/27);  new base G (2/2);  total (31/29)
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.50e-01

GAAGCGGCCCTGCTCGTCGCGCAGGTGCGCCAGGTCGTAGCCGGTCAGGCAACGGTTGAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑TTTTACCGTAGGTGCAGCTGCCCTGGCCACTGGCATCCGTGTTGATCA  >  NC_002947/5107353‑5107645                                                                                                                                                 |                                                                                                                                                      gAAGCGGCCCTGCTCGTCGCGCAGGTGCGCCAGGTCGTAGCCGGTCAGGCAACGGTTGAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCcggcgcgg                                                                                                                                                   >  6:186570/1‑149 (MQ=255) ctaccgtcccTGCTCGTCGCCCAGGAGCGCCAGGTCGAAGCCGGTCCGGCAACGGTTGAGCTTGGGGAAGCTGGTTTCAAACAGTTCGGCCTCGGGTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCcggcgcgg                                                                                                                                                   <  7:284144/142‑1 (MQ=255)    gCGGCCCTGCTCGTCGCGCAGGTGCGCCAGGTCGTAGCCGGTCAGGCAACGGTTGAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACCCATCGCCGCACCGGCCcggcgcgg                                                                                                                                                   >  8:284144/1‑146 (MQ=255)    gCGGCCCTGCTCGTCGCGCAGGTGCGCCAGGTCGTAGCCGGTCAGGCAACGGTTGAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGCGCCGGGCcgggtcggcga                                                                                                                                                <  7:198785/149‑2 (MQ=255)     cGGCCCTGCTCGTCGCGCAGGTGCGCCAGGTCGTAGCCGGTCAGGCAACGGTTGAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCa                                                                                                                                               <  7:301349/149‑1 (MQ=255)     cGGCCCTGCTCGTCGCGCAGGTGCGCCAGGTCGTAGCCGGTCAGGCAACGGTTGAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCa                                                                                                                                               <  6:112156/149‑1 (MQ=255)        cccTGCTCGTCGCGCAGGTGCGCCAGGTCGTAGCCGGTCAGGCAACGGTTGAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGc                                                                                                                                            <  4:259202/149‑1 (MQ=255)          cTGCTCGTCGCGCAGGTGCGCCAGGTCGTAGCCGGTCAGGCAACGGTTGAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGc                                                                                                                                            >  8:44785/1‑147 (MQ=255)          cTGCTCGTCGCGCAGGTGCGCCAGGTCGTAGCCGGTCAGGCAACGGTTGAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGc                                                                                                                                            <  7:44785/147‑1 (MQ=255)                         gTGCGCCAGGTCGTAGCCGGTCAGGCAACGGTTGAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCAt                                                                                                                           >  2:31366/1‑149 (MQ=255)                            cgcCAGGTCGTAGCCGGTCAGGCAACGGTTGAGCTTGGGGAAGGTGGTTTCATTCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCAtcgt                                                                                                                        <  7:23101/149‑1 (MQ=255)                               cAGGTCGTAGCCGGTCAGGCAACGGTTGAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGtt                                                                                                                                       >  6:242754/1‑131 (MQ=255)                               cAGGTCGTAGCCGGTCAGGCAACGGTTGAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGa                                                                                                                     >  5:228061/1‑149 (MQ=255)                                 ggTCGTAGCCGGTCAGGCAACGGTTGAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATc                                                                                                                   <  1:38944/149‑1 (MQ=255)                                    cGTAGCCGGTCAGGCAACGGTTGAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCgg                                                                                                                 >  8:243848/1‑148 (MQ=255)                                     gTAGCCGGTCAGGCAACGGTTGAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCgg                                                                                                                 >  2:320556/1‑147 (MQ=255)                                     gTAGCCGGTCAGGCAACGGTTGAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCa                                                                                                               <  7:203993/149‑1 (MQ=255)                                     gTAGCCGGTCAGGCAACGGTTGAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCa                                                                                                               <  7:121176/149‑1 (MQ=255)                                      tAGCCGGTCAGGCAACGGTTGAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAg                                                                                                              >  5:35785/1‑149 (MQ=255)                                        gCCGGTCAGGCAACGGTTGAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGc                                                                                                            <  1:203050/149‑1 (MQ=255)                                          cGGTCAGGCAACGGTTGAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCtg                                                                                                          <  4:161458/149‑1 (MQ=255)                                            gTCAGGCAACGGTTGAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCtgtg                                                                                                        >  7:310593/1‑149 (MQ=255)                                               aGGCAACGGTTGAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAg                                                                                                     <  1:320556/149‑1 (MQ=255)                                               aGGCAACGGTTGAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAg                                                                                                     <  7:243848/149‑1 (MQ=255)                                                        ttGAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGcc                                                                                            <  7:131792/149‑1 (MQ=255)                                                          gAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCa                                                                                           <  8:310593/148‑1 (MQ=255)                                                           aGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGc                                                                                                       <  1:49517/135‑1 (MQ=255)                                                           aGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGc                                                                                                       >  2:49517/1‑135 (MQ=255)                                                           aGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCAc                                                                                          <  2:82237/148‑1 (MQ=255)                                                            gCTTGGGGAAGGTTGTTTCAATCAGTTCGGCCTGGGTTTCCTTAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCAcc                                                                                         <  7:281326/148‑1 (MQ=255)                                                                       gtggtTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCt                                                                             >  2:326517/1‑149 (MQ=255)                                                                         ggtTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCtgt                                                                           >  6:164706/1‑149 (MQ=255)                                                                           tttCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCAGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCtgtg                                                                          >  7:136775/1‑148 (MQ=255)                                                                                   aGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAg                                                                 <  5:186570/149‑1 (MQ=255)                                                                                    gTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAg                                                                 >  8:146920/1‑148 (MQ=255)                                                                                    gTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAg                                                                 >  3:35586/1‑148 (MQ=255)                                                                                         gCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTg                                                           >  4:19526/1‑149 (MQ=255)                                                                                               gTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTcgcg                                                     >  3:138532/1‑149 (MQ=255)                                                                                                    cTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑ttt                                               >  5:166540/1‑149 (MQ=255)                                                                                                    cTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑ttt                                               >  4:266417/1‑149 (MQ=255)                                                                                                    cTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑ttt                                               >  8:10857/1‑149 (MQ=255)                                                                                                    cTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑ttt                                               >  3:220734/1‑149 (MQ=255)                                                                                                     tGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑ttt                                               >  2:148077/1‑148 (MQ=255)                                                                                                     tGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑ttt                                               >  3:273362/1‑148 (MQ=255)                                                                                                     tGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑ttt                                               >  8:43927/1‑148 (MQ=255)                                                                                                     tGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGAGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCAACGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑ttt                                               >  7:116820/1‑148 (MQ=255)                                                                                                       aaTTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑TTTTAc                                            >  2:331365/1‑149 (MQ=255)                                                                                                            cccGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑TTTTACCGTa                                        <  5:164706/148‑1 (MQ=255)                                                                                                             ccGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCCGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCCCCGAGCAGCACGCTGTGCAGCTCCAGCACGGGGTCGCGGG‑TTTTACCGTAgg                                      >  2:52781/1‑149 (MQ=255)                                                                                                             ccGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑TTTTACCGTAg                                       <  6:35785/148‑1 (MQ=255)                                                                                                                   cATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGCTCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑TTTTACCGTAGGTGCAGc                                <  8:138466/149‑1 (MQ=255)                                                                                                                       cGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑TTTTACCGTAGGTGCAGCTGcc                            <  7:10857/149‑1 (MQ=255)                                                                                                                       cGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑TTTTACCGTAGGTGCAGCTGc                             <  4:220734/148‑1 (MQ=255)                                                                                                                            ggctgcgccgcgccgGGCCGGGTCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCGGTGCAGCTCCAGCACGTGGTCGCGGG‑TTTTACCGTAGGTGCAGCTGCCCTGGc                       >  5:340505/6‑149 (MQ=255)                                                                                                                            gcctgcgccgcgccgGGCCGGGTCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑TTTTACCGTAGGTGCAGCTGCCCTgg                        <  6:340505/143‑1 (MQ=255)                                                                                                                             acCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑TTTTACCGTAGGTGCAGCTGCCCTGGcc                      <  1:326517/149‑1 (MQ=255)                                                                                                                              ctgcgccgcgccgGGCCGGGTCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑TTTTACCGTAGGTGCAGCTGCCCTg                         >  2:254102/4‑145 (MQ=255)                                                                                                                              ctgcgccgcgccgGGCCGGGTCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑TTTTACCGTAGGTGCAGCTGCCCTg                         <  1:254102/142‑1 (MQ=255)                                                                                                                              ccTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑TTTTACCGTAGGTGCAGCTGCCCTGGCCa                     <  3:266417/149‑1 (MQ=255)                                                                                                                               cTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑TTTTACCGTAGGTGCAGCTGCCCTGGCCAc                    >  2:236319/1‑149 (MQ=255)                                                                                                                                 cGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAg                                                                       >  6:262199/1‑97 (MQ=255)                                                                                                                                 cGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAg                                                                       <  5:262199/97‑1 (MQ=255)                                                                                                                                 cGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑TGTTACCGTAGGTGCAGCTGCCCTGGCCACTg                  >  8:231432/1‑149 (MQ=255)                                                                                                                                  gCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGGTTTTTACCGTAGGTGCAGCTGCCCTGGCCACt                   <  1:31366/148‑1 (MQ=255)                                                                                                                                        ccggcccggcGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑TTTTACCGTAGGTGCAGCTGCCCTGGCCACTGGCATCCg           >  3:174227/1‑149 (MQ=255)                                                                                                                                             ccggcGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑TTTTACCGTAGGTGCAGCTGCCCTGGCCACTGGCATCCGTGtt       >  2:159524/1‑148 (MQ=255)                                                                                                                                             ccggcGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑TTTTACCGTAGGTGCAGCTGCCCTGGCCACTGGCATCCGTGTTg      <  8:116820/149‑1 (MQ=255)                                                                                                                                              cggcgcggcgCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑TTTTACCGTAGGTGCAGCTGCCCTGGCCACTGGCATCCGTGTTGa     <  6:166540/149‑1 (MQ=255)                                                                                                                                                 cgcggcgcAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑TTTTACCGTAGGTGCAGCTGCCCTGGCCACTGGCATCCGTGTTGATCa  <  6:228061/149‑1 (MQ=255)                                                                                                                                                 |                                                                                                                                                      GAAGCGGCCCTGCTCGTCGCGCAGGTGCGCCAGGTCGTAGCCGGTCAGGCAACGGTTGAGCTTGGGGAAGGTGGTTTCAATCAGTTCGGCCTGGGTTTCCTGAATTTCCCGGGCCATGCGGTACACCTCGCCGACCCGGCCCGGCGCGGCGCAGGCCTGTTCCAGCGCAGCATCGTCGATCGGCAGGCTGTGCAGGCGTTCGCCACCGAGCAGCACGCTGTGCAGCTCCAGCACGTGGTCGCGGG‑TTTTACCGTAGGTGCAGCTGCCCTGGCCACTGGCATCCGTGTTGATCA  >  NC_002947/5107353‑5107645

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