Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_002947	5,674,753	G→C	intergenic (+115/+25)	rhlE‑II → / ← yceI	ATP‑dependent RNA helicase/polyisoprene‑binding acidic stress response factor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	5,674,753	0	G	C	97.9%	312.3 / ‑4.8	94	intergenic (+115/+25)	rhlE‑II/yceI	ATP‑dependent RNA helicase/polyisoprene‑binding acidic stress response factor
Reads supporting (aligned to +/- strand):  ref base G (0/2);  new base C (51/40);  total (51/43)
Fisher's exact test for biased strand distribution p-value = 2.01e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGC‑G‑CG‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCA  >  NC_002947/5674615‑5674892                                                                                                                                             |                                                                                                                                             gcgcTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGtt                                                                                                                                    >  6:286871/1‑149 (MQ=255)   gctgctGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGtttt                                                                                                                                  <  3:306496/149‑1 (MQ=255)       ctgACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCg                                                                                                                                       <  8:224172/140‑1 (MQ=255)       ctgACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCg                                                                                                                                       >  7:224172/1‑140 (MQ=255)       ctgACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGc                                                                                                                              <  1:113387/149‑1 (MQ=255)         gACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTg                                                                                                                            >  6:300039/1‑149 (MQ=255)         gACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAACGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTg                                                                                                                            <  5:245879/149‑1 (MQ=255)          aCCCGCAAGGGCTGATTGAAGTACCGAGGGGGCTTCTTCGCAGCCCTATCACGACACAAGGGCGCTCCTGCAGAGATTTTCTTGATCTCTACAGGAGCGTCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTGAGc                                                                                                                              <  3:307217/146‑1 (MQ=25)            ccGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTgg                                                                                                                           >  1:159117/1‑147 (MQ=255)            ccGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTgg                                                                                                                           >  1:89435/1‑147 (MQ=255)            ccGCAAGGGCTGATAGAAGTACCGAGGGGGCTCCTTCGCAGCCCTATCACGACAAAAGGCCGCTCCTGAAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGAGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCAAGGTTTTCAGCTGGc                                                                                                                          >  5:18641/1‑148 (MQ=255)             cGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTgg                                                                                                                           >  2:297388/1‑146 (MQ=255)              gCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAg                                                                                                                        >  3:161969/1‑148 (MQ=255)                    gCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTAc                                                                                                                 <  4:161969/149‑1 (MQ=255)                     cTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACt                                                                                                                <  6:425104/149‑1 (MQ=255)                        aTAGAAGTACCGAGGGGGCGGCATCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACtt                                                                                                               <  8:212103/147‑1 (MQ=37)                            aaGTACCGAGGGGGCTGCTTCGCAGGCCTATGACGACACATGGCCTCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGC‑GTCC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCtt                                                                                                         <  6:418184/149‑1 (MQ=38)                            aaGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCtt                                                                                                         >  5:234652/1‑149 (MQ=255)                            aaGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCtt                                                                                                         >  5:65601/1‑149 (MQ=255)                                 ccGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCacac                                                                                                    >  7:269331/1‑149 (MQ=255)                                 ccGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCATTGTGTCGTGATAGTGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCacac                                                                                                    >  2:278221/1‑149 (MQ=38)                                  cGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACAcc                                                                                                   >  1:382253/1‑149 (MQ=255)                                   gatagggCTGCGAAGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCt                                                                                                  >  8:343598/5‑149 (MQ=37)                                     gggggCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACATGAGCGGACTTGTGTAGTGATAGGGCCGCAAAGCGG‑CG‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTc                                                                                                <  1:297388/149‑1 (MQ=38)                                          ctgctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGacac                                                                                           <  8:269331/149‑1 (MQ=255)                                           tgctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACt                                                                                          >  2:279908/1‑149 (MQ=255)                                           tgctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCAGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACt                                                                                          >  3:243334/1‑149 (MQ=255)                                            gctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCg                                                                                               >  2:29990/1‑143 (MQ=255)                                            gctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCg                                                                                               <  1:29990/143‑1 (MQ=255)                                             ctTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGa                                                                                        <  7:375991/149‑1 (MQ=255)                                               tCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGata                                                                                      >  6:89831/1‑149 (MQ=255)                                                     cccTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGat                                                                                       <  3:400914/142‑1 (MQ=255)                                                     cccTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGat                                                                                       >  4:400914/1‑142 (MQ=255)                                                        tATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGa                                                                                              <  6:351895/132‑1 (MQ=255)                                                        tATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGa                                                                                              >  5:351895/1‑132 (MQ=255)                                                           cACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCtt                                                                                                 <  4:99264/126‑1 (MQ=255)                                                           cACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCtt                                                                                                 >  3:99264/1‑126 (MQ=255)                                                            aCGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGa                                                                                              >  5:110345/1‑128 (MQ=255)                                                            aCGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CA‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGa                                                                                              <  6:110345/128‑1 (MQ=255)                                                                   aaGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTAATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCg                                                                   >  2:258406/1‑148 (MQ=255)                                                                    aGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGaa                                                                 <  6:386389/149‑1 (MQ=255)                                                                     ggCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAg                                                                <  6:65601/149‑1 (MQ=255)                                                                     ggCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCATTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACTTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAg                                                                <  1:278221/149‑1 (MQ=255)                                                                         gCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGc                                                            <  5:300039/149‑1 (MQ=255)                                                                              tGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCt                                                                                                          <  1:202234/98‑1 (MQ=255)                                                                              tGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCt                                                                                                          >  2:202234/1‑98 (MQ=255)                                                                                cAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCAGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCgg                                                      <  4:243334/148‑1 (MQ=255)                                                                                  gagaATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGcc                                                   <  2:147412/149‑1 (MQ=255)                                                                                      aTTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTAAGCTGGCAGAAATTACTTCTGc                                                                                                           >  3:187631/1‑89 (MQ=25)                                                                                         tCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATg                                            >  5:422219/1‑149 (MQ=255)                                                                                             gATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGt                                                                        <  3:259584/117‑1 (MQ=255)                                                                                             gATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGt                                                                        >  4:259584/1‑117 (MQ=255)                                                                                             gATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCg                                         <  6:422219/148‑1 (MQ=255)                                                                                             gATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCAGGCAGAAAATACTTCAGCTTCAAACCTTCGACACAGATATCCAGGTCCAGGGTCTGc                                                                    >  4:372603/1‑121 (MQ=37)                                                                                             gATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGC‑G‑GGCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGc                                                                    <  3:372603/121‑1 (MQ=255)                                                                                                 tctACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCTTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTc                                    >  8:272980/1‑149 (MQ=255)                                                                                                 tctACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTc                                    >  3:119436/1‑149 (MQ=255)                                                                                                 tctACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTc                                    >  4:350959/1‑149 (MQ=255)                                                                                                  ctACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCtt                                                 >  7:15801/1‑135 (MQ=255)                                                                                                  ctACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCtt                                                 <  8:15801/135‑1 (MQ=255)                                                                                                   tACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTc                                    <  2:382253/147‑1 (MQ=255)                                                                                                   tACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgt                                  >  4:136230/1‑149 (MQ=255)                                                                                                    aCAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGACGTAGCACCCGGGCCTTTGATGCCGAAGTCgtt                                 <  6:234652/149‑1 (MQ=255)                                                                                                     ctgtaGCGGCCTTTTTTAGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCTTTTTTCAGCTTGCAGAAATTACTTCTGCTTCACACCTTCGAGACTGATATCCAGGTCCAGGGTCTTCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttc                                <  6:18641/145‑1 (MQ=255)                                                                                                     cAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACCGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgtt                                 >  2:8983/1‑148 (MQ=255)                                                                                                      aGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCGCCCGGGCCttt                                                <  7:226398/132‑1 (MQ=255)                                                                                                      aGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCGCCCGGGCCttt                                                >  8:226398/1‑132 (MQ=255)                                                                                                        gAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCTTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCTCCCGGGCCTTTGATGCCGAAGTCgttcagg                             <  7:272980/149‑1 (MQ=255)                                                                                                        gAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGc                                                                                                           >  3:256394/1‑71 (MQ=25)                                                                                                        gAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGc                                                                                                           <  4:256394/71‑1 (MQ=25)                                                                                                        gAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcagg                             <  5:286871/149‑1 (MQ=255)                                                                                                            ggCCTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcaggttca                         <  4:119436/149‑1 (MQ=255)                                                                                                               cTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCg                      >  3:200019/1‑149 (MQ=255)                                                                                                               cTTGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCg                      >  5:19444/1‑149 (MQ=255)                                                                                                                ttGTGTCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgt                     >  8:200740/1‑149 (MQ=255)                                                                                                                  gtgtCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtgg                   >  7:347201/1‑149 (MQ=255)                                                                                                                    gtCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtg                 >  7:334265/1‑149 (MQ=255)                                                                                                                     tCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtgg                >  4:19804/1‑149 (MQ=255)                                                                                                                     tCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtgg                >  2:206674/1‑149 (MQ=255)                                                                                                                     tCGTGATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggtg                 <  7:343598/148‑1 (MQ=255)                                                                                                                         gATAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGc            >  1:403647/1‑149 (MQ=255)                                                                                                                          aTAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGcc                                          <  4:36158/118‑1 (MQ=255)                                                                                                                          aTAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGcc                                          >  3:36158/1‑118 (MQ=255)                                                                                                                           tAGGGCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGt          >  6:30842/1‑149 (MQ=255)                                                                                                                             gggCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGtt         <  5:89831/148‑1 (MQ=255)                                                                                                                              ggCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGt                                                               >  5:95495/1‑93 (MQ=255)                                                                                                                              ggCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGt                                                               <  6:95495/93‑1 (MQ=255)                                                                                                                              ggCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTg        >  3:137132/1‑148 (MQ=255)                                                                                                                               gCCGCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGaa      <  3:258298/149‑1 (MQ=255)                                                                                                                                  gCAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGc    >  5:154568/1‑148 (MQ=255)                                                                                                                                   cAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcaggttcag                        <  7:389804/127‑1 (MQ=255)                                                                                                                                   cAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCgttcaggttcag                        >  8:389804/1‑127 (MQ=255)                                                                                                                                   cAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGc    >  3:401918/1‑147 (MQ=255)                                                                                                                                   cAAAGCGG‑CC‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCa  <  4:401918/149‑1 (MQ=255)                                                                                                                                             |                                                                                                                                             GCGCTGCTGACCCGCAAGGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGC‑G‑CG‑CCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCA  >  NC_002947/5674615‑5674892

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