Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	5,706,587	T→A	6.9%	L202Q (CTG→CAG)	PP_5009 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,706,587	0	T	A	6.9%	212.8 / 5.3	72	L202Q (CTG→CAG)	PP_5009	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/33);  new base A (2/3);  total (35/37)
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.88e-01

GCTCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCATGTC  >  NC_002947/5706441‑5706732                                                                                                                                                    |                                                                                                                                                  gCTCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTgg                                                                                                                                                 <  8:286313/149‑1 (MQ=255)     gacgaAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTgg                                                                                                                                                 >  3:39925/1‑145 (MQ=255)     gacgaAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTgg                                                                                                                                                 <  4:39925/145‑1 (MQ=255)     gacgaAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCt                                                                                                                                               <  1:367607/147‑1 (MQ=255)      acgaAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTc                                                                                                                                            <  2:47874/149‑1 (MQ=255)      acgaAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCAGGCTGTc                                                                                                                                            <  3:246545/149‑1 (MQ=255)         aaGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAg                                                                                                                                         <  5:163521/149‑1 (MQ=255)          aGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAgg                                                                                                                                        >  6:153793/1‑149 (MQ=255)          aGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCAGGCTGTCgtgt                                                                                                                                        <  4:242541/149‑4 (MQ=255)                ccAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAgg                                                                                                                                        <  7:300229/143‑1 (MQ=255)                ccAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAgg                                                                                                                                        >  8:300229/1‑143 (MQ=255)                ccAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATa                                                                                                                                   >  3:93477/1‑148 (MQ=255)                  aaGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACaa                                                                                                                                <  8:343203/149‑1 (MQ=255)                           gcATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTc                                                                                                                       >  4:359107/1‑149 (MQ=255)                           gcATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAAACTCCCTCGCAAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTc                                                                                                                       >  4:179226/1‑149 (MQ=255)                              tGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCAGGCTGTCgtgt                                                                                                                                        <  5:162804/129‑4 (MQ=255)                              tGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCAGGCTGTCgtgt                                                                                                                                        >  6:162804/1‑126 (MQ=255)                                       tggAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGc                                                                                                           >  6:66283/1‑149 (MQ=255)                                       tggAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGc                                                                                                           >  8:183980/1‑149 (MQ=255)                                          aGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCggg                                                                                                        <  5:29259/149‑1 (MQ=255)                                           gCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCgggg                                                                                                       <  6:92779/149‑1 (MQ=255)                                               acCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGa                                                                                                               >  2:28275/1‑137 (MQ=255)                                               acCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGa                                                                                                               <  1:28275/137‑1 (MQ=255)                                                 cATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCggg                                                                                                        <  3:316961/142‑1 (MQ=255)                                                 cATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCggg                                                                                                        >  4:316961/1‑142 (MQ=255)                                                 cATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCg                                                                                                          <  5:208247/140‑1 (MQ=255)                                                 cATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCg                                                                                                          >  6:208247/1‑140 (MQ=255)                                                   tCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATCCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGCCGGGGCACACAcc                                                                                               <  2:222007/149‑1 (MQ=255)                                                         cTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTAcc                                                                                         <  4:93477/149‑1 (MQ=255)                                                          tGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTAccc                                                                                        >  4:112590/1‑149 (MQ=255)                                                                  tCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGa                                                                                 >  6:244906/1‑148 (MQ=255)                                                                    gagaACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATc                                                                               >  5:369182/1‑148 (MQ=255)                                                                    gagaACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCt                                                                              >  8:289750/1‑149 (MQ=255)                                                                     agaACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTa                                                                             <  7:183980/149‑1 (MQ=255)                                                                     agaACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTa                                                                             >  3:365916/1‑149 (MQ=255)                                                                      gaACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTAt                                                                            <  5:66283/149‑1 (MQ=255)                                                                       aaCCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTAtt                                                                           <  4:365916/149‑1 (MQ=255)                                                                       aaCCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCGAtt                                                                           >  8:23294/1‑149 (MQ=255)                                                                        aCCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTg                                                                          <  5:244906/149‑1 (MQ=255)                                                                            ccctcGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGACGGGGGTACCCGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCt                                                                        >  1:101850/1‑148 (MQ=255)                                                                                 gCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCtcatca                                                                   <  7:317246/147‑1 (MQ=255)                                                                                 gCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCtcatca                                                                   <  3:215136/147‑1 (MQ=255)                                                                                  ccAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATcacca                                                                >  2:26152/1‑149 (MQ=255)                                                                                        ccTGGCAGATGATGCAACAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTg                                                          <  6:369182/149‑1 (MQ=255)                                                                                        ccTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTg                                                          <  7:289750/149‑1 (MQ=255)                                                                                        catcgCAGATGAGGCAACAGCTCCAATGTACGCCCTGCTCGAC‑GGGCTCCTCGCCAGCCGGGCTGTCGAGGCGATCCACGCCACGCTCTGAACGGCTGCCGGGGCACACACCTTTACCCTACCGGATCTATTGCTCATCACCACAAGTg                                                          <  4:277528/145‑1 (MQ=255)                                                                                            gCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGACGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGgcaa                                                      >  5:216021/1‑149 (MQ=255)                                                                                                 tgatgCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGgcaaccgc                                                  <  3:179226/148‑1 (MQ=255)                                                                                                   atgCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGgcaaccgcaac                                               <  8:126164/149‑1 (MQ=255)                                                                                                     gCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCt                                             >  2:171078/1‑149 (MQ=255)                                                                                                        aaaGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGgcaaccgcagcttact                                          >  1:132789/1‑142 (MQ=255)                                                                                                        aaaGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGt                                          <  3:359107/149‑1 (MQ=255)                                                                                                             tGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCtg                                      >  4:232284/1‑148 (MQ=255)                                                                                                                    tACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCtgcaatgc                                >  4:362750/1‑147 (MQ=255)                                                                                                                    tACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGACGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCtgcaatgcaa                              <  6:216021/149‑1 (MQ=255)                                                                                                                      cGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAg                             <  5:153793/148‑1 (MQ=255)                                                                                                                       gCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGcc                           >  4:123401/1‑149 (MQ=255)                                                                                                                        cccTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGccc                          >  3:286154/1‑149 (MQ=255)                                                                                                                        cccTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGccc                          >  2:248851/1‑149 (MQ=255)                                                                                                                             cTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGtt                                         <  7:142478/129‑1 (MQ=255)                                                                                                                             cTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGtt                                         >  8:142478/1‑129 (MQ=255)                                                                                                                             cTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGcgg                     >  2:192840/1‑149 (MQ=255)                                                                                                                                gAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGa                                                                                 >  4:249610/1‑86 (MQ=255)                                                                                                                                gAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGa                                                                                 <  3:249610/86‑1 (MQ=255)                                                                                                                                 aC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGcggttcg                 <  5:122463/149‑1 (MQ=255)                                                                                                                                      ggTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCggttcggttc             <  1:26152/149‑1 (MQ=255)                                                                                                                                          cacGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTtcgtc          <  3:112590/148‑1 (MQ=255)                                                                                                                                           tcGACAGCCAGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGc        >  4:328688/2‑149 (MQ=255)                                                                                                                                           acGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGc        >  6:321198/1‑149 (MQ=255)                                                                                                                                             gACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCa      >  8:332979/1‑149 (MQ=255)                                                                                                                                               cAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCATg    <  7:23294/149‑1 (MQ=255)                                                                                                                                                  ccTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCATGTc  >  6:122742/1‑148 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GCTCGACGAAGACAGCCAAGTCTCGCGCATGCTCTGGCTGGAGCACACCATCAAGGCTGCATTGCTCGAGAACCTCCCTCGCCAAGGCCTGGCAGATGATGCAAAAGCTGCAATGTACGCCCTGCTCGAC‑GGGGTACACGACAGCCTGGCTGTCGAGGCGATACAAGTCACGCTCGGAACGGATGGCGGGGCACACACCGTTACCCTACCGGATCTATTGCTCATCACCACAAGTGGCAACCGCAACCTGGTTCTGCAATGCAAGCCCTGCGGTTCGGTTCGTCGCCATGTC  >  NC_002947/5706441‑5706732

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