Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,586,523	T→C	6.9%	T255A (ACC→GCC)	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,523	0	T	C	6.9%	247.1 / 6.7	87	T255A (ACC→GCC)	PP_4914	phytase‑like protein
Reads supporting (aligned to +/- strand):  ref base T (40/41);  new base C (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 = 7.83e-01

TCGATCACCAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCT  >  NC_002947/5586378‑5586669                                                                                                                                                    |                                                                                                                                                    tCGATCACCAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTaaa                                                                                                                                                  <  1:140919/149‑1 (MQ=255) tCGATCACCAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTaaa                                                                                                                                                  <  7:45804/149‑1 (MQ=255)    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                                                                                                                                  <  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                                                                           <  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                                                                                                      >  3:27134/1‑101 (MQ=255)                                                                                              cTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACa                                                                                                      <  4:27134/101‑1 (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                                                                           >  8:354415/1‑90 (MQ=255)                                                                                                                                    cGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCgg                                                                           <  7:354415/90‑1 (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                  <  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)                                                                                                                                                    |                                                                                                                                                    TCGATCACCAGTGCTTCGGCCAGGCCAAACGGCTGCGCGTAGCGCCGTGACTCGACCAGCGCA‑TCAGCAGCGAACGACCAGCAACGCTCGACCTTGCCACTGTCGATATCACGTCGGCACACCCGGTAGGCGTT‑ACGTTCGAGGGTAAACAGCTTGCCCTCGAACCAGGCCAGGTCAGCGAAGTCGCGCGACAGGGCCCTGGCATTGGGCATCTGCGCGGGTTGCATTTCGACCCCC‑GCCTCGCTCAGCAGCACGCAACCACGCCCGCAGGTCCAGACCGACTGCTGGCGCT  >  NC_002947/5586378‑5586669

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