Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,114,050	T→C	6.4%	W232R (TGG→CGG)	PP_3619 →	hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,114,050	0	T	C	6.4%	163.1 / 2.4	62	W232R (TGG→CGG)	PP_3619	hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (30/28);  new base C (2/2);  total (32/30)
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.46e-01

CACCGAACCCGGCGTCAAGCATGTCGCCGAAGTGTGGGAATACCTGCGCCACTTCATGGACCACGGCCCGGACAAACTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTCCCCTACAACCTGACCGTCGCTCTCTGTTGGTACTGTGTCTGCAAGCTTTTCAATGTCC  >  NC_002947/4113902‑4114197                                                                                                                                                     |                                                                                                                                                    cACCGAACCCGGCGTCAAGCATGTCGCCGAAGTGTGGGAATACCTGCGCCACTTCATGGACCACGGCCCGGACAAACTCCCCCCGCCCGGTGCACCCAACTGGTGGCACAAGCCACTTCACGCCATCGGCCTGCCCCCGGCAGACGCCt                                                                                                                                                     >  6:89636/1‑149 (MQ=255) cACCGAACCCGGCGTCAAGCATGTCGCCGAAGTGTGGGAATACCTGCGCCACTTCATGGACCACGGCCCGGACAAACTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCt                                                                                                                                                     <  8:104191/149‑1 (MQ=32)           ggCGTCAAGCATGTCGCCGAAGTGTGGGAATACCTGCGCCACTTCATGGACCACGGCCCGGACAAACTCCCCCCGCCCGGTGAA‑CCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTAt                                                                                                                                          <  5:73123/149‑1 (MQ=17)               tctaGCATGTCGCCGATGTGTGGGTATACCTGCGCCACTTCATGGACCTCTGCCCGGACAAACTCCCCCCGCCCTGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGc                                                                                                                                        <  4:217100/145‑1 (MQ=11)                 aaGCATGTCGCCGAAGTGTGGGAATACCTGCGCCACTTCATGGACCACGGCCCGGACAAACTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGccc                                                                                                                                      <  3:347998/148‑1 (MQ=32)                   gCATGTCGCCGAAGTGTGGGAATACCTGCGCCACTTCATGGACCACGGCCCGGACAAACTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCt                                                                                                                                   >  6:170928/1‑149 (MQ=33)                   gCATGTCGCCGAAGTGTGGGAATACCTGCGCCACTTCATGGACCACGGCCCGGACAAACTCCCCCCGCCCGGTGAACCCAACTGGAGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGccacat                                                                                                                                   >  3:217100/1‑145 (MQ=17)                          gCCGAAGTGTGGGAATACCTGCGCCACTTCATGGACCACGGCCCGGACAAACTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCCAAc                                                                                                                            >  3:43283/1‑149 (MQ=17)                              aaGTGTGGGAATACCTGCGCCACTTCATGGACCACGGCCCGGACAAACTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGc                                                                                                                        <  8:11154/149‑1 (MQ=32)                              aaGTGTGGGAATACCTGCGCCACTTCATGGACCACGGCCCGGACAAACTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGc                                                                                                                        <  6:331411/149‑1 (MQ=32)                                     ggAATACCTGCGCCACTTCATGGACCACGGCCCGGACAAACTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCg                                                                                                                 >  1:375288/1‑149 (MQ=32)                                      gAATACCTGCGCCACTTCATGGACCACGGCCCGGACAAACTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCgg                                                                                                                >  7:62820/1‑149 (MQ=32)                                       aaTACCTGCGCCACTTCATGGACCACGGCCCGGACAAACTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGt                                                                                                               <  2:271156/149‑1 (MQ=32)                                               cgcCACTTCATGGACCACGGCCCGGACAAACTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCa                                                                                                       <  4:228770/149‑1 (MQ=32)                                                gcCACTTCATGGACCACGGCCCGGACAAACTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAg                                                                                                      <  2:117920/149‑1 (MQ=32)                                                gcCACTTCATGGACCACGGCCCGGACAAACTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAg                                                                                                      <  8:62820/149‑1 (MQ=32)                                                  cACTTCATGGACCACGGCCCGGACAAACTCCCCCCGCCCGGGGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCGGGAGCAATGCAggg                                                                                                    >  4:223001/1‑149 (MQ=255)                                                        aTGGACCACGGCCCGGACACACTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCaaaaa                                                                                              <  5:89636/149‑1 (MQ=17)                                                        aTGGACCACGGCCCGGACAAACTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCaaaa                                                                                               >  7:313870/1‑148 (MQ=32)                                                             ccACGGCCCGGACAAACTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAggg                                                                                                    >  8:117672/1‑138 (MQ=32)                                                             ccACGGCCCGGACAAACTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAggg                                                                                                    <  7:117672/138‑1 (MQ=32)                                                             ccACGGCCCGGACAAACTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACt                                                                                         <  3:153671/149‑1 (MQ=32)                                                              cACGGCCCGGACAAACTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTg                                                                                        <  3:386638/149‑1 (MQ=32)                                                              cACGGCCCGGACAAACTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTg                                                                                        <  8:313870/149‑1 (MQ=32)                                                                   cccGGACAAACTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGCCCCCGGCCGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAAt                                                                                   >  7:341945/1‑149 (MQ=17)                                                                       gACAAACTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGcc                                                                               <  6:7583/149‑1 (MQ=32)                                                                         cAAACTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGCCCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGt                                                                             <  2:375288/149‑1 (MQ=17)                                                                          aaaCTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGCCCCCGGCAGAATCCCGGCGCCACTATGCCCCCTGGCTAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGt                                                                             >  5:310002/1‑148 (MQ=255)                                                                          aaaCTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAg                                                                                                      >  2:147416/1‑123 (MQ=32)                                                                          aaaCTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAg                                                                                                      <  1:147416/123‑1 (MQ=32)                                                                          aaaCTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCCGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGtt                                                                            >  5:315308/1‑149 (MQ=17)                                                                          aaaCTCCCCCCGCCCGGTGAACCAAACTGGTGGCCAAAGCCACATAACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCCAGATGCCCCCTGGCGAACCGGCGACCCCCCTGAAATGCAGGGCAAAAAGAAAAGGCAATTGCCGtt                                                                            >  1:170366/1‑149 (MQ=11)                                                                            aCTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCt                                                                          >  6:364589/1‑149 (MQ=32)                                                                             cTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTg                                                                         <  5:170928/149‑1 (MQ=32)                                                                                cccccGCCCTGGTAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTgg                                                                        <  1:354189/147‑1 (MQ=12)                                                                                  cccgcccgGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCATGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAg                                                                    >  8:343484/1‑149 (MQ=18)                                                                                  cccgcccgGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTg                                                                                 <  5:260264/136‑1 (MQ=32)                                                                                  cccgcccgGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTg                                                                                 >  6:260264/1‑136 (MQ=32)                                                                                               cccAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTcccc                                                         <  6:242677/147‑1 (MQ=37)                                                                                               cccAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTcccc                                                         >  5:242677/1‑147 (MQ=37)                                                                                                       gtggCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTCCCCTACAACCTGa                                               >  6:55893/1‑149 (MQ=38)                                                                                                        tggCACAAGCCACTTCACGCCATCTGCCTGCCCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTCCCCTACAACCTGAc                                              <  8:341945/149‑1 (MQ=255)                                                                                                            acaAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTCCCCTACAACCTGACCGTc                                          >  5:354990/1‑149 (MQ=255)                                                                                                              aaGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTCCCCTACAACCTGACCGTCGc                                        >  6:381160/1‑149 (MQ=255)                                                                                                                 ccACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGGGCTGTTCCCCTACAACCTGACCGTCGctct                                     <  4:40234/149‑1 (MQ=255)                                                                                                                     ttCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTCCCCTACAACCTGACCGTCGCTCTCTGt                                 >  3:182024/1‑149 (MQ=255)                                                                                                                       cACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCTATTGACGTTCTGGGCAGTGCTGTTCCCCTACAACCTGACCGTCGAGCTCTGTTg                               >  2:331448/1‑149 (MQ=255)                                                                                                                        aCGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTCCCCTACAACCTGACCGTCGCTCTCTGTTgg                              >  5:239139/1‑149 (MQ=255)                                                                                                                        aCGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTCCCCTACAACCTGACCGTCGCTCTCTGTTg                               <  6:354990/148‑1 (MQ=255)                                                                                                                        aCGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTCCCCTACAACCTGACCGTCGCTCTCTGTTg                               <  6:239139/148‑1 (MQ=255)                                                                                                                          gCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTCCCCTACAACCTGACCGTCGCTCTCTGTTGGTa                            >  8:296669/1‑149 (MQ=255)                                                                                                                            cATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTCCCCTACAACCTGACCGTCGCTCTCTGTTGGTACt                          >  3:20536/1‑149 (MQ=255)                                                                                                                             aTCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTCCCCTACAACCTGACCGTCGCTCTCTGTTGGTACtg                         >  8:77891/1‑149 (MQ=255)                                                                                                                                tgcctgccCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGTCGAACCGGCGAGCCCGGGGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTCCCCTACAACCTGACCGTCGCTCTCTGTTGGTACTgtgt                      <  2:170366/149‑1 (MQ=255)                                                                                                                                  ccTGACCCCGGCAGAAGCCTGGCGCAACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTCCACGACAACCTGTCCGTCGCGCTCTGTTGGTACTGTGTCt                    >  6:161673/1‑149 (MQ=14)                                                                                                                                   cTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTCCCCTACAACCTGACCGTCGCTCTCTGTTGGTACTGTGTCTg                   >  4:280397/1‑149 (MQ=255)                                                                                                                                      aCCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTCCCCTACAACCTGACCGTCGCTCTCTGTTGGTACTGTGTCTGCaa                <  7:343484/149‑1 (MQ=255)                                                                                                                                        cccGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTCCCCTACAACCTGACCGTCGCTCTCTGTTGGTACTGTGTCTGCAAGc              >  5:155604/1‑149 (MQ=255)                                                                                                                                         ccGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTCCCCTACAACCTGACCGTCGCTCTCTGTTGGTACTGTGTCTGCAAGCt             >  8:298918/1‑149 (MQ=255)                                                                                                                                          cGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTCCCCTACAACCTGACCGTCGCTCTCTGTTGGTACTGTGTCTGCAAGCtt            >  5:60409/1‑149 (MQ=255)                                                                                                                                            gCAGAAGCCCTGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCTGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTCCCCTACAACCTGACCGTCGCTCTCTGTTGGTACTGTGTCTGCAAGCtttt          <  6:310002/149‑1 (MQ=255)                                                                                                                                            gCAGAAGCCCGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTCCCCTACAACCTGACCGTCGCTCTCTGTTGGTACTGTGTCTGCAAGCtttt          <  6:315308/149‑1 (MQ=255)                                                                                                                                            gCAGAAGCATGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTCCCCTACAACCTGACCGTCGCTCTCTGTTGGTACCGTGTCTGCAAGCtttt          <  3:223001/149‑1 (MQ=255)                                                                                                                                                    cTGGCGCCACTATGCCCCCTGGCGAACCGCCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTCCCCTACAACCTGACCGTCGCTCTATGTTGGTACTGTGTCTGCAAGCTTTTCAATGTcc  >  7:92699/1‑149 (MQ=255)                                                                                                                                                     |                                                                                                                                                    CACCGAACCCGGCGTCAAGCATGTCGCCGAAGTGTGGGAATACCTGCGCCACTTCATGGACCACGGCCCGGACAAACTCCCCCCGCCCGGTGAACCCAACTGGTGGCACAAGCCACTTCACGCCATCTGCCTGACCCCGGCAGAAGCCTGGCGCCACTATGCCCCCTGGCGAACCGGCGAGCCCGGTGAAATGCAGGGCAAAAAGAACTGGCAATTGCCGTTCTGGGCAGTGCTGTTCCCCTACAACCTGACCGTCGCTCTCTGTTGGTACTGTGTCTGCAAGCTTTTCAATGTCC  >  NC_002947/4113902‑4114197

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