Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	2,776,115	G→T	5.3%	L154L (CTG→CTT)	sotB →	sugar efflux transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	2,776,115	0	G	T	5.3%	239.4 / 3.5	76	L154L (CTG→CTT)	sotB	sugar efflux transporter
Reads supporting (aligned to +/- strand):  ref base G (34/38);  new base T (2/2);  total (36/40)
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.67e-01

AGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGATCAAATCGCT  >  NC_002947/2775969‑2776261                                                                                                                                                       |                                                                                                                                                    aGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑Tgg                                                                                                                                                    <  6:127128/148‑1 (MQ=255)   ccACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCC‑‑‑GCCCGCCAAGGCGC‑Tggggc                                                                                                                                                 >  2:292535/1‑147 (MQ=255)     aCTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCctg                                                                                                                                               <  1:192160/149‑1 (MQ=255)      cTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCctgc                                                                                                                                              <  6:249164/149‑1 (MQ=255)         gCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCCAGCGCC‑TTGGCctgctg                                                                                                                                            <  4:250569/148‑1 (MQ=255)         gCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCctgctg                                                                                                                                            <  1:382778/148‑1 (MQ=255)             tCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCCGGCCAAGGCGC‑TGGGCCTGCAGGcccc                                                                                                                                       >  4:248930/1‑147 (MQ=255)                 ggCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑gg                                                                                                                                   >  8:255741/1‑148 (MQ=255)                  gCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGt                                                                                                                                  <  3:55618/148‑1 (MQ=255)                        gCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTAccc                                                                                                                           >  8:234179/1‑149 (MQ=255)                        gCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTAccc                                                                                                                           >  3:282877/1‑149 (MQ=255)                        gCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTAcc                                                                                                                            >  4:42319/1‑148 (MQ=255)                            aGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑CCAGGCCAAGGCGG‑GGGGCCTGCGGGCCCCC‑CCTACTACCCAGGc                                                                                                                       >  1:144689/1‑149 (MQ=255)                             gCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGcc                                                                                                                      <  8:335056/149‑1 (MQ=255)                              cTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCa                                                                                                                     <  2:222189/149‑1 (MQ=255)                               ttCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCctgc                                                                                                                                              >  1:306864/1‑124 (MQ=255)                               ttCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCctgc                                                                                                                                              <  2:306864/124‑1 (MQ=255)                                 cGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATgg                                                                                                                  >  3:65662/1‑149 (MQ=255)                                 cGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATgg                                                                                                                  >  2:276093/1‑149 (MQ=255)                                     aTGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑ACAGGCCAACTCGC‑TGGGCCAGCAGGCCACC‑CGTACTACCCTAGACAAGGTc                                                                                                                >  8:351978/1‑147 (MQ=255)                                      tgctgcTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTg                                                                                                             <  8:231238/149‑1 (MQ=255)                                       gctgctGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGGCTGCTGGCCACC‑GGGACTAACCTGGCCAGGGTCCTgg                                                                                                            >  6:205194/1‑149 (MQ=255)                                        ctgctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGTCCATGGTCCTgg                                                                                                            <  8:326896/148‑1 (MQ=255)                                        ctgctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTgg                                                                                                            <  2:350869/148‑1 (MQ=255)                                        ctgctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTgg                                                                                                            <  3:302063/148‑1 (MQ=255)                                         tgctgGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTggg                                                                                                           >  6:355109/1‑148 (MQ=255)                                              gTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGtactac                                                                                                                             >  6:322157/1‑125 (MQ=255)                                              gTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGtactac                                                                                                                             <  5:322157/125‑1 (MQ=255)                                                 aGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCt                                                                                                  >  4:135578/1‑149 (MQ=255)                                                  gCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTg                                                                                                 <  7:47857/149‑1 (MQ=255)                                                     gCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCTTCCCGCTGGGc                                                                                              >  2:34617/1‑149 (MQ=255)                                                            cATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCAAAGGCAGGCCAAAGCGCGAGGGCCTGCTGGCCACC‑GGTACTAACCTGGCCATGGTCCTGGGCATCCCGCTGGGCCgt                                                                                           >  1:122126/1‑149 (MQ=255)                                                             aTTGCCCTGGCGCATGCGGTGTTCATGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑CCATGCCAAGGCGC‑AGGGGCTGCTGGCCACC‑CGTACTAACCAGGCCAAGGAACATGGCAAACCGCAGGGACgtgtt                                                                                        >  6:211495/1‑146 (MQ=255)                                                                gCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCTTCCCGCTGGGCCGTGtggtgg                                                                                    <  1:34617/148‑1 (MQ=255)                                                                 cccTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACGACCCAGGCCATGGTCCGGGGCCTCCCGCTGGGCCGGGTGGTGGGc                                                                                  >  8:349760/1‑149 (MQ=255)                                                                 cccTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGCCCTGCTGGCCCCC‑GCGACAACCCTCGCCCTGGTCCTGGGCATACCGCAGGGCCCTCtggtggcc                                                                                  >  3:288052/1‑147 (MQ=255)                                                                   cTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGa                                                                                <  3:135578/149‑1 (MQ=255)                                                                   cTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGa                                                                                <  7:351978/149‑1 (MQ=255)                                                                      gcgcATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGta                                                                                                                                 <  3:70002/97‑1 (MQ=255)                                                                      gcgcATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGta                                                                                                                                 >  4:70002/1‑97 (MQ=255)                                                                      gcgcATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAg                                                                               <  3:42319/147‑1 (MQ=255)                                                                          aTGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTGCTGGGCATACCGCTGGGCCGTGTGGTGGGCGCGGCACTg                                                                         >  3:182473/1‑149 (MQ=255)                                                                              ggTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCt                                                                                                  <  5:235522/120‑1 (MQ=255)                                                                                       gTCCATCTCCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTa                                                               <  1:133916/146‑1 (MQ=255)                                                                                       gTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCGGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCa                                                            >  4:357946/1‑149 (MQ=255)                                                                                          cATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCCCGt                                                         >  5:277594/1‑149 (MQ=255)                                                                                          cATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGt                                                         <  4:65662/149‑1 (MQ=255)                                                                                           aTCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTggg                                                                       <  3:93371/134‑1 (MQ=255)                                                                                           aTCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGtt                                                        >  5:361612/1‑149 (MQ=255)                                                                                            tCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTc                                                       <  4:282877/149‑1 (MQ=255)                                                                                              aCCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCt                                                     <  5:355109/149‑1 (MQ=255)                                                                                                cGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGa                                                   <  7:349760/149‑1 (MQ=255)                                                                                                cGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGa                                                   <  2:174271/149‑1 (MQ=255)                                                                                                  ccTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGc                                                                  <  5:120104/132‑1 (MQ=255)                                                                                                      gCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCAtt                                              <  4:244952/148‑1 (MQ=255)                                                                                                          gCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCgg                                         <  2:122126/149‑1 (MQ=255)                                                                                                                  tgtgGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAgcgc                                 <  3:248930/149‑1 (MQ=255)                                                                                                                      gcgcCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACCGGGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAg                              >  1:86055/1‑149 (MQ=255)                                                                                                                      gcgcCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAg                              <  2:86055/148‑1 (MQ=255)                                                                                                                        gccgccGGGCCA‑‑‑GCAGGCCAACCCGC‑GGGGCCGGCTGGCCCCC‑CGTACGACCCTGGCCCAGGGCCGGGGGCTCCCCCGGGGGCGTGTGGTGGGCGCGGCACGGGGCGGGCGTTTCACGTTGCGGCGCATTGCCGGGGTAGGGGTAGccc                           >  6:6162/1‑148 (MQ=255)                                                                                                                          cgccgGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGccacc                         <  7:255741/149‑1 (MQ=255)                                                                                                                               ggCCA‑‑‑GCAGGCCCAGCGCC‑TTGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGGGGCACTGGGCTGGCGTATCACGTTCCt                                                     >  1:175897/1‑116 (MQ=38)                                                                                                                               ggCCA‑‑‑GCAGGCCCAGCGCC‑TTGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGGGGCACTGGGCTGGCGTATCACGTTCCt                                                     <  2:175897/116‑1 (MQ=38)                                                                                                                                  ca‑‑‑gcagGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGGATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTg                 >  6:259302/1‑149 (MQ=255)                                                                                                                                   a‑‑‑gcagGCCCAGCGCC‑TTGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGc                >  2:26696/1‑149 (MQ=255)                                                                                                                                   a‑‑‑gcagGCCAAGGCGC‑TGGGCCGGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATGACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGc                <  1:276093/149‑1 (MQ=255)                                                                                                                                       gcagGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCCTCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGcc               >  7:143893/1‑149 (MQ=255)                                                                                                                                       gcagGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGcc               <  5:205194/149‑1 (MQ=255)                                                                                                                                           gCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGc                                                                  >  3:151096/1‑94 (MQ=255)                                                                                                                                           gCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGc                                                                  <  4:151096/94‑1 (MQ=255)                                                                                                                                           gCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTCGCGCTAGCCACCATGC‑TTGCCTGAt           >  8:155987/1‑148 (MQ=255)                                                                                                                                           gCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGAt           <  7:234179/149‑1 (MQ=255)                                                                                                                                           gCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGAt           >  7:356142/1‑149 (MQ=255)                                                                                                                                            ccAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGAt           >  4:340539/1‑148 (MQ=255)                                                                                                                                                  cgc‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGATCAAATcg    <  5:299209/149‑1 (MQ=255)                                                                                                                                                    c‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGATCAAATcgct  >  5:152519/1‑149 (MQ=255)                                                                                                                                                      tGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGAt           >  2:195881/1‑139 (MQ=255)                                                                                                                                                      tGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGAt           <  1:195881/139‑1 (MQ=255)                                                                                                                                                       |                                                                                                                                                    AGCCACTTGCTGTCCTGGCTTTCGCAAAGCTTCGCCATGCTGCTGGTCAGCCGCATCGGCATTGCCCTGGCGCATGCGGTGTTCTGGTCCATCACCGCCTCGCTGGCGGTGCGTGTGGCGCCGCCGGGCCA‑‑‑GCAGGCCAAGGCGC‑TGGGCCTGCTGGCCACC‑GGTACTACCCTGGCCATGGTCCTGGGCATCCCGCTGGGCCGTGTGGTGGGCGAGGCACTGGGCTGGCGTATCACGTTCCTGAGCATTGCCGGGGTAGCGCTAGCCACCATGCTTTGCCTGATCAAATCGCT  >  NC_002947/2775969‑2776261

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