Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,809,194	A→C	10.0%	K82T (AAG→ACG)	eno →	enolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,809,194	0	A	C	10.0%	167.4 / 5.4	56	K82T (AAG→ACG)	eno	enolase
Reads supporting (aligned to +/- strand):  ref base A (25/25);  new base C (2/4);  total (27/29)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GGCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCCCAGGACCAGGACCTGCCGCTGTAC  >  NC_002947/1809046‑1809342                                                                                                                                                     |                                                                                                                                                     ggCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCaa                                                                                                                                                      <  5:223024/149‑1 (MQ=255) ggCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCaa                                                                                                                                                      <  5:215357/149‑1 (MQ=255) ggCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCaa                                                                                                                                                      <  5:186350/149‑1 (MQ=255) ggCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCaa                                                                                                                                                      <  4:170062/149‑1 (MQ=255)  gCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTCGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGCAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAg                                                                                                                                                     >  1:272574/1‑149 (MQ=255)       gcTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCTTGCCCAGCAGCAGGTCACGGATcg                                                                                                                                                <  3:151667/149‑2 (MQ=255)        cTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCTTGCCCAGCAGCAGGTCACGGATcgg                                                                                                                                               >  2:269606/1‑147 (MQ=255)         ttGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCtt                                                                                                                                              >  6:77220/1‑149 (MQ=255)           gcgcgcCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTc                                                                                                                                             <  7:189553/148‑1 (MQ=255)           gcgcgcCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGACAAGGATCCTTc                                                                                                                                             <  2:312817/148‑1 (MQ=255)            cgcgcCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCTTGCCCAGCAGCAGGTCACGGATcgggcc                                                                                                                                            <  1:269606/148‑6 (MQ=255)            cgcgcCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTcc                                                                                                                                            <  5:218870/148‑1 (MQ=255)              cgcCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGAt                                                                                                                                         <  5:77220/149‑1 (MQ=255)               gcCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTc                                                                                                                                             <  7:39391/144‑1 (MQ=255)               gcCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTc                                                                                                                                             >  8:39391/1‑144 (MQ=255)               gcCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATc                                                                                                                                        <  2:286624/149‑1 (MQ=255)                     cGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAgg                                                                                                                                  >  1:326410/1‑149 (MQ=255)                       gTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTTCGATCAGAAGGct                                                                                                                                >  2:139477/1‑149 (MQ=255)                       gTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGct                                                                                                                                >  7:6363/1‑149 (MQ=255)                           ttCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTg                                                                                                                             >  4:264126/1‑148 (MQ=255)                             cTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGa                                                                                                                           <  8:322205/148‑1 (MQ=255)                                ccGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACcgc                                                                                                                       >  3:206323/1‑149 (MQ=255)                                   gcgcgcgcgAAGCGCGGGAGCGGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAAAATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGcc                                                                                                                    <  6:285281/146‑1 (MQ=255)                                       gcgcgAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGa                                                                                                                <  3:264126/149‑1 (MQ=255)                                             aGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAAc                                                                                                          >  6:205996/1‑149 (MQ=255)                                               cgcTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTg                                                                                                        <  5:270813/149‑1 (MQ=255)                                                gcTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTgg                                                                                                       <  2:326410/149‑1 (MQ=255)                                                  tGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGAc                                                                                                     >  1:339626/1‑149 (MQ=255)                                                      gCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTa                                                                                                 >  7:87550/1‑149 (MQ=255)                                                            tGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCTTGCCCAGCAGCAGGTCACGGATcgg                                                                                                                                               >  8:335040/1‑95 (MQ=25)                                                            tGATGGCGACAAGAGCCGTTACCTGGGAAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCTTGCCCAGCAGCAGGTCACGGATcgg                                                                                                                                               <  7:335040/97‑3 (MQ=25)                                                                   gACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGc                                                                                    >  7:221317/1‑149 (MQ=255)                                                                             gTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGgc                                                                          <  8:221317/149‑1 (MQ=255)                                                                                 ccTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGcc                                                                       >  4:173655/1‑148 (MQ=255)                                                                                     ggCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGc                                                                  >  7:76523/1‑149 (MQ=255)                                                                                      gCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGcc                                                                 <  8:6363/149‑1 (MQ=255)                                                                                          gggCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATcc                                                             >  5:41588/1‑149 (MQ=255)                                                                                             cGTGCTGAAGGTCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTgg                                                          <  1:128275/149‑1 (MQ=255)                                                                                              gTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGc                                                         <  5:205996/149‑1 (MQ=255)                                                                                                    aaGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCAATCCTGGCTGTTTc                                                   >  8:52794/1‑149 (MQ=255)                                                                                                             accAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGc                                                                  >  4:103595/2‑125 (MQ=255)                                                                                                             accAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGc                                                                  <  3:103595/124‑1 (MQ=255)                                                                                                                        cGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAgctg                               <  5:171535/149‑1 (MQ=255)                                                                                                                         ggCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACAGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAgctgc                              >  4:241032/1‑149 (MQ=255)                                                                                                                          gCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTTCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAgctgct                             <  1:139477/149‑1 (MQ=255)                                                                                                                              gATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGccc                         <  6:41588/149‑1 (MQ=255)                                                                                                                              gATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGcc                          >  6:346227/1‑148 (MQ=255)                                                                                                                                    ggACCTGCTGCTGGGCACGGATCCTTCCGCTCAGCAGGCTCTGGCCCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCccaggacc                   <  2:272574/148‑1 (MQ=255)                                                                                                                                     gACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCccaggacca                  <  4:206323/149‑1 (MQ=255)                                                                                                                                      aCCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCATGGCTGCCGCCAAAGCTGCTGCccaggaccag                 >  7:278189/1‑149 (MQ=255)                                                                                                                                        ctgctgctgGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATcc                                                             >  7:296143/1‑103 (MQ=255)                                                                                                                                        ctgctgctgGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATcc                                                             <  8:296143/103‑1 (MQ=255)                                                                                                                                        ctgctgctgGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCccaggaccagga               <  8:87550/149‑1 (MQ=255)                                                                                                                                         tgctgctgGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCccagg                      >  2:159503/1‑141 (MQ=255)                                                                                                                                         tgctgctgGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCccagg                      <  1:159503/141‑1 (MQ=255)                                                                                                                                         tgctgctgGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCCcaggaccaggac              >  6:309572/1‑149 (MQ=255)                                                                                                                                                 ggCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCCCAGGACCAGGACCTGCCGCt      <  3:173655/149‑1 (MQ=255)                                                                                                                                                 ggCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCCCAGGACCAGGACCTGCCGCt      >  3:110274/1‑149 (MQ=255)                                                                                                                                                 ggCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCCCAGGACCAGGACCTGCCGCt      <  2:339626/149‑1 (MQ=255)                                                                                                                                                  gCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCCCAGGACCAGGACCTGCCGCTg     >  4:176290/1‑149 (MQ=255)                                                                                                                                                     aGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGAACCGAGATCAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCCCAGGACCAGGACCTGCCGCTGTAc  >  3:293663/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                     GGCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCCCAGGACCAGGACCTGCCGCTGTAC  >  NC_002947/1809046‑1809342

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