Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,377,831	T→C	5.8%	R116G (AGA→GGA)	PP_5638 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,377,831	0	T	C	5.8%	202.4 / 2.6	69	R116G (AGA→GGA)	PP_5638	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/32);  new base C (2/2);  total (35/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.37e-01

GCCTGCTCATATAGTCGTACTCCATGGTGCACTGGTCTCCCGCTATCTTGGCTCGGTCATACGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCAGCATCAGCAATGGTTCTTTCGTCCTGTGCATGGGCCCTCGCC  >  NC_002947/4377684‑4377976                                                                                                                                                    |                                                                                                                                                  gCCTGCTCATATAGTCGTACTCCATGGTGCACTGGTCTCCCGCTATCTTGGCTCGGTCATACGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTg                                                                                                                                                  >  5:298004/1‑149 (MQ=255) accTGCTCATATAGTCGTACTCCATGGTGCACTGGTCTCCCGCTATCTTGTCTCGGTCATACGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTg                                                                                                                                                  <  8:35982/148‑1 (MQ=255)          tataGTCGTACTCCATGGTGCACTGGTCTCCCGCTATCTTGGCTCGGTCATACGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAgcg                                                                                                                                         >  8:275170/1‑149 (MQ=255)          tataGTCGTACTACATGGTGCACTGGTCTCCCGCTATCTTGGCTCGGTCATACGCCGTCGCCAACTCTCACCCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCAGGCCTCTGGCGACAgcg                                                                                                                                         >  1:64047/1‑149 (MQ=255)             aGTCGTACTCCATGGTGCACTGGTCTCCCGCTATCTTGGCTCGGTCATACGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAgcggcg                                                                                                                                      >  2:163747/1‑149 (MQ=255)              gTCGTACTCCATGGTGCACTGGTCTCCCGCTATCTTGGCTCGGTCATACGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGAc                                                                                                                                             <  4:71163/141‑1 (MQ=255)              gTCGTACTCCATGGTGCACTGGTCTCCCGCTATCTTGGCTCGGTCATACGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGAc                                                                                                                                             >  3:71163/1‑141 (MQ=255)                  tACTGCATGGTGCAGTGGTCTCGCGCTATCTTGGCTCGGTCATACGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATc                                                                                                                                 <  2:64047/149‑1 (MQ=255)                     tCCATGGTGCACTGGTCTCCCGCTATCTTGGCTCGGTCATACGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATcc                                                                                                                                >  4:49568/1‑147 (MQ=255)                       cATGGTGCACTGGTCTCCCGCTATCTTGGCTCGGTCATACGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGc                                                                                                                            >  8:234457/1‑149 (MQ=255)                           gTGCACTGGTCTCCCGCTATCTTGGCTCGGTCATACGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCggg                                                                                                                        <  6:267713/149‑1 (MQ=255)                                 tGGTCTCCCGCTATCTTGGCTCGGTCATACGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACt                                                                                                                  >  1:41830/1‑149 (MQ=255)                                      tcCCGCTATCTTGGCTCGGTCATACGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAAc                                                                                                                   <  1:55725/143‑1 (MQ=255)                                      tcCCGCTATCTTGGCTCGGTCATACGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAAc                                                                                                                   >  2:55725/1‑143 (MQ=255)                                        ccGCTATCTTGGCTCGGTCATACGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGAtggtg                                                                                                           >  3:328793/1‑149 (MQ=255)                                        ccGCTATCTTGGCTCGGTCATACGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGAtgg                                                                                                             <  1:163747/147‑1 (MQ=255)                                        ccGCTATCTTGGCTCGGTCATACGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGAtg                                                                                                              <  1:2551/146‑1 (MQ=255)                                          ccTATCTTGGCTCGGTCATACGCCGTCGCCAACTCCCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCgg                                                                                                                         <  4:9091/132‑1 (MQ=255)                                           cTATCTTGGCTCGGTCATACGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTggcg                                                                                                        <  4:84939/149‑1 (MQ=255)                                                  ggCTCGGTCATACGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCt                                                                                                 >  1:33565/1‑149 (MQ=255)                                                    cTCGGTCATACGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCCGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTgg                                                                                               >  1:40959/1‑149 (MQ=255)                                                     tCGGTCATACGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTgg                                                                                               >  5:323412/1‑148 (MQ=255)                                                       ggTCATACGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTgg                                                                                               <  6:323412/146‑1 (MQ=255)                                                       ggTCATACGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTg                                                                                            <  1:334956/149‑1 (MQ=255)                                                           aTACGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCAt                                                                                        <  5:165299/149‑1 (MQ=255)                                                             aCGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGgcggcg                                                                                                     <  8:223823/134‑1 (MQ=255)                                                             aCGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGgcggcg                                                                                                     >  7:223823/1‑134 (MQ=255)                                                               gCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTg                                                                                            >  5:229153/1‑141 (MQ=255)                                                               gCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTg                                                                                            <  6:229153/141‑1 (MQ=255)                                                                ccGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGTGGCGAGCTGGGTg                                                                                            <  2:6832/140‑1 (MQ=255)                                                                ccGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGTGGCGAGCTGGGTg                                                                                            >  1:6832/1‑140 (MQ=255)                                                                 cGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTcc                                                                                  <  1:81572/149‑1 (MQ=255)                                                                       cAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTccgcagcc                                                                            >  3:92560/1‑149 (MQ=255)                                                                       cAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTccgcagcc                                                                            >  4:46038/1‑149 (MQ=255)                                                                             tcCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCCGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGAc                                                                      >  3:300858/1‑149 (MQ=255)                                                                              cccGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGAc                                                                      >  1:210364/1‑148 (MQ=255)                                                                                cGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGAGGCGCGGGTGGCTGCCTGGCCTCTGGCGCCAGCCGCGGTATCCGAGCCGGGGCAACTGATGGTGGCCGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACaag                                                                   >  6:171785/1‑149 (MQ=21)                                                                                  cTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACcagca                                                                 <  7:234457/149‑1 (MQ=255)                                                                                     gAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGCCAGCGGCGGTATCCGCGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAgcagca                                                              >  5:87373/1‑149 (MQ=255)                                                                                       gCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCCGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATc                                                            <  2:40959/149‑1 (MQ=255)                                                                                        cATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCg                                                           <  4:92560/149‑1 (MQ=255)                                                                                         aTCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCgg                                                          >  5:253998/1‑149 (MQ=255)                                                                                            aGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTGTGGCGACAGCGGCGGTATCCGAGCGGGGGCAACTGATGGTGGCGGCGAGCTGGGTGTCATCACTCCGGAGCCGCAGACCAGCAGCATCGGCGt                                                       <  5:171785/149‑1 (MQ=255)                                                                                                 gTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTcag                                                    >  1:259789/1‑147 (MQ=255)                                                                                                       gCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCa                                            <  3:46038/149‑1 (MQ=255)                                                                                                        cAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTccgcagcc                                                                            <  6:257603/116‑1 (MQ=255)                                                                                                        cAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCAg                                           >  1:58797/1‑149 (MQ=255)                                                                                                           gTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCAGCAt                                        >  8:82462/1‑149 (MQ=255)                                                                                                           gTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCAGCAt                                        >  5:85667/1‑149 (MQ=255)                                                                                                              ggAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCAGCATCa                                      <  7:275170/148‑1 (MQ=255)                                                                                                               gagaGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAgcggc                                                                                                                                       >  2:97950/1‑50 (MQ=38)                                                                                                               gagaGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCAGCATCAg                                     >  5:237307/1‑148 (MQ=255)                                                                                                               gagaGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCCGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCAGCATCAGc                                    <  4:300858/149‑1 (MQ=255)                                                                                                                   gCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCAGCATCAGCAAt                                 <  2:210364/148‑1 (MQ=255)                                                                                                                      ccATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCAGCATCAGCAATGGtt                             <  7:82462/149‑1 (MQ=255)                                                                                                                      ccATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCAGCATCAGCAATGGtt                             <  3:49568/149‑1 (MQ=255)                                                                                                                        aTGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCAGCATCAGCAATGGTTCt                           >  3:7668/1‑149 (MQ=255)                                                                                                                                   ggTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCAGCATCAGCAAt                                 <  4:67092/132‑1 (MQ=255)                                                                                                                                    gTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCAGCATCAGCAATGGTTCTTTCGTCCtgtg                >  8:269432/1‑148 (MQ=255)                                                                                                                                       gCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCAGCATCAGCAATGGTTCTTTCGTCCTGTGCAt             <  3:271082/148‑1 (MQ=255)                                                                                                                                       gCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCAGCATCAGCAATGGTTCTTTCGTCCTGTGCAt             <  2:259789/148‑1 (MQ=255)                                                                                                                                        cTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCAGCATCAGCAATGGTTCTTTCGTCCTGTGCATgg           <  2:33565/149‑1 (MQ=255)                                                                                                                                        cTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCAGCATCAGCAATGGTTCTTTCGTCCTGTGCATgg           <  4:328793/149‑1 (MQ=255)                                                                                                                                           ccTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCAGCATCAGCAATGGTTCtt                          >  6:3800/1‑131 (MQ=255)                                                                                                                                           ccTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCAGCATCAGCAATGGTTCtt                          <  5:3800/131‑1 (MQ=255)                                                                                                                                            cTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCAGCATCAGCAATGGTTCTTTCGTCCTGTGCATGGGccc       >  1:158472/1‑149 (MQ=255)                                                                                                                                             tGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCAGCATCAGCAATGGTTCTTTCGTCCTGTGCATGGGCCCt      >  5:39995/1‑149 (MQ=255)                                                                                                                                              ggCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCAGCATCAGCAATGGTTCTTTCGTCCTGTGCATgg           >  6:221452/1‑143 (MQ=255)                                                                                                                                              ggCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCAGCATCAGCAATGGTTCTTTCGTCCTGTGCATgg           <  5:221452/143‑1 (MQ=255)                                                                                                                                                 ctctGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCAGCATCAGCAATGGTTCTTTCGTCCTGTGCATGGGCCCTCGcc  <  2:41830/149‑1 (MQ=255)                                                                                                                                                    |                                                                                                                                                  GCCTGCTCATATAGTCGTACTCCATGGTGCACTGGTCTCCCGCTATCTTGGCTCGGTCATACGCCGTCGCCAACTCTCCCGCTCGAGCATCAGCCCGTGCGAGCAGGTCGGAGAGCACCATGGCGGCGCGGGTGGCTGCCTGGCCTCTGGCGACAGCGGCGGTATCCGAGCCGGGGCAACTGATGGTGGCGGCGAGCTGGGTGGCATCACTCCGCAGCCGCAGACCAGCAGCATCGGCGTCAGCAGCGCCAGCATCAGCAATGGTTCTTTCGTCCTGTGCATGGGCCCTCGCC  >  NC_002947/4377684‑4377976

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