Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,428,786	G→T	6.3%	S98* (TCG→TAG)	yedA ←	phenylalanine/threonine exporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,428,786	0	G	T	6.3%	203.0 / 3.2	63	S98* (TCG→TAG)	yedA	phenylalanine/threonine exporter
Reads supporting (aligned to +/- strand):  ref base G (35/24);  new base T (2/2);  total (37/26)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.40e-01

AGCAGGCCGATGCCCACCAGCCCAAGGGCTATCCCTGCCCATTCCAGCTTCGAATTGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCAGGTAGGTGCCGGCACGCCGCGCCAACGCAGGAAGCCGTACAGCAGCCCACCG  >  NC_002947/5428657‑5428932                                                                                                                                  |                                                                                                                                                   agcagGCCGATGCCCACCAGCCCAAGGGCTATCCATGCCCATTACAGCTTCGAATTGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTc                                                                                                                                 >  1:50504/1‑149 (MQ=255)    agGCCGATGCCCACCAGCCCAAGGGCTATCCCTGCCCATTCCAGCTTCGAATTGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCg                                                                                                                                >  4:319904/1‑147 (MQ=255)    agGCCGATGCCCACCAGCCCAAGGGCTATCCCTGCCCATTCCAGCTTCGAATTGCGGGGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCg                                                                                                                                <  3:319904/147‑1 (MQ=255)      gCCGATGCCCACCAGCCCAAGGGCTATCCCTGCCCATTCCAGCTTGGAATTGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCa                                                                                                                            >  1:258677/1‑149 (MQ=255)              ccaccaGCCCAAGGGCTATCCCTGCCCATTCCAGCTTCGAATTGCGGTGCCCGAACAACAGGCCGACCAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCAcc                                                                                                                    >  5:58329/1‑149 (MQ=255)                      cccaGGGCTATCCCTGCCCATTCCAGCTTCGAATTGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCTCCACGCCCGAAGCCACGTCGGCATGCTCGGCCACGGTCACCCCGCCGtt                                                                                                            <  4:267961/146‑1 (MQ=255)                      ccAAGGGCTATCCCTGCCCATTCCAGCTTCGAATTGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGtt                                                                                                            <  5:26607/149‑1 (MQ=255)                       cAAGGGCTATCCCTGCCCATTCCAGCTTCGAATTGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTg                                                                                                           >  6:171566/1‑149 (MQ=255)                             cTATCCCTGCCCATTCCAGCTTCGAATTGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAg                                                                                                     <  8:193420/149‑1 (MQ=255)                             cTATCCCTGCCCATTCCAGCTTCGAATTGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAg                                                                                                     <  2:258677/149‑1 (MQ=255)                                   cTGCCCATTCCAGCTTCGAATTGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAg                                                                                                <  5:356632/148‑1 (MQ=255)                                   cTGCCCATTCCAGCTTCGAATTGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAg                                                                                                >  6:194484/1‑148 (MQ=255)                                   cTGCCCATTCCAGCTTCGAATTGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAg                                                                                                >  2:124921/1‑148 (MQ=255)                                   cTGCCCATTCCAGCTTCGAATTGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAg                                                                                                >  2:335442/1‑148 (MQ=255)                                   cTGCCCATTCCAGCTTCGAATTGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAg                                                                                                >  2:361673/1‑148 (MQ=255)                                   cTGCCCATTCCAGCTTCGAATTGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAg                                                                                                >  6:356632/1‑148 (MQ=255)                                      cccATTCCAGCTTCGAATTGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCTACAGCTGAGcaa                                                                                             >  2:173284/1‑148 (MQ=255)                                             cAGCTTCGAATTGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGc                                                                                                                                       >  8:272883/1‑99 (MQ=255)                                             cAGCTTCGAATTGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGc                                                                                                                                       <  7:272883/99‑1 (MQ=255)                                             cAGCTTCGAATTGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAAc                                                                                     >  4:335258/1‑149 (MQ=255)                                             cAGCTTCGAATTGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCACGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGaaa                                                                                      <  1:261392/148‑1 (MQ=255)                                                 ttCGAATTGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGc                                                                                                          >  5:337255/1‑124 (MQ=255)                                                 ttCGAATTGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGc                                                                                                          <  6:337255/124‑1 (MQ=255)                                                    gAATTGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGc                                                                              <  3:91998/149‑1 (MQ=255)                                                       ttGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAAc                                                                                     >  7:241589/1‑139 (MQ=255)                                                       ttGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAAc                                                                                     <  8:241589/139‑1 (MQ=255)                                                           ggTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTgcc                                                                       >  5:339348/1‑149 (MQ=255)                                                               cccGGACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGccgccc                                                                   >  4:57276/1‑149 (MQ=255)                                                                  gAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCcgcccgcc                                                                <  3:335258/149‑1 (MQ=255)                                                                    acaacaGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCAt                                                              <  2:50504/149‑1 (MQ=255)                                                                      aacaGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTg                                                            >  6:205175/1‑149 (MQ=255)                                                                       acaGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGc                                                           >  3:41002/1‑149 (MQ=255)                                                                         aGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCgg                                                         <  5:194484/149‑1 (MQ=255)                                                                         aGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCgg                                                         <  1:361673/149‑1 (MQ=255)                                                                          ggCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGc                                                        <  7:213784/149‑1 (MQ=255)                                                                             cGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCag                                                     >  7:129311/1‑149 (MQ=255)                                                                             cGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCag                                                     >  6:341405/1‑149 (MQ=255)                                                                                      aGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCTACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCAGGTAGGTGcc                                            <  1:173284/149‑1 (MQ=255)                                                                                        gTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCAGGTAGGTGCCgg                                          >  1:109272/1‑149 (MQ=255)                                                                                        gTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCAGGTAGGTGCCgg                                          <  3:57276/149‑1 (MQ=255)                                                                                          gAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCTAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCAGGTAGGTGCCGGCa                                        >  6:78558/1‑149 (MQ=255)                                                                                          gAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCTAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCAGGTAGGTGCCGGCa                                        >  6:83879/1‑149 (MQ=255)                                                                                          gAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCAGGTAGGTGCCGGCa                                        >  7:66085/1‑149 (MQ=255)                                                                                           aaaaGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCAGGTAGGTGCCGGCAc                                       >  6:213345/1‑149 (MQ=255)                                                                                           aaaaGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCAGGTAGGTGCCGGCAc                                       >  4:29337/1‑149 (MQ=255)                                                                                                  gCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCaggt                                                   <  1:324283/130‑1 (MQ=255)                                                                                                  gCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCaggt                                                   >  2:324283/1‑130 (MQ=255)                                                                                                        tCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCAGGTAGGTGCCGGCACGCCGCGCCAACg                           >  5:243722/1‑148 (MQ=255)                                                                                                         cGCCACCGCCAAAGCCGCCACGCCCTAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCAGGTAGGTGCCGGCACGCCGCGCCAACGCa                         <  5:78558/149‑1 (MQ=255)                                                                                                          gCCACCGCCAAAGCCGCCACGCCCTAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCAGGTAGGTGCCGGCACGCCGCGCCAACGCa                         <  5:83879/148‑1 (MQ=255)                                                                                                          gCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCTGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCAGGTAGGTGCCGGCACGCCGCGCCAACGCAg                        <  5:205175/149‑1 (MQ=255)                                                                                                          gCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCAGGTAGGTGCCGGCACGCCGCGCCAACGCa                         <  6:339348/148‑1 (MQ=255)                                                                                                           ccACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCAGGTAGGTGCCGGCACGCCGCGCCAACGCAg                        >  6:375327/1‑148 (MQ=255)                                                                                                           ccACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCAGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCAGGTACGTGCCGGCACGCCGCGCCAACGCAg                        >  1:83413/1‑148 (MQ=255)                                                                                                             aCCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCAtt                                                             >  7:100502/1‑109 (MQ=255)                                                                                                             aCCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCAtt                                                             <  8:100502/109‑1 (MQ=255)                                                                                                                 ccAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCAGGTAGGTGCCGGCACGCCGCGCCAACGCAGGAAGCCg                 >  3:13960/1‑149 (MQ=255)                                                                                                                 ccAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCAGGTAGGTGCCGGCACGCCGCGCCAACGCAGGAAGCCg                 <  2:109272/149‑1 (MQ=255)                                                                                                                 ccAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCAGGTAGGTGCCGGCACGCCGCGCCAACGCAGGAAGCCg                 <  1:335442/149‑1 (MQ=255)                                                                                                                  cAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCAGGTAGGTGCCGGCACGCCGCGCCAACGCAGGAAGCCg                 >  7:94770/1‑148 (MQ=255)                                                                                                                  cAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCAGGTAGGTGCCGGCACGCCGCGCCAACGCAGGAAGCCg                 >  4:356119/1‑148 (MQ=255)                                                                                                                      gccgccACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCAGGTAGGTGCCGGCACGCCGCGCCAACGCAGGAAGCCGTAcag            >  4:305249/1‑149 (MQ=255)                                                                                                                                ccGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCAGGTAGGTGCCGGCACGCCGCGCCAACGCAGGAAGCCGTACAGCAGCCCCCCg  >  6:372927/1‑149 (MQ=255)                                                                                                                                  |                                                                                                                                                   AGCAGGCCGATGCCCACCAGCCCAAGGGCTATCCCTGCCCATTCCAGCTTCGAATTGCGGTGCCCGAACAACAGGCCGAACAGCAAGGTGAAAAGCGGCACCGTCGCCACCGCCAAAGCCGCCACGCCCGAAGCCACGCCGGCATGCTCGGCCACGGTCACCCCGCCGTTGCCACAGCTGAGCAACAAGAAACCGATCGCCCCTGCCGCCCGCCATTGCGGCCAGGTAGGTGCCGGCACGCCGCGCCAACGCAGGAAGCCGTACAGCAGCCCACCG  >  NC_002947/5428657‑5428932

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