Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,586,528	A→T	6.8%	L253Q (CTG→CAG)	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,528	0	A	T	6.8%	268.8 / 7.2	88	L253Q (CTG→CAG)	PP_4914	phytase‑like protein
Reads supporting (aligned to +/- strand):  ref base A (42/40);  new base T (3/3);  total (45/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 = 6.18e-01

ATCACCAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTCGATG  >  NC_002947/5586381‑5586674                                                                                                                                                      |                                                                                                                                                    aTCACCAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAg                                                                                                                                                    >  6:389300/1‑149 (MQ=255) aTCACCAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAg                                                                                                                                                    >  5:382690/1‑149 (MQ=255)    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                                                                                                                                         >  6:391949/1‑143 (MQ=255)                  caggccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGa                                                                                                                                         <  5:391949/143‑1 (MQ=255)                   aggccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGTGTAAACAGCTTGCCCTCGAAccaggc                                                                                                                                  <  5:311765/149‑1 (MQ=255)                     gccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAAccaggc                                                                                                                                  >  7:398539/1‑147 (MQ=255)                     gccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGCAAGCTGTTTACCCTCGAAc                                                                                                                                       >  2:242658/1‑142 (MQ=255)                     gccaAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGCAAGCTGTTTACCCTCGAAc                                                                                                                                       <  1:242658/142‑1 (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                                                                                <  4:271186/149‑1 (MQ=255)                                                                      gAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCgg                                                                                >  7:196076/1‑149 (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)                                                                                                                                                      |                                                                                                                                                    ATCACCAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCTCGATG  >  NC_002947/5586381‑5586674

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