Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,433,325	C→T	6.8%	E40E (GAG→GAA)	PP_5648 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,433,325	0	C	T	6.8%	176.3 / 2.0	59	E40E (GAG→GAA)	PP_5648	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (23/32);  new base T (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 = 9.61e-01

CCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAA  >  NC_002947/4433177‑4433469                                                                                                                                                     |                                                                                                                                                    ccgACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGc                                                                                                                                                     >  5:49901/1‑149 (MQ=255)    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                                                                                                                                   <  8:149444/137‑1 (MQ=255)                            cTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGc                                                                                                                                   >  7:149444/1‑137 (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                                                                                                                               >  8:282745/1‑133 (MQ=255)                                    tCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGt                                                                                                                               <  7:282745/133‑1 (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                                                                                                         >  2:24974/1‑143 (MQ=255)                                                acGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTgcc                                                                                                         <  1:24974/143‑1 (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)                                                                                                                                                     |                                                                                                                                                    CCGACTGCGGGTAAACGTCCACGAACGCTTCGGCCTCACGGGTAAACACGGTGAATCCGCCCTGCTTGCGCGGCTCCAGGTAGCCGCTGGCGTCATCACCGAAAACCTTCTCTTCCTGCCAGCTGAACTGGATGCACTGGGCAACCAGCT‑‑‑CTGGCGCCTTGTGCGAGTCGAGCTGCGCAGTCGGCTTGCCGCCACGCGCCGTTTCCATGGTCGCAGAAGCACAACCGGCCAGCATCAGCACCGCAGCCAACGGCAGAATCGCACGCATGTTCCATCCTCGGGGTCGCAAAAAA  >  NC_002947/4433177‑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: ‑