Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,799,744	A→T	6.6%	T27S (ACA→TCA)	PP_3360 →	membrane protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,799,744	0	A	T	6.6%	221.4 / 4.9	76	T27S (ACA→TCA)	PP_3360	membrane protein
Reads supporting (aligned to +/- strand):  ref base A (39/32);  new base T (2/3);  total (41/35)
Fisher's exact test for biased strand distribution p-value = 6.57e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.12e-01

ATCGACCTGTTGCGCCTCTTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGC  >  NC_002947/3799596‑3799888                                                                                                                                                      |                                                                                                                                                 aTCGACCTGTTGCGCCTCTTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGa                                                                                                                                                  >  3:365643/1‑149 (MQ=255)    gACCTGTTGCGCCTCTTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGATGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACa                                                                                                                                                <  3:161960/148‑1 (MQ=255)         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)                ctctTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGa                                                                                                                                                  <  5:65500/134‑1 (MQ=255)                ctctTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGa                                                                                                                                                  >  6:65500/1‑134 (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                                                                                                                                    <  8:18466/138‑1 (MQ=255)                          tCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTg                                                                                                                                    <  8:23992/138‑1 (MQ=255)                          tCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTg                                                                                                                                    >  7:18466/1‑138 (MQ=255)                          tCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTg                                                                                                                                    >  7:23992/1‑138 (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                                                                                                                     >  1:252562/1‑143 (MQ=255)                                    gCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGcc                                                                                                                     <  2:252562/143‑1 (MQ=255)                                     cACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTgccg                                                                                                              <  2:332687/149‑1 (MQ=255)                                      aCCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGt                                                                                                                                  <  6:98616/128‑1 (MQ=255)                                      aCCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGt                                                                                                                                  >  5:98616/1‑128 (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                                                                            >  1:20281/1‑138 (MQ=255)                                                                                  catcatCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTg                                                                            <  2:20281/138‑1 (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:58346/1‑149 (MQ=255)                                                                                                                            cATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCa                       >  5:64047/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)                                                                                                                                                      |                                                                                                                                                 ATCGACCTGTTGCGCCTCTTTGACGTCCCAAGAAAGCACCTTGGCAAATCAACCCAGAGGCGAGGTACTCATGAACGACCGCATCATCGCGTTTTTCACCGAAACCACCGTATTCGGCATCTCCATTCTCAATCTGCTGGTGGCG‑CTGACAGTCGCCACCTTGGTCTTCCTGGTCGCCCGTGCCGCCATCAGTTTCCTGATGCGCCGGGTACGGCGCTGGTCCGAACACGAAGGCCCACTGAGCCAGATAATGGCCAAGGTGCTTTCCGGCACCAGCAACTTTCTGTTGTTGC  >  NC_002947/3799596‑3799888

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