Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,196,292	G→A	9.7%	intergenic (+134/+235)	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,292	0	G	A	9.7%	165.3 / 7.4	62	intergenic (+134/+235)	PP_5618/PP_5619	Cro/CI family transcriptional regulator/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (27/29);  new base A (2/4);  total (29/33)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.58e-01

CAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACGAGACAGCGAAGAAGCTA  >  NC_002947/4196146‑4196438                                                                                                                                                   |                                                                                                                                                    cAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGgt                                                                                                                                                   <  6:274173/149‑1 (MQ=255)              cGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGACGTCGACACCGTCGACGGCTCgg                                                                                                                                      <  4:150273/149‑1 (MQ=255)              cGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGACGTCGACACCGTCGACGGCTCgg                                                                                                                                      <  1:333255/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                                                                                                                             <  3:312894/142‑1 (MQ=255)                              gCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCggg                                                                                                                             >  4:312894/1‑142 (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                                                  >  1:315892/1‑149 (MQ=255)                                                                                                 tAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCAt                                                  >  7:80791/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                                    >  4:201681/4‑123 (MQ=255)                                                                                                                                         agacGTCGACACCGTCGACGGCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGttt                                    <  3:201681/120‑1 (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     >  6:142713/1‑149 (MQ=255)                                                                                                                                              tCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACGAGACAGCGaagaag     >  7:282671/1‑149 (MQ=255)                                                                                                                                                 aCGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACGAGACAGCGAAGAAGCTa  <  5:51314/149‑1 (MQ=255)                                                                                                                                                 aCGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACGAGACAGCGAAGAAGCTa  <  4:218483/149‑1 (MQ=255)                                                                                                                                                 aCGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACGAGACAGCGAAGAAGCTa  <  1:9730/149‑1 (MQ=255)                                                                                                                                                 aCGGTGGCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACGAGACAGCGAAGAAGCTa  <  7:138726/149‑1 (MQ=255)                                                                                                                                                   |                                                                                                                                                    CAGGGCCTGGTGACGGTAAGCCACTGCTCGCTTCTGCTCTACGCGAGGAGGGTAATGGCCCGACAAGCTCGCATGAGATCGCCTAGCAGATCTTGCTAAGGACTGCGTGCCCTATGGGCCCGGCTGTCGTGAACCGAGCCGTCGACGGTGTCGACGTCTCGGCCCTCCGGGCAGCCATCAGTCACGC‑TTGCGCGCCCCACTTGCCTGCAGCAGAAATTCTGAACATGCTCTGTTCGACCGCTCATGGATTGCTGCGTTTGGGTATGACTCAGAACGAGACAGCGAAGAAGCTA  >  NC_002947/4196146‑4196438

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