Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,400,872	C→A	6.2%	intergenic (+35/‑134)	ybgF → / → queE	TolQRA transport system periplasmic protein/7‑carboxy‑7‑deazaguanine synthase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,400,872	0	C	A	6.2%	209.0 / 3.8	65	intergenic (+35/‑134)	ybgF/queE	TolQRA transport system periplasmic protein/7‑carboxy‑7‑deazaguanine synthase
Reads supporting (aligned to +/- strand):  ref base C (26/35);  new base A (2/2);  total (28/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.11e-01

GAGCGCCGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGACACATTA  >  NC_002947/1400724‑1401020                                                                                                                                                     |                                                                                                                                                      gAGCGCCGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAA‑AAACCCgcgct                                                                                                                                                      >  1:291085/1‑148 (MQ=255)    cgcCGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCATGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAg                                                                                                                                                    >  7:153479/1‑149 (MQ=255)             ggCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGtt                                                                                                                                           >  6:199533/1‑149 (MQ=255)                 aTACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGttt‑ttt                                                                                                                                      <  1:326150/149‑1 (MQ=255)                    ccGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAA‑AAACCCGCGACTGGCGCGGGTTTCTTTCg                                                                                                                                    <  2:291085/148‑1 (MQ=255)                      gACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTa                                                                                                                                 <  6:183043/149‑1 (MQ=255)                       aCAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGG‑CCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAg                                                                                                                                >  5:250745/1‑148 (MQ=255)                            gTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGG‑CCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCAc                                                                                                                          <  6:250745/149‑1 (MQ=255)                             tAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCAc                                                                                                                          >  1:76077/1‑149 (MQ=255)                             tAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCAc                                                                                                                          >  2:178920/1‑149 (MQ=255)                             tAAAAGGCATCCTTCAGCAGGTCATGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCAc                                                                                                                          <  8:153479/149‑1 (MQ=255)                               aaaGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACt                                                                                                                         <  7:174695/148‑1 (MQ=255)                                   gCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCAc                                                                                                                          >  8:113857/1‑143 (MQ=255)                                   gCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCAc                                                                                                                          <  7:113857/143‑1 (MQ=255)                                     aTCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCttt                                                                                                                  >  8:144924/1‑149 (MQ=255)                                         ttCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCATTAGAATCACTGCCCTTTTTTc                                                                                                              >  8:246490/1‑149 (MQ=255)                                           cagcagGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCg                                                                                                             <  2:76077/148‑1 (MQ=255)                                                ggTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAaca                                                                                                       <  3:103252/149‑1 (MQ=255)                                                    gTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACgct                                                                                                   <  6:260985/149‑1 (MQ=255)                                                      gACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACgctgc                                                                                                 >  2:241215/1‑149 (MQ=255)                                                        cccAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGctgctg                                                                                               >  8:243937/1‑149 (MQ=255)                                                         ccAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCtgctgc                                                                                              <  7:1743/149‑1 (MQ=255)                                                         ccAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGGACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCtgctgc                                                                                              >  2:296225/1‑149 (MQ=255)                                                           aGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCgg                                                                                            <  6:96741/149‑1 (MQ=255)                                                               ccccGGCACCTCTGGCGCTCTACTGCCCCAGCGTCAGCTGCAGAATCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCgg                                                                                                                                              <  2:198909/96‑1 (MQ=255)                                                                   ggCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCt                                                                                    <  1:178920/149‑1 (MQ=255)                                                                       ccTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGt                                                                                <  7:144924/149‑1 (MQ=255)                                                                        ctctGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTc                                                                               >  6:200353/1‑149 (MQ=255)                                                                         tctGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTc                                                                               >  7:195732/1‑148 (MQ=255)                                                                                cgcAACTGGCGCAGGGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAAGCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAAGACGCTGCTGCGGATAACGCTATGTCGAGTCcgc                                                                       <  8:229876/147‑1 (MQ=255)                                                                                   aaCTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGa                                                                     >  8:310621/1‑148 (MQ=255)                                                                                   aaCTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGa                                                                     <  5:199533/148‑1 (MQ=255)                                                                                        gCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTg                                                                >  4:78382/1‑148 (MQ=255)                                                                                        gCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGc                                                               <  5:200353/149‑1 (MQ=255)                                                                                          ccAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCt                                                             <  7:243937/149‑1 (MQ=255)                                                                                          ccAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCt                                                             >  7:103801/1‑149 (MQ=255)                                                                                          ccAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCt                                                             >  6:278134/1‑149 (MQ=255)                                                                                          ccAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCt                                                             >  8:217418/1‑149 (MQ=255)                                                                                               acgACCTGCTGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACg                                                         <  4:243915/146‑1 (MQ=255)                                                                                                tGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCATTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGa                                                       <  7:246490/149‑1 (MQ=255)                                                                                                  aCCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTCCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAgg                                                     <  1:70128/149‑1 (MQ=255)                                                                                                   ccTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGc                                                    <  8:195732/149‑1 (MQ=255)                                                                                                   ccTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGc                                                    <  7:310621/149‑1 (MQ=255)                                                                                                       cAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGa                                                 >  2:200894/1‑148 (MQ=255)                                                                                                         gAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAg                                              >  3:100305/1‑149 (MQ=255)                                                                                                             ctctAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTg                                          <  3:78382/149‑1 (MQ=255)                                                                                                                  ggCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGttt                                       <  3:158297/147‑1 (MQ=255)                                                                                                                  ggCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGtt                                        <  5:278134/146‑1 (MQ=255)                                                                                                                    cGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCt                                   <  4:100305/149‑1 (MQ=255)                                                                                                                     gTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGc                                                    <  5:192526/131‑1 (MQ=255)                                                                                                                     gTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGc                                                    >  6:192526/1‑131 (MQ=255)                                                                                                                        gTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCa                               >  6:21442/1‑149 (MQ=255)                                                                                                                         tCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACg                                                         <  4:179484/122‑1 (MQ=255)                                                                                                                               tACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGt                        <  8:193095/149‑1 (MQ=255)                                                                                                                                    cTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGa                   <  1:296225/149‑1 (MQ=255)                                                                                                                                    cTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGa                   <  1:135741/149‑1 (MQ=255)                                                                                                                                       gaaagaaaCCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGc                 <  3:162859/148‑1 (MQ=255)                                                                                                                                       gaaagaaaCCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGc                 >  4:162859/1‑148 (MQ=255)                                                                                                                                        gaaaaaaCCCGCGACTGGCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCAt               >  4:173700/5‑149 (MQ=255)                                                                                                                                         aaaaaaCCCGCGACTGGCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGa                                                 <  3:248424/111‑1 (MQ=255)                                                                                                                                         aaaaaaCCCGCGACTGGCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGa                                                 >  4:248424/4‑114 (MQ=255)                                                                                                                                           gaaaCCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGc                 <  2:277299/144‑1 (MQ=255)                                                                                                                                           gaaaCCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGc                 >  1:277299/1‑144 (MQ=255)                                                                                                                                            aaaCCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGt                                                                                >  8:134983/1‑80 (MQ=255)                                                                                                                                              aCCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATTCAAGa         <  5:137844/149‑1 (MQ=255)                                                                                                                                              aCCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGa         >  4:300278/1‑149 (MQ=255)                                                                                                                                                  gcgcCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGAcaca     >  5:256997/1‑149 (MQ=255)                                                                                                                                                   cgcCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGACACAt    <  1:200894/149‑1 (MQ=255)                                                                                                                                                    gcCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTg                       >  6:342980/1‑129 (MQ=255)                                                                                                                                                     ccAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGACACATTa  <  7:217418/149‑1 (MQ=255)                                                                                                                                                     ccAGTCGCCGGTTT‑TTTCCTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGACACATTa  >  2:320388/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                      GAGCGCCGTATGGGCCATACCGACAAGGTAAAAGGCATCCTTCAGCAGGTCGTGACCCAGTACCCCGGCACCTCTGCCGCTCAACTGGCCCAGCGTGACCTGCAGAAGCTCTAGGCGTTGTCGCTGTACCGCTCGAAAGAAACCCGCGCCAGTCGCGGGTTT‑TTTCGTTAGAATCACTGCCCTTTTTTCGTAAACACGCTGCTGCGGATAACGCTATGTCGAGTCCGCGACAGTGCCTGACGGAGGCGGACAGCCTGTTTAGCTGTCACGCCCGTGGCGAGCATGCAAGACACATTA  >  NC_002947/1400724‑1401020

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