Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,367,360	G→C	9.4%	L222L (CTG→CTC)	PP_2964 →	transposase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,367,360	0	G	C	9.4%	194.1 / 7.8	66	L222L (CTG→CTC)	PP_2964	transposase
Reads supporting (aligned to +/- strand):  ref base G (27/31);  new base C (3/3);  total (31/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AAGTTACCTTCAACCTAGTCGTGACTCGAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAGGACTGCTCGCTGCCG  >  NC_002947/3367215‑3367507                                                                                                                                                  |                                                                                                                                                    aaGTTACCTTCAACCTAGTCGTGACTCGAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAg                                                                                                                                                  <  6:60411/149‑1 (MQ=255) aaGTTACCTTCAACCTAGTCGTGACTCGAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAg                                                                                                                                                  <  6:157627/149‑1 (MQ=255) aaGTTACCTTCAACCTAGTCGTGACTCGAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAg                                                                                                                                                  <  5:354936/149‑1 (MQ=255)      aCCTTCAACCTAGTCGTGACTCGAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGc                                                                                                                                              <  1:77453/148‑1 (MQ=255)      aCCTTCAACCTAGTCGTGACTCGAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCAACACGCAAGCACCTGGAGGTGCt                                                                                                                                             >  7:121580/1‑149 (MQ=255)        cTTCAACCTAGTCGTGACTCGAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGa                                                                                                                                           >  6:220263/1‑149 (MQ=255)           cAACCTAGTCGTGACTCGAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTCCAGGTGCttgc                                                                                                                                          <  3:1769/147‑4 (MQ=255)            aaCCTAGTCGTGACTCGAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGa                                                                                                                                       <  8:330099/149‑1 (MQ=255)            aaCCTAGTCGTGACTCGAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGa                                                                                                                                       <  4:261436/149‑1 (MQ=255)             aCCTAGTCGTGACTCGAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGAc                                                                                                                                      <  1:169921/149‑1 (MQ=255)               cTAGTCGTGACTCGAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAgg                                                                                                                                                 >  2:34799/1‑136 (MQ=255)                   tCGTGACTCGAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTAGAGGTGCTGATCGACCACTTg                                                                                                                                >  1:80635/1‑149 (MQ=255)                     gTGACTCGAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGc                                                                                                                              >  6:332648/1‑149 (MQ=255)                        aCTCGAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTCCAGGTGCttg                                                                                                                                           >  8:76948/1‑131 (MQ=255)                        aCTCGAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTCCAGGTGCttg                                                                                                                                           <  7:76948/133‑3 (MQ=255)                                gcCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTagag                                                                                                                    <  7:149007/148‑1 (MQ=255)                                  cAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACATGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGAc                                                                                                                                      <  8:195808/128‑1 (MQ=255)                                  cAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACATGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGAc                                                                                                                                      >  7:195808/1‑128 (MQ=255)                                       gAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTTGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCg                                                                                                             <  2:80635/148‑1 (MQ=255)                                       gAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCgg                                                                                                            >  4:53915/1‑149 (MQ=255)                                       gAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCg                                                                                                             <  4:119605/148‑1 (MQ=255)                                        aGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCggg                                                                                                           <  1:159131/149‑1 (MQ=255)                                        aGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCggg                                                                                                           <  8:224007/149‑1 (MQ=255)                                                tACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTg                                                                                                     >  5:32277/1‑147 (MQ=255)                                                  cAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAg                                                                                                  <  5:301751/148‑1 (MQ=255)                                                         aaGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCgg                                                                                                            <  3:153293/131‑1 (MQ=255)                                                         aaGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCgg                                                                                                            >  4:153293/1‑131 (MQ=255)                                                           ggAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCttt                                                                                         >  5:91929/1‑148 (MQ=255)                                                           ggAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTg                                                                                        >  7:294407/1‑149 (MQ=255)                                                           ggAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTg                                                                                        <  3:53915/149‑1 (MQ=255)                                                             aGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAg                                                                                      <  7:211597/149‑1 (MQ=255)                                                             aGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAg                                                                                      <  2:315649/149‑1 (MQ=255)                                                                tgctCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAggg                                                                                    >  4:274187/1‑148 (MQ=255)                                                                          cGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTg                                                                                        <  2:188587/134‑1 (MQ=255)                                                                          cGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTg                                                                                        >  1:188587/1‑134 (MQ=255)                                                                          cGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATg                                                                         >  7:41208/1‑149 (MQ=255)                                                                            aCTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCGGGATAGAGAGCTACGGGGCCGTGGAGCCCATCTTTGAGGGGATCGACGATGc                                                                        >  6:321459/1‑148 (MQ=255)                                                                              ttCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGcc                                                                       >  8:196016/1‑147 (MQ=255)                                                                                cggccggcGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCaa                                                                     <  5:220263/147‑1 (MQ=255)                                                                                cggccggcGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGAGCGTCGATGcc                                                                       >  8:225876/1‑145 (MQ=255)                                                                                  gccggcGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGAtt                                                                 <  8:294407/149‑1 (MQ=255)                                                                                     ggcGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGc                                                              <  5:332648/149‑1 (MQ=255)                                                                                           agagCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACttt                                                        >  7:253202/1‑149 (MQ=255)                                                                                                          cgcTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCcgcg                                         <  8:41208/149‑1 (MQ=255)                                                                                                             tGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCgcggcg                                      >  2:280864/1‑149 (MQ=255)                                                                                                              gAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCgcggcgg                                     >  5:363606/1‑149 (MQ=255)                                                                                                               aaaCAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCT‑‑AGAGCTTCGGGGACCTGGAg                                                                                                  <  8:165282/85‑1 (MQ=255)                                                                                                               aaaCAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCT‑‑AGAGCTTCGGGGACCTGGAg                                                                                                  >  7:165282/1‑85 (MQ=255)                                                                                                                        cccAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGtctc                           >  7:222796/1‑149 (MQ=255)                                                                                                                         ccAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTcc                          <  6:91929/149‑1 (MQ=255)                                                                                                                          cAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCg                         <  7:196016/149‑1 (MQ=255)                                                                                                                           aagaagCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCAAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTcc                          >  5:60786/1‑147 (MQ=255)                                                                                                                             gaagCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAg                       <  7:240672/148‑1 (MQ=255)                                                                                                                             gaagCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAg                       >  8:240672/1‑148 (MQ=255)                                                                                                                               agCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCt                     <  6:363606/148‑1 (MQ=255)                                                                                                                                 ccATCTCGCAACCACCTGGCCGGGCGGTTCGCCCCCTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGaa                  <  5:321459/149‑1 (MQ=255)                                                                                                                                  cATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGTAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAg                 <  6:32277/149‑1 (MQ=255)                                                                                                                                  cATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAg                 <  3:274187/149‑1 (MQ=255)                                                                                                                                     ctcGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAGGAc              >  4:89424/1‑149 (MQ=255)                                                                                                                                      tcGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAGGACt             >  8:2342/1‑149 (MQ=255)                                                                                                                                      tcGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAGGACt             >  3:78576/1‑149 (MQ=255)                                                                                                                                      tcGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAGGACt             >  5:219146/1‑149 (MQ=255)                                                                                                                                         tcaGCACCTCCAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAg                 <  8:210171/140‑1 (MQ=255)                                                                                                                                         tcaGCACCTCCAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAg                 >  7:210171/3‑142 (MQ=255)                                                                                                                                         tcaGCACCTCCAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAGGACTGCt          >  6:354316/3‑149 (MQ=255)                                                                                                                                         cAAGCACGTGGAGGTGCTGATAGACGACGTGGCCTGTCTAGAGAGCTTTGGGGACCTGGAGGCCATCTTTGAGGGGATCTTCGATGCAAAGATTCGCCACTTTGCTGCTCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAGGACTGCt          <  8:121580/149‑1 (MQ=255)                                                                                                                                                  ggAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAGGACTGCTCGCTGCCg  >  5:207185/1‑148 (MQ=255)                                                                                                                                                  |                                                                                                                                                    AAGTTACCTTCAACCTAGTCGTGACTCGAGCGCCAATCGAGGTGATCTACAAGCTGAAGGAGCTGCTCGACACCGACTTCGGCCGGCGCCAGAGCGATTTCAATGCGCTGAAACAGGCACCCAAGAAGCCATCTCGCAAGCACCTGGAGGTGCTGATCGACCACTTGGCCTGGCTAGAGAGCTTCGGGGACCTGGAGGCCATCTTTGAGGGGATCGTCGATGCCAAGATTCGCCACTTTGCTGCCCAAGCCGCGGCGGCGGACGTCTCCGAGCTGAAGGACTGCTCGCTGCCG  >  NC_002947/3367215‑3367507

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