Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,922,701	T→C	5.6%	K26E (AAG→GAG)	PP_4331 ←	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,922,701	0	T	C	5.6%	209.3 / 2.0	71	K26E (AAG→GAG)	PP_4331	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (29/38);  new base C (2/2);  total (31/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.41e-01

ACCGCCGCACGCAGCTCATGGATGGCCTCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCAGCTCGGACCACTCTTCCTCGGTCATCATCTTGTCCA  >  NC_002947/4922556‑4922827                                                                                                                                                   |                                                                                                                               aCCGCCGCACGCAGCTCATGGATGGCCTCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCt                                                                                                                             >  5:341018/1‑149 (MQ=255)       gcacgcaGCTCATGGATGGCCTCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTg                                                                                                                               <  3:13841/141‑1 (MQ=255)       gcacgcaGCTCATGGATGGCCTCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCATTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTg                                                                                                                               >  4:13841/1‑141 (MQ=255)         gagcaGCTCATGGATGGCCTCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCGGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGa                                                                                                                     <  4:173847/147‑1 (MQ=255)         acgcaGCTCATGGATGGCCTCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGa                                                                                                                     <  7:74801/149‑1 (MQ=255)         acgcaGCTCATGGATGGCCTCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGa                                                                                                                     <  1:281362/149‑1 (MQ=255)         acgcaGCTCATGGATGGCCTCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGa                                                                                                                     <  2:129066/149‑1 (MQ=255)         acgcaGCTCATGGATGGCCTCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGa                                                                                                                     <  2:134101/149‑1 (MQ=255)           gcaGCTCATGGATGGTCTCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTg                                                                                                                      >  4:261661/1‑146 (MQ=255)           gcaGCTCATGGATGGTCTCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTg                                                                                                                      <  3:261661/146‑1 (MQ=255)             aGCTCATGGATGGCCTCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGaa                                                                                                                 <  3:126190/149‑1 (MQ=255)                tCATGGATGGCCTCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGaaaaa                                                                                                              >  3:135444/1‑149 (MQ=255)                   tggatggCCTCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACca                                                                                                           <  2:6748/149‑1 (MQ=255)                     gatggCCTCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACcaca                                                                                                         <  6:18076/149‑1 (MQ=255)                           cTCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCAGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGa                                                                                                                     >  8:71155/1‑131 (MQ=255)                           cTCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCAGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGa                                                                                                                     <  7:71155/131‑1 (MQ=255)                            tCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAg                                                                                                  >  2:148617/1‑149 (MQ=255)                                          cGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAgcg                                                                                    >  6:128895/1‑149 (MQ=255)                                            tGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAgcgc                                                                                   >  4:275750/1‑148 (MQ=255)                                                gTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAgg                                                                               >  7:284929/1‑148 (MQ=255)                                                    tGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGAt                                                                          >  5:302337/1‑149 (MQ=255)                                                      gTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAg                                                                                           <  6:193273/130‑1 (MQ=255)                                                      gTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAg                                                                                           >  5:193273/1‑130 (MQ=255)                                                      gTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGa                                                                           <  5:128895/146‑1 (MQ=255)                                                       tAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCACAGCGCCCAGGCCAGCGCCAGGAAGATGAc                                                                       <  1:29864/149‑1 (MQ=255)                                                        aGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGa                                                                        >  4:36860/1‑147 (MQ=255)                                                        aGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACa                                                                      >  8:1176/1‑149 (MQ=255)                                                         gTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACa                                                                      <  3:36860/148‑1 (MQ=255)                                                         gTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGGTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACa                                                                      <  6:96117/148‑1 (MQ=255)                                                            aTCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCaa                                                                  >  2:87874/1‑149 (MQ=255)                                                            aTCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCaa                                                                  >  8:148250/1‑149 (MQ=255)                                                             tCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAAc                                                                 <  1:87874/149‑1 (MQ=255)                                                                  aCGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTa                                                            <  3:275750/149‑1 (MQ=255)                                                                   cGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAg                                                           <  8:284929/149‑1 (MQ=255)                                                                         gCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAAcc                                                                <  3:243238/138‑1 (MQ=255)                                                                         gCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAAcc                                                                >  4:243238/1‑138 (MQ=255)                                                                          ccAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAAc                                                    >  7:205474/1‑149 (MQ=255)                                                                            aaGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGc                                                                                                      >  1:264073/1‑97 (MQ=255)                                                                            aaGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGc                                                                                                      <  2:264073/97‑1 (MQ=255)                                                                               tCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTa                                                            >  7:133810/1‑136 (MQ=255)                                                                               tCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTa                                                            <  8:133810/136‑1 (MQ=255)                                                                               tCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGa                                               <  4:135444/149‑1 (MQ=255)                                                                                 cTTGATGCGCTGGTCACGCAGCGCATCGCCGTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACa                                                                      >  1:29791/1‑124 (MQ=255)                                                                                 cTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACa                                                                      <  2:29791/124‑1 (MQ=255)                                                                                    gATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGtt                                           <  7:148250/148‑1 (MQ=255)                                                                                    gATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTtc                                          >  7:300471/1‑149 (MQ=255)                                                                                     aTGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTtct                                         <  8:205474/149‑1 (MQ=255)                                                                                        cgcTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTtctc                                        >  1:21608/1‑147 (MQ=255)                                                                                        cgcTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCa                                      <  6:341018/149‑1 (MQ=255)                                                                                         gcTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCTGCGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCa                                                                                            >  1:73176/1‑95 (MQ=255)                                                                                         gcTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCa                                                                                            <  2:73176/94‑1 (MQ=255)                                                                                            ggTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTCGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGAt                                              <  3:66327/137‑1 (MQ=255)                                                                                            ggTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTCGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGAt                                              <  3:71020/137‑1 (MQ=255)                                                                                            ggTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTCGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGAt                                              >  4:71020/1‑137 (MQ=255)                                                                                            ggTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTCGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGAt                                              >  4:66327/1‑137 (MQ=255)                                                                                                aCGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCAGCTCGGa                               <  2:320490/148‑1 (MQ=255)                                                                                                aCGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCAGCTCGGa                               <  7:1176/148‑1 (MQ=255)                                                                                                    aGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCAGCTCGGACCAct                          >  3:153500/1‑149 (MQ=255)                                                                                                         aTCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCAGCTCGGACCACTCTTcc                     >  6:94033/1‑149 (MQ=255)                                                                                                         aTCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCAGCTCGGACCACTCTTcc                     <  5:94033/149‑1 (MQ=255)                                                                                                             cccTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCAGCTCGGACCACTCTTCCTCgg                 <  6:302337/149‑1 (MQ=255)                                                                                                                tCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGAt                                              >  7:265797/1‑117 (MQ=255)                                                                                                                tCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGAt                                              <  8:265797/117‑1 (MQ=255)                                                                                                                     ctgctgCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCAGCTCGGACCACTCTTCCTCGGTCatcatc         <  2:21608/149‑1 (MQ=255)                                                                                                                           cgcAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCAGCTCGGACCACTCTTCCTCGGTCATCATCTTGTc    <  8:300471/148‑1 (MQ=255)                                                                                                                                ccagcT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAAc                                                    <  5:340759/95‑1 (MQ=255)                                                                                                                                ccagcT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAAc                                                    >  6:340759/1‑95 (MQ=255)                                                                                                                                        gcgcTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCAGCTCGGACCACTCTTCCTCGGTCATCATCTTGTCCa  >  4:14950/1‑138 (MQ=255)                                                                                                                                        gcgcTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCAGCTCGGACCACTCTTCCTCGGTCATCATCTTGTCCa  <  3:14950/138‑1 (MQ=255)                                                                                                                                        gcacTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCAGCTCGGACCACTCTTCCTCGGTCATCATCTTGTCCa  >  4:23096/4‑138 (MQ=255)                                                                                                                                        gcacTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCAGCTCGGACCACTCTTCCTCGGTCATCATCTTGTCCa  <  3:23096/135‑1 (MQ=255)                                                                                                                                                   |                                                                                                                               ACCGCCGCACGCAGCTCATGGATGGCCTCACCCATGCGCACCGTGCCGTTCTGGTAGTCATCCAGACGCTTGGCCAAGTCCTTGATGCGCTGGTCACGCAGCGCATCGCCCTCGGCCTGCTGCGCAGCCAGCT‑CGCGCTGGCGCTTGCTGTAGTTGAGGAAAAACCACAGGCTCAGCGCCCAGGCCAGCGCCAGGAAGATGACAGCAACCTCTAGGATCAACTCAGATGTTCTCCAGCTCGGACCACTCTTCCTCGGTCATCATCTTGTCCA  >  NC_002947/4922556‑4922827

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑