Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,586,533	G→A	7.1%	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	7.1%	251.9 / 5.5	84	G251G (GGC→GGT)	PP_4914	phytase‑like protein
Reads supporting (aligned to +/- strand):  ref base G (38/40);  new base A (3/3);  total (41/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTCGATGGCCAC  >  NC_002947/5586386‑5586679                                                                                                                                                      |                                                                                                                                                    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                                                           <  6:376108/149‑1 (MQ=255)                                                                                          ccACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCt                                                           <  7:19333/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                                 >  1:274184/1‑149 (MQ=255)                                                                                                                    ccGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAgg                                 >  3:143990/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                                                                                                <  7:350762/75‑1 (MQ=255)                                                                                                                                 aCGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTggg                                                                                                >  8:350762/1‑75 (MQ=255)                                                                                                                                 aCGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGAc                      <  5:288907/148‑1 (MQ=255)                                                                                                                                   gTTCGAGGGCAAGCTGTTTACCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAg                            >  4:358796/1‑140 (MQ=255)                                                                                                                                   gTTCGAGGGCAAGCTGTTTACCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAg                            <  3:358796/140‑1 (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                                              <  7:187858/112‑1 (MQ=255)                                                                                                                                             aaaCAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCaa                                              >  8:187858/1‑112 (MQ=255)                                                                                                                                               aCAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTCGATg       >  3:247257/1‑149 (MQ=255)                                                                                                                                                  gCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTCGATGGcc    <  5:154936/149‑1 (MQ=255)                                                                                                                                                      gCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTCGATGGccac  <  4:247689/147‑1 (MQ=255)                                                                                                                                                      |                                                                                                                                                    CAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTCGATGGCCAC  >  NC_002947/5586386‑5586679

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