Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	468,461	C→T	6.4%	A832A (GCG→GCA)	PP_0386 ←	sensory box protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	468,461	0	C	T	6.4%	280.3 / 6.7	94	A832A (GCG→GCA)	PP_0386	sensory box protein
Reads supporting (aligned to +/- strand):  ref base C (40/48);  new base T (3/3);  total (43/51)
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.31e-01

GGTGTCAACGTCGCGTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCACGAAGTGGGT  >  NC_002947/468314‑468607                                                                                                                                                     |                                                                                                                                                    ggTGTCAACGTCGCGTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTgcg                                                                                                                                                    <  5:233150/149‑1 (MQ=255)     tCAACGTCGCGTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGca                                                                                                                                                  >  5:274690/1‑147 (MQ=255)       aaCGTCGCGTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAgg                                                                                                                                              >  6:306246/1‑149 (MQ=255)        aCGTCGCGTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGt                                                                                                                                             <  4:141509/149‑1 (MQ=255)        aCGTCGCGTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGt                                                                                                                                             <  8:262293/149‑1 (MQ=255)          gTCGCGTACCTGCTCTTTCAGGCGCTGGGCGATGGCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCTCTGTGCCCTGGCGGTTGGACACCTGCGCACAGGTGa                                                                                                                                           <  6:124164/149‑1 (MQ=255)          gTCGCGTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGa                                                                                                                                           <  5:123467/149‑1 (MQ=255)          gTCGCGTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGa                                                                                                                                           <  6:309139/149‑1 (MQ=255)            cgcgTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGt                                                                                                                                         <  7:27910/149‑1 (MQ=255)             gcgTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGCAAGGACCCCGCCCGGGTGGCCAAGGCTGACGTTGATGTGCTTGAGGCGGTCGAGGTCCAGGAACAGCACTCCGCCCTGGCGGTTGGACACCTGCGCACAGGTCAg                                                                                                                                          >  1:73299/1‑148 (MQ=255)             gcgTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGTGCGCAGGTgtcc                                                                                                                                         <  4:322505/148‑4 (MQ=255)             gcgTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGt                                                                                                                                         <  5:76857/148‑1 (MQ=255)             gcgTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGt                                                                                                                                         <  1:339722/148‑1 (MQ=255)             gcgTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTa                                                                                                                                        >  1:71054/1‑149 (MQ=255)              cgTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGa                                                                                                                                           <  6:61840/145‑1 (MQ=255)              cgTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGa                                                                                                                                           >  5:61840/1‑145 (MQ=255)                  ccTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGa                                                                                                                                           <  6:29440/141‑1 (MQ=255)                  ccTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGa                                                                                                                                           >  5:29440/1‑141 (MQ=255)                  ccTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGccc                                                                                                                                   <  7:278519/149‑1 (MQ=255)                   cTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCt                                                                                                                                  >  5:183399/1‑149 (MQ=255)                     gCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTg                                                                                                                                 >  5:85009/1‑148 (MQ=255)                     actctTTCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCCCAGGTGAGTACGCCCTg                                                                                                                                 >  5:224544/2‑148 (MQ=255)                     actctTTCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTg                                                                                                                                 <  6:224544/147‑1 (MQ=255)                       tctTTCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAg                                                                                                                              >  8:35962/1‑149 (MQ=255)                         tttCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGcc                                                                                                                            >  6:230838/1‑149 (MQ=255)                           tCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGc                                                                                                                             <  2:97017/146‑1 (MQ=255)                                cgcTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTc                                                                                                                     <  6:183399/149‑1 (MQ=255)                                cgcTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTc                                                                                                                     <  5:230838/149‑1 (MQ=255)                                        gATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAgg                                                                                                             <  5:306246/149‑1 (MQ=255)                                             ccTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACaa                                                                                                               >  2:35470/1‑142 (MQ=255)                                             ccTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACaa                                                                                                               <  1:35470/142‑1 (MQ=255)                                             ccTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGa                                                                                                        >  6:109425/1‑149 (MQ=255)                                             ccTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGAGGTTGATGTCCTTGAAGCGGTCGAGGTCCAGGAACAGCACTTCGCCCTGGCGGTTGGACAACTGCGAAAACGTGACTACGAGCTGGAGCCGGTGTTCGACCAATGTGCGa                                                                                                        >  8:155306/1‑149 (MQ=255)                                              cTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGAt                                                                                                       >  8:339904/1‑149 (MQ=255)                                               ttGAGCAGC‑AGGTCGCAGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGAtt                                                                                                      >  2:259404/1‑149 (MQ=255)                                                  agcagc‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTgg                                                                                                    <  4:201441/148‑1 (MQ=255)                                                    cagc‑agGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTggg                                                                                                   >  1:217041/1‑147 (MQ=255)                                                    cagc‑agGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTggg                                                                                                   >  5:129245/1‑147 (MQ=255)                                                    cagc‑agGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCa                                                                                                 <  5:109425/149‑1 (MQ=255)                                                       c‑agGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGc                                                                                              <  1:362673/149‑1 (MQ=255)                                                                cGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCTCTGGTGTGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGc                                                                                     >  1:312677/1‑149 (MQ=255)                                                                 gACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCg                                                                                    >  5:90979/1‑149 (MQ=255)                                                                     ggggggCCACGGCTGTCGTTGCTGTGCTTGTCGCGGTCGAGGTCCCGGAACAGCTCTGCGCCCTGGCGGTTGGACTCCTTCGCACAGGGGAGTACGCCCGGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTc                                                                                <  8:265794/149‑1 (MQ=255)                                                                       gTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCgt                                                                              >  4:245877/1‑149 (MQ=255)                                                                        tGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCgtg                                                                             <  3:274041/149‑1 (MQ=255)                                                                           ccAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGcc                                                                                                                            >  5:265047/1‑100 (MQ=255)                                                                           ccAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGcc                                                                                                                            <  6:265047/100‑1 (MQ=255)                                                                            cAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGGCGtggtgg                                                                         <  6:274690/149‑1 (MQ=255)                                                                              aGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGc                                                                       <  3:348806/149‑1 (MQ=255)                                                                              aGGCTGTCGTTGATGTGCTGGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGc                                                                       <  2:73299/149‑1 (MQ=255)                                                                                gCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCtt                                                                     <  2:217041/149‑1 (MQ=255)                                                                                  tGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGa                                                                   <  3:116247/149‑1 (MQ=255)                                                                                   gTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTCCGCCCTGGAGCCGGGTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGGTTGAt                                                                  >  8:266763/1‑149 (MQ=255)                                                                                    tCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGa                                                                   >  3:172081/1‑147 (MQ=255)                                                                                    tCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATa                                                                 <  7:35962/149‑1 (MQ=255)                                                                                      gTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTc                                                                                                                     >  5:9592/1‑96 (MQ=255)                                                                                      gTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTc                                                                                                                     <  6:9592/96‑1 (MQ=255)                                                                                      gTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTc                                                                                                                     <  8:221462/96‑1 (MQ=255)                                                                                      gTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTc                                                                                                                     >  7:221462/1‑96 (MQ=255)                                                                                       ttGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTc                                                              >  3:261516/1‑149 (MQ=255)                                                                                          aTGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTCAAGCGGGTCGTGGTGGGCTTGATAGTCGa                                                            <  1:11802/149‑1 (MQ=255)                                                                                          aTGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTCGAg                                                           <  5:150972/149‑1 (MQ=255)                                                                                            gtgCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGc                                                                                     >  3:26844/1‑121 (MQ=255)                                                                                            gtgCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGc                                                                                     <  4:26844/121‑1 (MQ=255)                                                                                               cTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTg                                                       <  6:85009/148‑1 (MQ=255)                                                                                               cTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTg                                                       <  1:259404/148‑1 (MQ=255)                                                                                                 tGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCg                                                                               <  8:132998/122‑1 (MQ=255)                                                                                                 tGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCg                                                                               >  7:132998/1‑122 (MQ=255)                                                                                                 tGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGcc                                                    >  8:210356/1‑149 (MQ=255)                                                                                                 tGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGcc                                                    >  4:118716/1‑149 (MQ=255)                                                                                                   aaGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTg                                                  <  6:90979/149‑1 (MQ=255)                                                                                                      cGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGc                                               <  2:71054/149‑1 (MQ=255)                                                                                                        gTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGt                                             >  7:92284/1‑149 (MQ=255)                                                                                                         tCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGGACAGGTGAGTACGCCCTGGAGCCGGTTTTCGATCAAGCTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTAGAGTTTGGCCTGGGCGtc                                            >  6:107489/1‑148 (MQ=255)                                                                                                           gAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCt                                          <  7:339904/149‑1 (MQ=255)                                                                                                             ggTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGGGTGCtg                                         >  2:342049/1‑147 (MQ=255)                                                                                                               tCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGa                                                                   <  5:292058/120‑1 (MQ=255)                                                                                                               tCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGa                                                                   >  6:292058/1‑120 (MQ=255)                                                                                                                  aGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCt                                   <  7:210356/149‑1 (MQ=255)                                                                                                                        aGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGAt                             <  2:312677/149‑1 (MQ=255)                                                                                                                          cACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGTGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGt                           >  8:49854/1‑149 (MQ=255)                                                                                                                            cTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGTGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCt                         <  7:49854/149‑1 (MQ=255)                                                                                                                                   cTGGTGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCaaaaa                  >  2:156893/1‑149 (MQ=255)                                                                                                                                     ggCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAAcc                >  5:254084/1‑149 (MQ=255)                                                                                                                                          ttGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCAc           <  4:172081/149‑1 (MQ=255)                                                                                                                                            actcacCTGTGCGCAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTTCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCACg          >  5:260039/4‑148 (MQ=255)                                                                                                                                             gACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCACGaa        <  4:261516/149‑1 (MQ=255)                                                                                                                                              tcacCTGTGCGCAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGc                                                                                     <  6:306830/70‑1 (MQ=255)                                                                                                                                              tcacCTGTGCGCAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGc                                                                                     >  5:306830/2‑71 (MQ=255)                                                                                                                                              tcacCTGTGCGCAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCaaaa                   <  5:68116/136‑1 (MQ=255)                                                                                                                                              tcacCTGTGCGCAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCaaaa                   >  6:68116/2‑137 (MQ=255)                                                                                                                                              acacCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGgag                               <  4:26051/125‑1 (MQ=255)                                                                                                                                                  cTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGga                                <  7:33602/120‑1 (MQ=255)                                                                                                                                                  cTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGga                                >  8:33602/1‑120 (MQ=255)                                                                                                                                                   tGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCACGAAGGGGGt  >  4:141522/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    GGTGTCAACGTCGCGTACCTGCTCTTTCAGGCGCTGGGCGATGCCCTTGAGCAGC‑AGGTCGCCGACCGGGTGGCCAAGGCTGTCGTTGATGTGCTTGAAGCGGTCGAGGTCCAGGAACAGCACTGCGCCCTGGCGGTTGGACACCTGCGCACAGGTGAGTACGCCCTGGAGCCGGTTTTCGAACAAGGTGCGATTGGGCAGGCCGGTC‑AGCGGGTCGTGGTGGGCTTGATAGTCGAGTTTGGCCTGGGCGTGCTTGATGCTGGAGATGTCTGCAAAAACCGCCACGAAGTGGGT  >  NC_002947/468314‑468607

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