Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,433,326:2	+A	6.8%	coding (119/351 nt)	PP_5648 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,433,326	2	.	A	6.8%	231.3 / 10.1	59	coding (119/351 nt)	PP_5648	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base . (23/32);  new base A (2/2);  total (25/34)
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.50e-01

ACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAA  >  NC_002947/4433180‑4433469                                                                                                                                                     |                                                                                                                                                 aCTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGTTGCCCAGTGCATCCAGtt‑‑‑ca                                                                                                                                               <  4:199512/149‑2 (MQ=255)  cTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTg                                                                                                                                              <  1:193320/149‑1 (MQ=255)     cGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGgcg                                                                                                                                           >  6:167256/1‑149 (MQ=255)       ggTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGcc                                                                                                                                         >  3:8719/1‑149 (MQ=255)       ggTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGcc                                                                                                                                         <  5:167256/149‑1 (MQ=255)           aaCGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCCGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTtgt                                                                                                                                     <  2:80889/149‑1 (MQ=255)             cGTGCACGAACGCTTCGGCCTGAGGGGTAATCACGGTGAATCCGCCCTGCTTGCGCGGTGCCAGGTAGCCGCTGGCGTCATCACCGAACACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCGTTGTGc                                                                                                                                   <  6:252553/149‑1 (MQ=255)             cGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGc                                                                                                                                   <  6:49901/149‑1 (MQ=255)               tCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGa                                                                                                                                 <  2:151095/149‑1 (MQ=255)                         cTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGc                                                                                                                                   >  7:149444/1‑137 (MQ=255)                         cTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGc                                                                                                                                   <  8:149444/137‑1 (MQ=255)                         cTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTgc                                                                                                                       <  5:84667/149‑1 (MQ=255)                              ccgatctgggTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGt                                                                                                                  <  1:122422/142‑1 (MQ=255)                                cTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCg                                                                                                                >  3:168789/1‑149 (MQ=255)                                 tCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGt                                                                                                                               <  7:282745/133‑1 (MQ=255)                                 tCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGt                                                                                                                               >  8:282745/1‑133 (MQ=255)                                 tCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCgg                                                                                                               >  7:97221/1‑149 (MQ=255)                                    cGGGTAAACACGGTGATTCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCtt                                                                                                            <  1:15324/149‑1 (MQ=255)                                    cGGGTAAACACCGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCtt                                                                                                            <  3:196868/149‑1 (MQ=255)                                      ggTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGt                                                                                                                  >  2:122422/1‑141 (MQ=255)                                            cacGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCAc                                                                                                    >  1:845/1‑149 (MQ=255)                                             acGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTgcc                                                                                                         <  1:24974/143‑1 (MQ=255)                                             acGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTgcc                                                                                                         >  2:24974/1‑143 (MQ=255)                                                gTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTg                                                                                                           >  4:85483/1‑138 (MQ=255)                                                gTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTg                                                                                                           <  3:85483/138‑1 (MQ=255)                                                     tCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGttt                                                                                           >  5:255050/1‑149 (MQ=255)                                                       cGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTc                                                                                          <  6:264573/148‑1 (MQ=255)                                                        gCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCa                                                                                        <  8:97221/149‑1 (MQ=255)                                                                   gcgGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAg                                                                               >  1:225415/1‑147 (MQ=255)                                                                   gcgGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGaa                                                                             <  6:255050/149‑1 (MQ=255)                                                                        tCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGTTCAGCTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAg                                                                            >  7:288500/1‑148 (MQ=255)                                                                        tCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGTTCAGCTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAg                                                                            >  1:228229/1‑148 (MQ=255)                                                                          cAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCg                                                                                  >  6:43813/1‑137 (MQ=255)                                                                          cAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCg                                                                                  <  5:43813/137‑1 (MQ=255)                                                                               aGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGcc                                                                 >  6:156695/1‑149 (MQ=255)                                                                                gCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCa                                                                <  4:168789/149‑1 (MQ=255)                                                                                 ccGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGTTCAGCTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGc                                                                  <  2:228229/149‑1 (MQ=255)                                                                                 ccGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGTTCAGCTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGc                                                                  <  8:288500/149‑1 (MQ=255)                                                                                       gCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAg                                                         >  1:190679/1‑149 (MQ=255)                                                                                             tcaCCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCg                                                    >  1:261107/1‑148 (MQ=255)                                                                                               aCCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCa                                                  <  5:156695/148‑1 (MQ=255)                                                                                               aCCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAg                                                 <  4:8719/149‑1 (MQ=255)                                                                                                ccGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCa                                                          >  8:270160/1‑139 (MQ=255)                                                                                                ccGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCa                                                          <  7:270160/139‑1 (MQ=255)                                                                                                    aaaCCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAAc                                            >  1:265443/1‑149 (MQ=255)                                                                                                     aaCCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACg                                           <  7:312163/149‑1 (MQ=255)                                                                                                     aaCCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACg                                           <  2:223442/149‑1 (MQ=255)                                                                                                       ccTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGc                                         <  2:225415/149‑1 (MQ=255)                                                                                                        cTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACgg                                          >  2:105599/1‑147 (MQ=255)                                                                                                           ctctTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGa                                      <  1:105599/148‑1 (MQ=255)                                                                                                                    ccAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATcgcacg                              <  2:845/147‑1 (MQ=255)                                                                                                                           aaCTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATCTTc                       <  2:265443/147‑1 (MQ=255)                                                                                                                              tGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGtt                        >  1:316615/1‑143 (MQ=255)                                                                                                                              tGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGtt                        <  2:316615/143‑1 (MQ=255)                                                                                                                               ggATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCt                 >  3:252477/1‑149 (MQ=255)                                                                                                                                gATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATcgc                                 >  6:293030/1‑132 (MQ=255)                                                                                                                                gATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATcgc                                 <  5:293030/132‑1 (MQ=255)                                                                                                                                 aTGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCg               >  5:314828/1‑149 (MQ=255)                                                                                                                                              ccAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCaaaaaa  <  6:108188/149‑1 (MQ=255)                                                                                                                                              ccAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCaaaaaa  <  4:252477/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                 ACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAA  >  NC_002947/4433180‑4433469

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