Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,809,182	T→A	11.1%	L78Q (CTG→CAG)	eno →	enolase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,809,182	0	T	A	11.1%	148.0 / 8.4	54	L78Q (CTG→CAG)	eno	enolase
Reads supporting (aligned to +/- strand):  ref base T (21/27);  new base A (3/3);  total (24/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

ATCATCGGCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCCCAGGACCAGGAC  >  NC_002947/1809040‑1809330                                                                                                                                               |                                                                                                                                                     atcatcGGCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACctgctgct                                                                                                                                                >  1:312817/1‑149 (MQ=255) atcatcGGCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACctgctgct                                                                                                                                                >  6:171535/1‑149 (MQ=255)    atcGGCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCAGTGACCTGCTGCTgg                                                                                                                                              <  1:366836/148‑1 (MQ=255)    atcGGCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGACTCGCGCGAAGCGGTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCATCATCAACGGCCCGATCCGTGACCTGCTGCTggg                                                                                                                                             <  1:169212/149‑1 (MQ=255)       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)       ggCAGCGCTTGCGCGCCGTCCGGTGCTTCTACAGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGCCAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGc                                                                                                                                            >  5:285281/1‑147 (MQ=255)        gCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTCGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGCAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAg                                                                                                                                         >  1:272574/1‑149 (MQ=255)         cAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCTTGCCCAGCAGCAGGTCacg                                                                                                                                         >  4:151667/1‑146 (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                                                                                                                                 >  8:39391/1‑144 (MQ=255)                     gcCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTc                                                                                                                                 <  7:39391/144‑1 (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)                                                                                                                                               |                                                                                                                                                     ATCATCGGCAGCGCTTGCGCGCCGTCCGGTGCTTCTACCGGCTCGCGCGAAGCGCTGGAGCTGCGTGATGGCGACAAGAGCCGTTACCTGGGCAAGGGCGTGCTGAAGGCCGTCGCCAACATCAACGGCCCGATCCGTGACCTGCTGCTGGGCAAGGATCCTTCCGATCAGAAGGCTCTGGACCGCGCCATGATCGAACTGGACGGTACCGAGAACAAGGCCAAGCTGGGCGCCAACGCCATCCTGGCTGTTTCCCTGGCTGCCGCCAAAGCTGCTGCCCAGGACCAGGAC  >  NC_002947/1809040‑1809330

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