Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,586,530	C→T	6.9%	K252K (AAG→AAA)	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,530	0	C	T	6.9%	262.9 / 5.7	87	K252K (AAG→AAA)	PP_4914	phytase‑like protein
Reads supporting (aligned to +/- strand):  ref base C (40/41);  new base T (3/3);  total (43/44)
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.09e-01

ACCAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTCGATGGCC  >  NC_002947/5586384‑5586677                                                                                                                                                     |                                                                                                                                                     aCCAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGCCCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGCCCTTGCCACTGTCGCTATCACGTCGGCACACCCGGTAGCCGTT‑ACGTTCGAGGGTAAACAGCtt                                                                                                                                                    >  1:67017/1‑149 (MQ=255) aCCAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCtt                                                                                                                                                    <  4:313209/149‑1 (MQ=255)   cAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCAGTCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTg                                                                                                                                                   >  1:337693/1‑149 (MQ=255)     gTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGccc                                                                                                                                                <  3:61610/149‑1 (MQ=255)      tGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCt                                                                                                                                               <  1:377242/149‑1 (MQ=255)        cTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTTCACGTTCGAGGGTAAACAGCTTGCCCt                                                                                                                                               <  1:6157/148‑1 (MQ=255)        cTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTc                                                                                                                                              >  8:19333/1‑148 (MQ=255)         ttCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCg                                                                                                                                             >  2:177662/1‑148 (MQ=255)           cGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAAc                                                                                                                                          >  6:154936/1‑149 (MQ=255)           cGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGCAAGCTGTTTACCCTCGaa                                                                                                                                           <  6:116338/148‑1 (MQ=255)             gccaggccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAAc                                                                                                                                          >  3:271186/1‑147 (MQ=255)              ccaggccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAAccag                                                                                                                                       <  6:382690/149‑1 (MQ=255)               caggccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGa                                                                                                                                            <  5:391949/143‑1 (MQ=255)               caggccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGa                                                                                                                                            >  6:391949/1‑143 (MQ=255)                aggccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGTGTAAACAGCTTGCCCTCGAAccaggc                                                                                                                                     <  5:311765/149‑1 (MQ=255)                  gccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAAccaggc                                                                                                                                     >  7:398539/1‑147 (MQ=255)                  gccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGCAAGCTGTTTACCCTCGAAc                                                                                                                                          <  1:242658/142‑1 (MQ=255)                  gccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGCAAGCTGTTTACCCTCGAAc                                                                                                                                          >  2:242658/1‑142 (MQ=255)                   ccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAAccaggccag                                                                                                                                  >  5:191121/1‑149 (MQ=255)                       aCGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGt                                                                                                                                >  3:50145/1‑147 (MQ=255)                        cGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAg                                                                                                                             <  1:219333/149‑1 (MQ=255)                            tGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCg                                                                                                                           >  1:437467/1‑147 (MQ=255)                               gcgTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTc                                                                                                                      >  2:222476/1‑149 (MQ=255)                                  tAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGATCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTcgcg                                                                                                                   <  2:261886/149‑1 (MQ=255)                                      gcCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGAc                                                                                                               <  7:266004/149‑1 (MQ=255)                                        cGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAg                                                                                                             >  3:124273/1‑149 (MQ=255)                                             cTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGccc                                                                                                        <  1:289378/149‑1 (MQ=255)                                              tCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCt                                                                                                       >  5:376108/1‑149 (MQ=255)                                                 aCCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGc                                                                                                    <  2:437467/149‑1 (MQ=255)                                                      cgga‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTgg                                                                                               <  8:314973/146‑1 (MQ=255)                                                      cgcA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTgg                                                                                               <  4:50145/149‑1 (MQ=255)                                                      cgcA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTgg                                                                                               <  6:191121/149‑1 (MQ=255)                                                           tCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCAt                                                                                           >  1:414428/1‑149 (MQ=255)                                                              gcagcGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTg                                                                                        >  3:247689/1‑149 (MQ=255)                                                                   gAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCGGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCgg                                                                                   >  7:190371/1‑149 (MQ=255)                                                                   gAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCgg                                                                                   >  7:196076/1‑149 (MQ=255)                                                                   gAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCgg                                                                                   <  4:271186/149‑1 (MQ=255)                                                                       gACCAGCATCGGTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTATACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCTCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTg                                                                               <  7:157859/149‑1 (MQ=255)                                                                         ccAGCAACGCTCGACCTTGCCACTGTAGAAATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGATGGTAACCAGCTTGCCCTCCAACCAGGCCAGGTCAGCGAAGTCGCGCGCCAGGGCCCTGGCATTGTGCATCTCCGCCGGTTGCa                                                                             >  2:132331/1‑149 (MQ=255)                                                                           aGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCAtt                                                                           >  5:337217/1‑149 (MQ=255)                                                                             cAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGGT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTc                                                                         <  7:8483/149‑1 (MQ=255)                                                                                cGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGAc                                                                      <  2:337693/149‑1 (MQ=255)                                                                                   tCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCgg                                                                                   <  4:172816/133‑1 (MQ=255)                                                                                   tCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCgg                                                                                   >  3:172816/1‑133 (MQ=255)                                                                                       ccTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑Gcc                                                              >  7:62079/1‑149 (MQ=255)                                                                                        cTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACa                                                                                                              <  4:27134/101‑1 (MQ=255)                                                                                        cTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACa                                                                                                              >  3:27134/1‑101 (MQ=255)                                                                                           gCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCg                                                           <  1:222476/148‑1 (MQ=255)                                                                                            ccACTGTCGATATCACGTCGGCTCACGGGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGGGAAGTCGCGCGAGAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCt                                                         <  2:67017/149‑1 (MQ=255)                                                                                            ccACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCt                                                         <  7:19333/149‑1 (MQ=255)                                                                                            ccACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCt                                                         <  6:376108/149‑1 (MQ=255)                                                                                                  tCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGc                                                                                                                            >  3:144212/1‑77 (MQ=255)                                                                                                  tCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGc                                                                                                                            <  4:144212/77‑1 (MQ=255)                                                                                                  tCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTcagca                                                    <  3:247833/148‑1 (MQ=255)                                                                                                   cGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTcagcag                                                   >  3:62584/1‑148 (MQ=255)                                                                                                   cGATATCACGGCGGCACACCCGGTAGGCGTT‑ACCGTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGGCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCagcagc                                                  <  6:403308/149‑1 (MQ=255)                                                                                                     atatCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCAc                                                <  4:124273/149‑1 (MQ=255)                                                                                                     atatCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCACTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCAc                                                >  1:138783/1‑149 (MQ=255)                                                                                                        tCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGc                                              >  5:199100/1‑148 (MQ=255)                                                                                                          aCGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCTGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCaa                                            <  8:190371/148‑1 (MQ=255)                                                                                                          aCGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTcag                                                      >  6:109685/1‑138 (MQ=255)                                                                                                          aCGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTcag                                                      <  5:109685/138‑1 (MQ=255)                                                                                                          aCGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCaa                                            <  8:196076/148‑1 (MQ=255)                                                                                                            gTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCa                                         >  5:200121/1‑149 (MQ=255)                                                                                                             tCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCAc                                        <  4:62584/149‑1 (MQ=255)                                                                                                                      ccGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAgg                               >  3:143990/1‑149 (MQ=255)                                                                                                                      ccGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAgg                               >  1:274184/1‑149 (MQ=255)                                                                                                                       cGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGt                              >  6:332802/1‑149 (MQ=255)                                                                                                                         gTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCCAGCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTc                             <  1:71829/149‑1 (MQ=255)                                                                                                                             gCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGAc                        >  6:288907/1‑149 (MQ=255)                                                                                                                              cGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCgg                                                                                   <  7:354415/90‑1 (MQ=255)                                                                                                                              cGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCgg                                                                                   >  8:354415/1‑90 (MQ=255)                                                                                                                                   tggTTCGAGGGCAAGCTGTTTACCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGAct                   >  3:338055/3‑149 (MQ=255)                                                                                                                                   aCGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTggg                                                                                              >  8:350762/1‑75 (MQ=255)                                                                                                                                   aCGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTggg                                                                                              <  7:350762/75‑1 (MQ=255)                                                                                                                                   aCGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGAc                    <  5:288907/148‑1 (MQ=255)                                                                                                                                     gTTCGAGGGCAAGCTGTTTACCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAg                          <  3:358796/140‑1 (MQ=255)                                                                                                                                     gTTCGAGGGCAAGCTGTTTACCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAg                          >  4:358796/1‑140 (MQ=255)                                                                                                                                         gAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTgg              <  2:274184/148‑1 (MQ=255)                                                                                                                                          aGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGgcg            >  6:170017/1‑149 (MQ=255)                                                                                                                                          aGGGTAAACAGCTTGCCCTCGAACAAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGgcg            >  4:227944/1‑149 (MQ=255)                                                                                                                                           gggTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGgcgc           >  1:403185/1‑149 (MQ=255)                                                                                                                                            ggTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCt          <  6:200121/149‑1 (MQ=255)                                                                                                                                               aaaCAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCaa                                            >  8:187858/1‑112 (MQ=255)                                                                                                                                               aaaCAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCaa                                            <  7:187858/112‑1 (MQ=255)                                                                                                                                                 aCAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTCGATg     >  3:247257/1‑149 (MQ=255)                                                                                                                                                    gCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTCGATGGcc  <  5:154936/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                     ACCAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTCGATGGCC  >  NC_002947/5586384‑5586677

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