Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,799,748	Δ1 bp	6.5%	coding (83/1125 nt)	PP_3360 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,799,748	0	T	.	6.5%	237.2 / 10.4	77	coding (83/1125 nt)	PP_3360	membrane protein
Reads supporting (aligned to +/- strand):  ref base T (40/32);  new base . (2/3);  total (42/35)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.94e-01

GTTGCGCCTCTTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCC  >  NC_002947/3799604‑3799896                                                                                                                                                  |                                                                                                                                                     gTTGCGCCTCTTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCGAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGc                                                                                                                                                   >  2:232207/1‑148 (MQ=255)     cgcCTCTTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCt                                                                                                                                              >  1:13870/1‑149 (MQ=255)      gcCTCTTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCGACTGTCAG‑CGCCAccag                                                                                                                                             <  5:2950/149‑3 (MQ=255)      gcCTCTTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCGACTGTCAG‑CGCCAccag                                                                                                                                             <  4:333398/149‑3 (MQ=255)         tctTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGt                                                                                                                                          >  6:337580/1‑149 (MQ=255)             tGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTc                                                                                                                                      >  1:332687/1‑149 (MQ=255)             tGACGTCCAAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAAACTGCTGGTGCCG‑CTGACAGTCGCCACCTTGGACTTc                                                                                                                                      >  8:78036/1‑149 (MQ=255)               aCGTCCGAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCTTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTcc                                                                                                                                     <  5:178068/148‑1 (MQ=255)                  tCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTg                                                                                                                                            >  7:23992/1‑138 (MQ=255)                  tCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTg                                                                                                                                            <  8:18466/138‑1 (MQ=255)                  tCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTg                                                                                                                                            >  7:18466/1‑138 (MQ=255)                  tCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTg                                                                                                                                            <  8:23992/138‑1 (MQ=255)                    ccAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCg                                                                                                                               >  2:63570/1‑149 (MQ=255)                     cAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGc                                                                                                                              >  1:249999/1‑149 (MQ=255)                     cAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCAGAATAGCGTTTTTAACCCAAACCACCGTTTTAGGCCTAAGCATTATAAATAAGCTTGTTGCG‑CAGACAGTAGGCACCTTGGTGTTACTTGTAGc                                                                                                                              >  1:293611/1‑149 (MQ=255)                       aGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGccc                                                                                                                            <  1:63570/149‑1 (MQ=255)                          aaGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTg                                                                                                                         >  2:192644/1‑149 (MQ=255)                           aGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTgc                                                                                                                        >  5:301002/1‑149 (MQ=255)                            gCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGcc                                                                                                                             <  2:252562/143‑1 (MQ=255)                            gCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGcc                                                                                                                             >  1:252562/1‑143 (MQ=255)                             cACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTgccg                                                                                                                      <  2:332687/149‑1 (MQ=255)                              aCCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGt                                                                                                                                          >  5:98616/1‑128 (MQ=255)                              aCCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGt                                                                                                                                          <  6:98616/128‑1 (MQ=255)                                cTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCa                                                                                                                   >  8:294425/1‑149 (MQ=255)                                   ggCAAATCAACCCAGAGGCGAGGTACTCTTGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCAt                                                                                                                  >  6:145955/1‑147 (MQ=255)                                   ggCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCAt                                                                                                                  >  2:155446/1‑147 (MQ=255)                                   ggCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCAt                                                                                                                  <  2:13870/147‑1 (MQ=255)                                     cAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGGTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGt                                                                                                              <  7:67089/149‑1 (MQ=255)                                      aaaTCAACCCAGAGGCGAGGGACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGtt                                                                                                             <  7:78036/149‑1 (MQ=255)                                        aTCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTc                                                                                                           <  6:28107/149‑1 (MQ=255)                                            aCCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGa                                                                                                       >  6:336053/1‑149 (MQ=255)                                             cccAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGAt                                                                                                      >  5:363286/1‑149 (MQ=255)                                                  agGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGcc                                                                                                 <  6:301002/149‑1 (MQ=255)                                                    gCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCgg                                                                                               <  5:145955/149‑1 (MQ=255)                                                    gCGACGTACTGATGAACGACCGCATGATCTCGTTTTTCACCGATAGCACCGTATTCTGCATCTGTATTCTCAATGTTCTGGTGTAG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCTCCATCAGTTTCCTGATGCGCCgt                                                                                               <  4:227712/149‑2 (MQ=255)                                                            cTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGgcg                                                                                       >  4:22229/1‑149 (MQ=255)                                                            cTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCAGCAGTTTCCTGATGCGCCGGGTACGgc                                                                                        <  2:249999/148‑1 (MQ=255)                                                             tCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGgcgc                                                                                      >  5:117433/1‑149 (MQ=255)                                                               aTGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGgcgc                                                                                      <  1:192644/147‑1 (MQ=255)                                                               aTGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTg                                                                                    <  1:232207/149‑1 (MQ=255)                                                                      aCCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGaa                                                                             <  3:22229/149‑1 (MQ=255)                                                                          catcatCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTg                                                                                    <  2:20281/138‑1 (MQ=255)                                                                          catcatCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTg                                                                                    >  1:20281/1‑138 (MQ=255)                                                                              atcGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAgg                                                                     >  1:8544/1‑149 (MQ=255)                                                                               tcGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGc                                                                    >  1:280444/1‑149 (MQ=255)                                                                               tcGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGc                                                                    >  2:274761/1‑149 (MQ=255)                                                                                 gcgTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGccc                                                                  <  1:155446/149‑1 (MQ=255)                                                                                    tttttCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACt                                                               <  2:293611/149‑1 (MQ=255)                                                                                          aCCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCa                                                         <  6:117433/149‑1 (MQ=255)                                                                                           ccGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGc                                                           >  4:137025/1‑146 (MQ=255)                                                                                              aaaCCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATa                                                     <  4:365643/149‑1 (MQ=255)                                                                                                  caccGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATgg                                                 >  7:18088/1‑149 (MQ=255)                                                                                                   accGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGc                                                >  8:87598/1‑149 (MQ=255)                                                                                                   accGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGc                                                >  5:144654/1‑149 (MQ=255)                                                                                                     cGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGcc                                               >  4:176924/1‑148 (MQ=255)                                                                                                         ttCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAac                                                                            <  1:180533/115‑1 (MQ=255)                                                                                                         ttCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAac                                                                            >  2:180533/1‑115 (MQ=255)                                                                                                            ggCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGc                                        <  5:336053/148‑1 (MQ=255)                                                                                                            cgCATCTCCATTCTCATTCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGaa                                                                             <  8:168079/110‑1 (MQ=255)                                                                                                              cATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCttt                                     >  5:263903/1‑149 (MQ=255)                                                                                                                  tcCATTCTCTATCTGGTGGTGGCG‑CTGACAGTCGCCAGCTTTGTCTTCCTGGTCTCCCGTGCCGACATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCgg                                 <  8:297214/149‑1 (MQ=255)                                                                                                                   ccATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCgg                                 >  7:122151/1‑148 (MQ=255)                                                                                                                    cATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAATGTGCTTTCCGGc                                <  2:280444/148‑1 (MQ=255)                                                                                                                    cATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCa                               >  5:64047/1‑149 (MQ=255)                                                                                                                    cATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCa                               >  5:58346/1‑149 (MQ=255)                                                                                                                               tgctgGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTcc                                   <  1:60394/134‑1 (MQ=255)                                                                                                                               tgctgGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTcc                                   >  2:60394/1‑134 (MQ=255)                                                                                                                               tgctgGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACttt                    >  1:194854/1‑149 (MQ=255)                                                                                                                                 ctgGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCt                  >  2:59031/1‑149 (MQ=255)                                                                                                                                 aaggtggCGACTGTCAG‑CGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCt                  >  2:119856/3‑149 (MQ=255)                                                                                                                                 aaggtggCGACTGTCAG‑CGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCt                  <  1:119856/147‑1 (MQ=255)                                                                                                                                  aggtggCGACTGTCAG‑CGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCtg                 >  5:227148/2‑149 (MQ=255)                                                                                                                                    gtggCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCtgtt               <  6:263903/149‑1 (MQ=255)                                                                                                                                    gtggCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCtg                 <  7:294425/147‑1 (MQ=255)                                                                                                                                         g‑cTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGc          >  1:119083/1‑149 (MQ=255)                                                                                                                                              aCAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCAGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCttt                                     >  6:349741/1‑118 (MQ=255)                                                                                                                                              aCAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCAGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCttt                                     <  5:349741/118‑1 (MQ=255)                                                                                                                                                  tCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGcctcc  <  2:8544/149‑1 (MQ=255)                                                                                                                                                  tCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGcctcc  >  1:331094/1‑149 (MQ=255)                                                                                                                                                  |                                                                                                                                                     GTTGCGCCTCTTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGCTGGCCTCC  >  NC_002947/3799604‑3799896

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