Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,196,302	T→G	9.6%	intergenic (+144/+225)	PP_5618 → / ← PP_5619	Cro/CI family transcriptional regulator/hypothetical protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,196,302	0	T	G	9.6%	178.4 / 8.6	62	intergenic (+144/+225)	PP_5618/PP_5619	Cro/CI family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (28/28);  new base G (2/4);  total (30/32)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

CGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACGAGACAGCGAAGAAGCTAGGACTGCCCACA  >  NC_002947/4196159‑4196450                                                                                                                                                |                                                                                                                                                      cGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGACGTCGACACCGTCGACGGCTCgg                                                                                                                                                  <  1:333255/149‑1 (MQ=255) cGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGACGTCGACACCGTCGACGGCTCgg                                                                                                                                                  <  4:150273/149‑1 (MQ=255)         cACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCgg                                                                                                                                          <  8:150185/149‑1 (MQ=255)          aCTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCggg                                                                                                                                         >  7:154980/1‑149 (MQ=255)            tGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCggg                                                                                                                                         >  3:127690/1‑147 (MQ=255)              cTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAg                                                                                                                                      <  6:201918/148‑1 (MQ=255)              cTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCCGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGc                                                                                                                                     <  1:91331/149‑1 (MQ=255)                 gCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCggg                                                                                                                                         >  4:312894/1‑142 (MQ=255)                 gCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCggg                                                                                                                                         <  3:312894/142‑1 (MQ=255)                     cTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCa                                                                                                                                <  4:110545/147‑1 (MQ=255)                         tctACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGCGCCGTCGACGGTGTCGCCGTCTCGGCCCTCCGGGCCGCCAGCATTCCCg                                                                                                                          >  8:138726/1‑149 (MQ=255)                         tctACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACg                                                                                                                          >  4:1080/1‑149 (MQ=255)                         tctACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACg                                                                                                                          >  6:51314/1‑149 (MQ=255)                            aCGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑tt                                                                                                                      >  2:239085/1‑149 (MQ=255)                              gcgAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTgc                                                                                                                    <  3:74759/149‑1 (MQ=255)                                 aggaggGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTgcgcg                                                                                                                 >  2:311175/1‑149 (MQ=255)                                          aTGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACtt                                                                                                         <  8:160359/148‑1 (MQ=255)                                              cccGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTg                                                                                                    <  4:127690/149‑1 (MQ=255)                                               ccGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTgc                                                                                                   <  7:140042/149‑1 (MQ=255)                                               ccGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTg                                                                                                    >  6:16915/1‑148 (MQ=255)                                                      gCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGaa                                                                                            <  4:137971/149‑1 (MQ=255)                                                      gCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGa                                                                                             >  6:182483/1‑148 (MQ=255)                                                       cTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGaa                                                                                            >  2:340474/1‑148 (MQ=255)                                                           cATGAGATCGCCTAGCATATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCt                                                                                       <  3:94746/149‑1 (MQ=255)                                                            aTGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATAAGTCACGC‑TTGAGCGACCCACTTGCCTGCAGCAGAAATTCt                                                                                       >  6:224348/1‑148 (MQ=255)                                                             tGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGGGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGa                                                                                     <  8:277279/149‑1 (MQ=255)                                                               agaTCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGaa                                                                                    <  2:310962/148‑1 (MQ=255)                                                                    gCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGct                                                                              <  3:1080/149‑1 (MQ=255)                                                                                  gCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCt                                                                 >  3:182007/1‑148 (MQ=255)                                                                                    tAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCAt                                                              >  7:80791/1‑149 (MQ=255)                                                                                    tAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCAt                                                              >  1:315892/1‑149 (MQ=255)                                                                                    tAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTGCGCCCGCTCAt                                                              >  6:270750/1‑149 (MQ=255)                                                                                      aGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCAt                                                              <  1:239085/147‑1 (MQ=255)                                                                                      aGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATgg                                                            <  8:80791/149‑1 (MQ=255)                                                                                                      tGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGGCCGCTCATGGATTGCTGCGTTTgg                                              <  5:23979/147‑1 (MQ=255)                                                                                                      tGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTgg                                              <  5:16915/147‑1 (MQ=255)                                                                                                      tGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTgg                                              <  5:224348/147‑1 (MQ=255)                                                                                                        ggCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCAGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTa                                           <  8:154980/148‑1 (MQ=255)                                                                                                         gCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTAt                                          >  4:287447/1‑148 (MQ=255)                                                                                                         gCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTAt                                          >  2:9730/1‑148 (MQ=255)                                                                                                              gCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACt                                      >  3:218483/1‑147 (MQ=255)                                                                                                                    gTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGt                                            >  4:266210/1‑135 (MQ=255)                                                                                                                    gTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGt                                            <  3:266210/135‑1 (MQ=255)                                                                                                                    gTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAAc                                >  1:364829/1‑147 (MQ=255)                                                                                                                     tGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCCCCCCACTTGCCTGCAGCAGAAAATCTGAACATGATCAGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAAc                                >  1:75010/1‑146 (MQ=255)                                                                                                                        aCCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGTCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACGAGACa                          <  1:311175/149‑1 (MQ=255)                                                                                                                           gagaCGTCGACACCGTCGACGGCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCa                                    >  4:360476/1‑136 (MQ=255)                                                                                                                           gagaCGTCGACACCGTCGACGGCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCa                                    <  3:360476/136‑1 (MQ=255)                                                                                                                            agacGTCGACACCGTCGACGGCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGttt                                                <  3:201681/120‑1 (MQ=255)                                                                                                                            agacGTCGACACCGTCGACGGCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGttt                                                >  4:201681/4‑123 (MQ=255)                                                                                                                              ccGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACGAGACAGCgaa                     <  2:315892/148‑1 (MQ=255)                                                                                                                                gTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACGAGACAGCgaag                    >  7:186086/1‑147 (MQ=255)                                                                                                                                gTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACGAGACAGCgaag                    >  7:191329/1‑147 (MQ=255)                                                                                                                                 tCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGCGTTGCGGGCCCCACTTGCCTGCAGCAGAAATTATGAACATGGTCTGTTCGAc                                                                     >  1:260981/1‑98 (MQ=38)                                                                                                                                 tCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGAc                                                                     <  2:260981/97‑1 (MQ=255)                                                                                                                                 tCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACGAGACAGCGaagaag                 >  7:282671/1‑149 (MQ=255)                                                                                                                                 tCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACGAGACAGCGaagaag                 >  6:142713/1‑149 (MQ=255)                                                                                                                                    aCGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACGAGACAGCGAAGAAGCTa              <  1:9730/149‑1 (MQ=255)                                                                                                                                    aCGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACGAGACAGCGAAGAAGCTa              <  5:51314/149‑1 (MQ=255)                                                                                                                                    aCGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACGAGACAGCGAAGAAGCTa              <  4:218483/149‑1 (MQ=255)                                                                                                                                    aCGGTGGCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACGAGACAGCGAAGAAGCTa              <  7:138726/149‑1 (MQ=255)                                                                                                                                            gACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGGTGCGTTTGGGTATGACTCAGAACGAGACAGCGAAGAAGCTAGGACTGcc      >  4:240170/1‑149 (MQ=255)                                                                                                                                                tctcGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACGAGACAGCGAAGAAGCTAGGACTGCCcaca  <  1:340474/149‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                      CGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACGAGACAGCGAAGAAGCTAGGACTGCCCACA  >  NC_002947/4196159‑4196450

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