Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,586,533	G→A	6.3%	G251G (GGC→GGT)	PP_4914 ←	phytase‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,586,533	0	G	A	6.3%	190.1 / 2.3	63	G251G (GGC→GGT)	PP_4914	phytase‑like protein
Reads supporting (aligned to +/- strand):  ref base G (28/31);  new base A (2/2);  total (30/33)
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

CGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTCGATGGCC  >  NC_002947/5586394‑5586677                                                                                                                                            |                                                                                                                                                 cGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCTACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGCAAGCTGTTTACCCTCGaa                                                                                                                                          <  8:247016/148‑1 (MQ=255)    ccaggccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAAccag                                                                                                                                      <  7:121995/149‑1 (MQ=255)      aggccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAAccaggc                                                                                                                                    <  3:179190/149‑1 (MQ=255)      aggccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAAccaggc                                                                                                                                    <  4:361764/149‑1 (MQ=255)         ccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAAccaggccag                                                                                                                                 <  3:50647/149‑1 (MQ=255)                 cTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGa                                                                                                                         <  7:306747/149‑1 (MQ=255)                   gcgcGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGCAAGCTGTTTACCCTCGAAc                                                                                                                                         <  6:197597/131‑1 (MQ=255)                   gcgcGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGCAAGCTGTTTACCCTCGAAc                                                                                                                                         >  5:197597/1‑131 (MQ=255)                      cgTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTcg                                                                                                                    >  7:49853/1‑149 (MQ=255)                      cgTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTcg                                                                                                                    >  5:242266/1‑149 (MQ=255)                          gcgcCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGaa                                                                                                                        >  5:230943/1‑141 (MQ=255)                          gcgcCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGaa                                                                                                                        <  6:230943/141‑1 (MQ=255)                            gcCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGa                                                                                                               >  2:206033/1‑148 (MQ=255)                               gTGACTCGACCAGCGCATCGGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAgg                                                                                                           >  2:366614/1‑149 (MQ=255)                                   cTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTc                                                                                                                                             <  4:254880/111‑1 (MQ=255)                                   cTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTc                                                                                                                                             >  3:254880/1‑111 (MQ=255)                                         cAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGaa                                                                                                                                          >  7:257229/1‑108 (MQ=255)                                         cAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGaa                                                                                                                                          <  8:257229/108‑1 (MQ=255)                                          aGCACATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCAtt                                                                                                <  3:107815/149‑1 (MQ=255)                                            cgcATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTgg                                                                                              <  6:242266/149‑1 (MQ=255)                                             gcATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTggg                                                                                             <  7:160120/149‑1 (MQ=255)                                               aTCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCCATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACCGGGCCCTGGCATTGGGCa                                                                                           >  1:362802/1‑149 (MQ=255)                                                 cagcagCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATc                                                                                         >  6:23459/1‑149 (MQ=255)                                                 cagcagCAAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATc                                                                                         >  8:316397/1‑149 (MQ=255)                                                           cGACCAGCAACGCTCGACCTTGCCACTGTGGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGtt                                                                               <  5:47633/149‑1 (MQ=255)                                                               cAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTggg                                                                                             >  2:35952/1‑131 (MQ=255)                                                               cAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTggg                                                                                             <  1:35952/131‑1 (MQ=255)                                                                aGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTg                                                                              >  6:286193/1‑145 (MQ=255)                                                                aGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTg                                                                              <  5:286193/145‑1 (MQ=255)                                                                 gCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCAttt                                                                         >  6:54106/1‑149 (MQ=255)                                                                    aCGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGcgcg                                                                                   >  5:196961/1‑136 (MQ=255)                                                                    aCGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGcgcg                                                                                   <  6:196961/136‑1 (MQ=255)                                                                        tCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGtgcg                                                                                   <  7:147089/132‑1 (MQ=255)                                                                        tCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGtgcg                                                                                   >  8:147089/1‑132 (MQ=255)                                                                               tGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCg                                                           >  3:102834/1‑149 (MQ=255)                                                                                 ccACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGc                                                          <  1:206033/148‑1 (MQ=255)                                                                                   aCTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTca                                                       <  4:102834/149‑1 (MQ=255)                                                                                     tGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGcc                                                                                                        >  5:377482/1‑98 (MQ=255)                                                                                     tGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGcc                                                                                                        <  6:377482/98‑1 (MQ=255)                                                                                     tGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTcagc                                                     <  7:316397/149‑1 (MQ=255)                                                                                           tatCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACg                                               >  3:14842/1‑149 (MQ=255)                                                                                              cACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTcagc                                                     >  1:197750/1‑140 (MQ=255)                                                                                              cACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTcagc                                                     <  2:197750/140‑1 (MQ=255)                                                                                               aCGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCaa                                            <  4:356790/148‑1 (MQ=255)                                                                                               aCGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCaa                                            >  3:356790/1‑148 (MQ=255)                                                                                               aCGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAAc                                           <  1:366614/149‑1 (MQ=255)                                                                                                     gCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGcc                                                                                                        >  8:292996/1‑82 (MQ=255)                                                                                                     gCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGcc                                                                                                        <  7:292996/82‑1 (MQ=255)                                                                                                     gCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTca                                                       >  6:114176/1‑131 (MQ=255)                                                                                                     gCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTca                                                       <  5:114176/131‑1 (MQ=255)                                                                                                     gCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGcc                                     >  8:9322/1‑149 (MQ=255)                                                                                                     gCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGcc                                     >  8:3765/1‑149 (MQ=255)                                                                                                                 ggCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAg                          <  8:49853/148‑1 (MQ=255)                                                                                                                  gCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGAc                        <  5:23459/149‑1 (MQ=255)                                                                                                                   cGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTgcg                                                                                     <  3:51982/87‑1 (MQ=255)                                                                                                                   cGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTgcg                                                                                     >  4:51982/1‑87 (MQ=255)                                                                                                                        cGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGAct                   <  5:54106/148‑1 (MQ=255)                                                                                                                           tCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGActgctg               >  1:241091/1‑149 (MQ=255)                                                                                                                            cGAGGGCAAGCTGTTTACCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTAAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTgg              >  2:331921/1‑149 (MQ=255)                                                                                                                                ggTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCt          >  5:224629/1‑149 (MQ=255)                                                                                                                                     aCAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTCGATg     <  6:224629/149‑1 (MQ=255)                                                                                                                                      cAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTCGATgg    >  5:124508/1‑149 (MQ=255)                                                                                                                                        gCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTCGATGGcc  <  8:317852/149‑1 (MQ=255)                                                                                                                                            |                                                                                                                                                 CGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCATCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTCGATGGCC  >  NC_002947/5586394‑5586677

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑