Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,989,288	G→C	10.5%	G157G (GGC→GGG) ‡	PP_4398 ←	MFS transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,989,288	0	G	C	10.5%	192.2 / 12.3	67	G157G (GGC→GGG) ‡	PP_4398	MFS transporter
Reads supporting (aligned to +/- strand):  ref base G (28/32);  new base C (4/3);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 7.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.06e-01

CCGTACTGAGCGACTGCGGCACGCGCTTGATGAAGAACTTCGGCTCCATCGGTGCCGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCTTCACATCGGCCTGCTCATTGAGCCAGCTCTCGATGACCATGTATTG  >  NC_002947/4989140‑4989428                                                                                                                                                     |                                                                                                                                             ccGTACTGAGCGACTGCGGCACGCGCTTGATGAAGAACTTCGGCTCCATCGGTGCCGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAggg                                                                                                                                              <  5:55332/149‑1 (MQ=255)   gTACTGAGCGACTGCGGCACGCGCTTGATGAAGAACTTCGGCTCCATCGGTGCCGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGcc                                                                                                                                            >  7:14384/1‑149 (MQ=255)    tACTGAGCGACTGCGGCACGCGCTTGATGAAGAACTTCGGCTCCATCGGTGCCGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACAAACATCAGCAGCTCAGGGcc                                                                                                                                            >  2:290186/1‑148 (MQ=255)      cTGAGCGACTGCGGCACGCGCTTGATGAAGAACTTCGGCTCCATCGGTGCCGGGCGCAAGGGCGCCGGGGGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGa                                                                                                                                          <  6:234515/148‑1 (MQ=255)         aGCGACTGCGGCACGCGCTTGATGAAGAACTTCGGCTCCATCGGTGCCGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTg                                                                                                                                      >  7:155378/1‑149 (MQ=255)           cGACTGCGGCACGCGCTTGATGAAGAACTTCGGCTCCATCGGTGCCGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCg                                                                                                                                    >  7:238992/1‑149 (MQ=255)            gACTGCGGCACGCGCTTGATGAAGAACTTCGGCTCCATCGGTGCCGGGCGCAAGGGCGCCGGGTTAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCGGCCCTGAGCTGctg                                                                                                                                   <  7:129414/149‑3 (MQ=255)             aCTGCGGCACGCGCTTGATGAAGAACTTCGGCTCCATCGGTGCCGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCggg                                                                                                                                  >  7:336757/1‑149 (MQ=255)                 cGGCACGCGCTTGATGAAGAACTTCGGCTCCATCGGTGCCGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGca                                                                                                                              >  2:173560/1‑149 (MQ=255)                 cGGCACGCGCTTGATGAAGAACTTCGGCTCCATCGGTGCCGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGc                                                                                                                               <  8:47689/148‑1 (MQ=255)                 cGGCACGCGCTTGATGAAGAACTTCGGCTCCATCGGTGCCGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGc                                                                                                                               <  8:155378/148‑1 (MQ=255)                 cGGCACGCGCTTGATGAAGAACTTCGGCTCCATCGGTGCCGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGc                                                                                                                               <  8:238992/148‑1 (MQ=255)                     aCGCGCTTGATGAAGAACTTCGGCTCCATCGGTGCCGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGcaccac                                                                                                                          >  3:306327/1‑149 (MQ=255)                           ttGATGAAGAACTTCGGCTCCATCGGTGCCGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAAt                                                                                                                    >  6:389413/1‑149 (MQ=255)                             gatgaAGAACTTCGGCTCCATCGGTGCCGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGTCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAAt                                                                                                                    <  8:293393/147‑1 (MQ=255)                               tgaAGAACTTCGGCTCCATCGGTGCCGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAAc                                                                                                                <  7:217063/149‑1 (MQ=255)                                 aagaaCTTCGGCTCCATCGGTGCCGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGTCAAAGCACATCACCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTg                                                                                                              >  4:338903/1‑149 (MQ=255)                                 aagaaCTTCGGCTCCATCGGTGCCGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTg                                                                                                              <  5:410174/149‑1 (MQ=255)                                          ggCTCCATCGGTGCCGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCAc                                                                                                     >  4:241425/1‑149 (MQ=255)                                           gCTCCATCGGTGCCGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCAc                                                                                                     >  4:106659/1‑148 (MQ=255)                                             tCCATCGGTGCCGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTACGGGTGCACCACCAGAATCAACTGGc                                                                                                            >  7:310923/1‑139 (MQ=255)                                             tCCATCGGTGCCGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTACGGGTGCACCACCAGAATCAACTGGc                                                                                                            <  8:310923/139‑1 (MQ=255)                                                aTCGGTGCCGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGccc                                                                                               <  3:384899/149‑1 (MQ=255)                                                   ggTGCCGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAg                                                                                            <  3:106659/149‑1 (MQ=255)                                                    gTGCCGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCCGCCCCAgg                                                                                           >  8:368195/1‑149 (MQ=255)                                                       ccGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTAc                                                                                        <  4:306327/149‑1 (MQ=255)                                                          ggCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAg                                                                                     >  6:374926/1‑149 (MQ=255)                                                               aaGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCCCCACCAGAATCCACTGGCCAAGCACCAGCCCCAGGTACGAGGCAAt                                                                                >  2:272401/1‑149 (MQ=255)                                                                  ggCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCAACAGCCCCAGGTACGAGGCAATCAt                                                                             <  1:290186/149‑1 (MQ=255)                                                                     gcCGTGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTa                                                                          <  2:96251/149‑1 (MQ=255)                                                                          gTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGc                                                                     >  6:378510/1‑149 (MQ=255)                                                                              aTACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGaa                                                                 <  3:241425/149‑1 (MQ=255)                                                                                  gCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACAcc                                                             <  8:379914/149‑1 (MQ=255)                                                                                      ggTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCgcgc                                                         >  4:336806/1‑149 (MQ=255)                                                                                           tCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCgg                                                    <  8:336757/149‑1 (MQ=255)                                                                                             gCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTg                                                  <  5:378510/149‑1 (MQ=255)                                                                                               cACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGc                                                                                                 <  2:221332/100‑1 (MQ=255)                                                                                               cACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGc                                                                                                 >  1:221332/1‑100 (MQ=255)                                                                                                aCGGGCACCAGGCACAGGTCAAAGCACATCACCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCtt                                               <  3:338903/149‑1 (MQ=255)                                                                                                        cAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGa                                                                                      <  2:6993/102‑1 (MQ=255)                                                                                                        cAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGa                                                                                      >  1:6993/1‑102 (MQ=255)                                                                                                        cAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCtt                                               >  3:335969/1‑141 (MQ=255)                                                                                                        cAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCtt                                               <  4:335969/141‑1 (MQ=255)                                                                                                        cAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCAATGGACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCTTCACATCgg                                       >  7:149568/1‑149 (MQ=255)                                                                                                           gCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCTTCACATCGGCCt                                    <  8:14384/149‑1 (MQ=255)                                                                                                                gggCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCTTCACATCGGCCTGCTCa                               >  8:201072/1‑149 (MQ=255)                                                                                                                 ggCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAAc                                                                                                                <  7:98708/67‑1 (MQ=255)                                                                                                                 ggCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAAc                                                                                                                >  8:98708/1‑67 (MQ=255)                                                                                                                  gCAAAGCACATCGGCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCTTCACATCGGCCTGCTCAtt                             <  6:342723/149‑1 (MQ=255)                                                                                                                         acaTCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCTTCACATCGGCCTGCTCATTGAGc                         <  7:188346/146‑1 (MQ=255)                                                                                                                         acaTCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCTTCACATCGGCCTGCTCATTGAGc                         >  8:188346/1‑146 (MQ=255)                                                                                                                         acaTCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCTTCACATCGGCCTGCTCATTGAGCCAg                      <  7:201072/149‑1 (MQ=255)                                                                                                                                caccaACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCTTCACATCGGCCt                                    >  7:134139/1‑128 (MQ=255)                                                                                                                                caccaACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCTTCACATCGGCCt                                    <  8:134139/128‑1 (MQ=255)                                                                                                                                  ccaACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCTTCACATCGGCCTGCTCATTGAGCCAGCTCTCGATg             >  7:333166/1‑149 (MQ=255)                                                                                                                                   caACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCTTCACATCGGCCTGCTCATTGAGCCAGCTCTCGATGa            <  5:374926/149‑1 (MQ=255)                                                                                                                                     aCATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCTTCACATCGGCCTGCTCATTGAGCCAGCTCTCGATGAcc          <  4:250299/149‑1 (MQ=255)                                                                                                                                        tCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCTTCACATCGGCCTGCTCATTGAGCCAGCTCTCGATGACCAt        <  4:257890/148‑1 (MQ=255)                                                                                                                                        tCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCTTCACATCGGCCTGCTCATTGAGCCAGCTCTCGATGACCAt        >  3:257890/1‑148 (MQ=255)                                                                                                                                        cccgcagcTCGGCCCTGAGCTGCTGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCTTCACATCGGCCTGCTCATTGAGCCAGCTCTCGATGACCATg       >  1:96048/4‑149 (MQ=255)                                                                                                                                        cccgcagcTCGGCCCTGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCTTCACATCGGCCTGCTCATTGAGc                         <  1:398552/128‑1 (MQ=255)                                                                                                                                        cccgcagcTCGGCCCTGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCTTCACATCGGCCTGCTCATTGAGc                         >  2:398552/4‑131 (MQ=255)                                                                                                                                        cccgcagcTCGGCCCTGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCTTCACATCGGCCTGCTCATTGAGCCAGCtctc                 >  3:145554/4‑139 (MQ=255)                                                                                                                                        cccgcagcTCGGCCCTGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCTTCACATCGGCCTGCTCATTGAGCCAGCtctc                 <  4:145554/136‑1 (MQ=255)                                                                                                                                         ccgcagcTCGGCCCTGAGCTGCTGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCTTCACATCGGCCTGCTCATTGAGCCAGCTCTCGATGACCATg       >  2:97527/3‑148 (MQ=255)                                                                                                                                           gcagcTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCTTCACATCGGCCTGCTCATTGAGCCAGCTCTCGATGAcc          <  3:188172/143‑1 (MQ=255)                                                                                                                                           gcagcTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCTTCACATCGGCCTGCTCATTGAGCCAGCTCTCGATGAcc          >  4:188172/1‑143 (MQ=255)                                                                                                                                             agcTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCTTCACATCGGCCTGCTCATTGAGCCAGCTCTCGATGACCATGTATTg  <  5:389413/149‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                             CCGTACTGAGCGACTGCGGCACGCGCTTGATGAAGAACTTCGGCTCCATCGGTGCCGGGCGCAAGGGCGCCGGGTGAATACGCCGGGTCATCGCCACGGGCACCAGGCACAGGGCAAAGCACATCGCCACCAACATCAGCAGCTCAGGGCCGAGCTGCGGGTGCACCACCAGAATCAACTGGCCAAGCACCAGCCCCAGGTACGAGGCAATCATGTAGCCGCTGAACACCGCGCCCCGGTGCTTCACATCGGCCTGCTCATTGAGCCAGCTCTCGATGACCATGTATTG  >  NC_002947/4989140‑4989428

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