Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,249,690	A→G	6.2%	D196G (GAC→GGC)	bvgA →	virulence factors transcription regulator BvgA

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,249,690	0	A	G	6.2%	189.6 / 2.4	64	D196G (GAC→GGC)	bvgA	virulence factors transcription regulator BvgA
Reads supporting (aligned to +/- strand):  ref base A (28/32);  new base G (2/2);  total (30/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.65e-01

CTTTCTGCGCGAGAAATGACGGTGCTGCAACAACTGGCGCGAGGGCGTTCGAACAAGGAGATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑CCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCT  >  NC_002947/1249542‑1249829                                                                                                                                                     |                                                                                                                                             cTTTCTGCGCGAGAAATGACGGTGCTGCAACAACTGGCGCGAGGGCGTTCGAACAAGGAGATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGa                                                                                                                                              <  6:314617/149‑1 (MQ=255)        cgcgAGAAATGACGGTGCTGCAACAACTGGCGCGAGGGCGTTCGAACAAGGAGATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGa                                                                                                                                         <  7:10879/147‑1 (MQ=255)         gcgAGAAATGACGGTGCTGCAACAACTGGCGCGAGGGCGTTCGAACAAGGAGATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATc                                                                                                                                       >  4:187135/1‑148 (MQ=255)          cgAGAAATGACGGTGCTGCAACAACTGGCGCGAGGGCGTTCGAACAAGGAGATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATc                                                                                                                                       >  3:293479/1‑147 (MQ=255)            agaAATGACGGTGCTGCAACAACTGGCGCGAGGGCGTTCGAACAAGGAGATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGc                                                                                                                                   <  5:204719/149‑1 (MQ=255)               aaTGACGGTGCTGCAACAACTGGCGCGAGGTCGTTCGAACAAGGAGATTTCCGACAGCATGTACCTCAGTAACAAGACGGTCTGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTg                                                                                                                                 <  1:154882/148‑1 (MQ=255)                aTGACGGTGCTGCAACAACTGGTGCGAGGGCGTTCGAACAAGGAGATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGgc                                                                                                                               >  4:367557/1‑149 (MQ=255)                aTGACGGTGCTGCAACAACTGGCGCGAGGGCGTTCGAACAAGGAGATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGgc                                                                                                                               >  5:28059/1‑149 (MQ=255)                aTGACGGTGCTGCAACAACTGGCGCGAGGGCGTTCGAACAAGGAGATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGgc                                                                                                                               >  7:20558/1‑149 (MQ=255)                    cGGTGCTGCAACAACTGGCGCGAGGGCGTTCGAACAAGGAGATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCaa                                                                                                                           >  6:353535/1‑149 (MQ=255)                     ggTGCTGCAACAACTGGCGCGAGGGCGTTCGAACAAGGAGATTGCCGACAGCAATTACCTCAGTTTCAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAAc                                                                                                                          <  1:307235/149‑1 (MQ=255)                           gcAACAACTGGCGCGAGGGCGTTCGAACAAGGAGATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAg                                                                                                                    >  1:247567/1‑149 (MQ=255)                           gcAACAACTGGCGCGAGGGCGTTCGAACAAGGAGATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAg                                                                                                                    >  1:242261/1‑149 (MQ=255)                                  cTGGCGCGAGGGCGTTCGAACAAGGAGATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGaa                                                                                                             <  5:353535/149‑1 (MQ=255)                                  cTGGCGCGAGGGCGTTCGAACAAGGAGATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGaa                                                                                                             <  4:293479/149‑1 (MQ=255)                                  cTGGCGCGAGCGCGGTCGAACCAGGATAAATCCGACAGAATGTACCTCACTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGaa                                                                                                             <  8:241225/149‑1 (MQ=255)                                   tGGCTCGAGGTCGTTCGAACAAGGAGATGGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAATCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCTACGGCAAGGTCTGAAc                                                                                                            <  3:218424/149‑1 (MQ=255)                                             gCGTTCGAACAAGGAGATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCaa                                                                                                  <  1:119780/149‑1 (MQ=255)                                                   gAACAAGGAGATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCtt                                                                                             >  2:264939/1‑148 (MQ=255)                                                    aaCAAGGAGATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTg                                                                                            >  1:127770/1‑148 (MQ=255)                                                    aaCAAGGAGATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGa                                                                                           <  2:204957/149‑1 (MQ=255)                                                     aCAAGGAGATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGaa                                                                                          <  7:106611/149‑1 (MQ=255)                                                     aCAAGGAGATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGaa                                                                                          <  3:153593/149‑1 (MQ=255)                                                         ggAGATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCt                                                                                      <  2:247567/149‑1 (MQ=255)                                                         ggAGATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCt                                                                                      <  2:242261/149‑1 (MQ=255)                                                           agaTTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTc                                                                                     <  8:20558/148‑1 (MQ=255)                                                            gaTTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTcca                                                                                   >  8:92875/1‑149 (MQ=255)                                                            gaTTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTc                                                                                     >  8:160006/1‑147 (MQ=255)                                                                 ccGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTccagccag                                                                              <  7:92875/149‑1 (MQ=255)                                                                 ccGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTccagccag                                                                              >  1:296305/1‑149 (MQ=255)                                                                  cGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCagg                                                                             <  5:145438/149‑1 (MQ=255)                                                                       gCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCt                                                                        <  2:127770/149‑1 (MQ=255)                                                                            tACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTTAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCTGGCTGGAGGCttcaa                                                                    <  2:263441/148‑4 (MQ=255)                                                                             aCCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGGCCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑c                                                                   >  1:373835/1‑147 (MQ=255)                                                                               cTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTc                                                                                                                 >  6:339255/1‑100 (MQ=255)                                                                               cTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTc                                                                                                                 <  5:339255/100‑1 (MQ=255)                                                                                  aGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑CCAggg                                                              <  3:214780/147‑1 (MQ=255)                                                                                        aaGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGGCCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑CCAGGGCTATCGGt                                                      <  2:373835/149‑1 (MQ=255)                                                                                        aaGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑CCAGGGCTATCGGt                                                      >  6:209587/1‑149 (MQ=255)                                                                                               tCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑CCAGGGCTATCGGTCCATGAc                                               >  3:67277/1‑149 (MQ=255)                                                                                                 aGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCTGGCTGGAGGCttcaa                                                                    <  3:128845/127‑4 (MQ=255)                                                                                                 aGTACCTACAAGTCACGGCTACTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCTGGCTGGAGGCttcaa                                                                    >  4:128845/1‑124 (MQ=255)                                                                                                  gTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑CCAGGGCTATCGGTCCATGACTAc                                            <  2:296305/149‑1 (MQ=255)                                                                                                   tACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑CCAGGGCTATCGGTCCATGACTACg                                           >  5:140563/1‑149 (MQ=255)                                                                                                      cTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGGCCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑CCAGGGCTATCGGTCCATGACTACGGGa                                        >  8:56372/1‑149 (MQ=255)                                                                                                      cTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑CCAGGGCTATCGGTCCATGACTACGGGa                                        >  7:304674/1‑149 (MQ=255)                                                                                                       tACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑CCAGGGCTATCGGTCCATGACTACGGGAc                                       >  1:256565/1‑149 (MQ=255)                                                                                                         cAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGGCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTcca                                                                                   <  2:129148/104‑1 (MQ=255)                                                                                                          aaGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑CCAGGGCTATCGGTCCATGACTACGGGACCCg                                    >  1:246717/1‑149 (MQ=255)                                                                                                           aGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑CCAGGGCTATCGGTCCATGACTACGGGACCCg                                    >  2:274852/1‑148 (MQ=255)                                                                                                             tCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑CCAGGGCTATCGGTCCATGACTACGGGACCCGCAc                                 <  3:367557/149‑1 (MQ=255)                                                                                                               aCGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑CCAGGGCTATCGGTCCATGACTACGGGACCCGCACTa                               >  2:237962/1‑149 (MQ=255)                                                                                                                  gctgctgctCAAGCTCAATGCTCGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑CCAGGGCGATCGGGCCATGACTACGGGACCCGCACTGCAg                            <  6:144595/149‑1 (MQ=255)                                                                                                                      ctgctCAAGCTCAATGCACGGTCATTGGTTGGCCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑CCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTc                         <  7:56372/148‑1 (MQ=255)                                                                                                                      ctgctCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑CCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTc                         <  8:304674/148‑1 (MQ=255)                                                                                                                        gctcaagctcaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCttt                                                                      >  7:198954/1‑102 (MQ=255)                                                                                                                        gctcaagctcaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCttt                                                                      <  8:198954/102‑1 (MQ=255)                                                                                                                         ctcaagctcaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑CCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTaa                     >  2:343395/1‑149 (MQ=255)                                                                                                                            aagctcaaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTGCCCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTaa                     <  1:264939/147‑1 (MQ=255)                                                                                                                                  aaTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑CCAGGGCTATCGGTGCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCg            >  5:325740/1‑149 (MQ=255)                                                                                                                                   aTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑CCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCgg           <  1:343395/149‑1 (MQ=255)                                                                                                                                   aTGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑CCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCgg           <  1:237962/149‑1 (MQ=255)                                                                                                                                      cACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAAATGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑CCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGca         >  2:99128/1‑147 (MQ=255)                                                                                                                                       aCGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑CCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCa        >  3:281286/1‑148 (MQ=255)                                                                                                                                             aTTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑CCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCt  <  6:27761/148‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                             CTTTCTGCGCGAGAAATGACGGTGCTGCAACAACTGGCGCGAGGGCGTTCGAACAAGGAGATTGCCGACAGCATGTACCTCAGTAACAAGACGGTCAGTACCTACAAGTCACGGCTGCTGCTCAAGCTCAATGCACGGTCATTGGTTGACCTGATCGAGCTGGCGCAACGGCAAGGTCTGAACTGATCGCCAAGCCTTGAAGCCTCCAGCCAGGAGGCTTTG‑CCAGGGCTATCGGTCCATGACTACGGGACCCGCACTACAGGTCGTAATCGAAATCGGCCAACTGCT  >  NC_002947/1249542‑1249829

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