Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	1,156,608	C→A	6.7%	A32A (GCC→GCA)	yihS →	aldose‑ketose isomerase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	1,156,608	0	C	A	6.7%	213.1 / 4.3	68	A32A (GCC→GCA)	yihS	aldose‑ketose isomerase
Reads supporting (aligned to +/- strand):  ref base C (29/34);  new base A (2/3);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.77e-01

CCCGCTGGCTAGAATTCGCAAACGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCAGGGCATCCCCGGTAGCCTGGC  >  NC_002947/1156466‑1156754                                                                                                                                                |                                                                                                                                                     cccGCTGGCTAGAATTCGCAAACGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCAGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGcc                                                                                                                                                >  2:22124/1‑149 (MQ=255)    gCTGGCTAGAATTCGCAAACGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTGGCAAAGGccagc                                                                                                                                             <  3:143226/149‑4 (MQ=255)          tAGAATTCGCAAACGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTGGCAAAGGccagcag                                                                                                                                           >  1:244008/1‑140 (MQ=255)          tAGAATTCGCAAACGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTGGCAAAGGccagcag                                                                                                                                           <  2:244008/145‑6 (MQ=255)          tAGAATTCGCAAACGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTGGCAAAGGccagca                                                                                                                                            >  1:184351/1‑140 (MQ=255)          tAGAATTCGCAAACGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTGGCAAAGGccagca                                                                                                                                            <  2:184351/144‑5 (MQ=255)            gAATTCGCAAACGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAc‑c                                                                                                                                    <  7:194726/149‑1 (MQ=255)              aTTCGCAAACGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGCAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAc‑cc                                                                                                                                   >  4:40114/1‑148 (MQ=255)                  gCAAACGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACgg                                                                                                                              >  3:83111/1‑149 (MQ=255)                     aaCGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGttt                                                                                                                           <  1:343665/149‑1 (MQ=255)                       cGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCg                                                                                                                         >  6:80116/1‑149 (MQ=255)                        gACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTTGCCTAAGCACCAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCg                                                                                                                         <  1:86479/149‑1 (MQ=255)                               aaGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTAGGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTg                                                                                                                 >  4:195591/1‑149 (MQ=255)                                  gATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGAt                                                                                                              <  3:170367/149‑1 (MQ=255)                                          tGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATgg                                                                                                      >  3:156465/1‑149 (MQ=255)                                              ccATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCaa                                                                                                  >  7:29796/1‑149 (MQ=255)                                              ccATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCaa                                                                                                  >  3:307099/1‑149 (MQ=255)                                                  gACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCt                                                                                                                  >  8:298105/1‑129 (MQ=255)                                                  gACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCt                                                                                                                  <  7:298105/129‑1 (MQ=255)                                                    ccccccTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CACCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGcc                                                                                            >  4:118028/1‑149 (MQ=255)                                                     cccccTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGccc                                                                                           <  1:269188/149‑1 (MQ=255)                                                            aaCCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTACTCCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTg                                                                                     <  1:86826/149‑1 (MQ=255)                                                               cTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCg                                                                                 >  1:363745/1‑149 (MQ=255)                                                                tGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGt                                                                                >  3:141004/1‑149 (MQ=255)                                                                 gTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTg                                                                               >  4:232752/1‑149 (MQ=255)                                                                 gTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAAAGGGCAACTCGCCCCCGGTGCCCGTg                                                                               >  4:237599/1‑149 (MQ=255)                                                                  tATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGt                                                                                <  5:80116/147‑1 (MQ=255)                                                                  tATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑c                                                                             <  5:36321/149‑1 (MQ=255)                                                                    tgtgAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGGGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑cc                                                                            <  5:321175/148‑1 (MQ=255)                                                                      tgAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGa                                                                          >  6:120735/1‑148 (MQ=255)                                                                       gAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAAAGGGGAACTCCGCCCCGGACCCCCTC‑CCGaat                                                                        >  2:147358/1‑148 (MQ=255)                                                                        aaCAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCg                                                                                        >  4:101118/1‑133 (MQ=255)                                                                        aaCAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCg                                                                                        <  3:101118/133‑1 (MQ=255)                                                                          cAGTTGGCTCAATGCCCCGGCCCACACCTCCTGGCGCGTGGCTGAAGCA‑CAGCGCCTGCT‑GCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCAt                                                                    <  6:211033/149‑1 (MQ=255)                                                                            gTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATc                                                                   <  3:195591/149‑1 (MQ=255)                                                                            gTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATc                                                                   <  4:156465/149‑1 (MQ=255)                                                                                  tCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGATACCATCAACAcc                                                             >  7:49553/1‑149 (MQ=255)                                                                                  tCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCCTCaccacc                                                             >  7:43890/1‑149 (MQ=255)                                                                                    aaagCCCCGGCCCACACCGCCAGGCGCCTGCCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCCCGAAACTAC‑CCGACGGTTACGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CAGAAACCATCAACAccgc                                                           >  8:193812/4‑149 (MQ=255)                                                                                        ccccGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACAccgcccgc                                                       <  5:238052/149‑1 (MQ=255)                                                                                              cccACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTGCCCGAAACCATCAACACCGCCCGCATGAc                                                  <  1:22124/149‑1 (MQ=255)                                                                                                  cacCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCAGGTGCCCGTG‑CCGCAACCATCAACACCGCCCGCATGACCCAtt                                             >  1:113328/1‑149 (MQ=255)                                                                                                   acCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCG‑CCCCGGTGCCCGTG‑CCGAAACCATCAACAccgcccgc                                                       <  2:311888/137‑1 (MQ=255)                                                                                                     cgcctggcgcctggCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCt                                          >  3:360430/1‑149 (MQ=255)                                                                                                     cgcctggcgcctggCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCt                                          >  5:15963/1‑149 (MQ=255)                                                                                                       cctggcgcctggcTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCttt                                        <  5:82703/149‑1 (MQ=255)                                                                                                       cctggcgcctggcTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCttt                                        <  1:147358/149‑1 (MQ=255)                                                                                                          ggcgcctggcCGAAGCA‑CCGCGCCGGCTGGCCTGTGCCCAGGCCGCGAAACTAC‑CCGACGGGTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGcc                                     <  8:43890/149‑1 (MQ=255)                                                                                                          gccgcgtggcTGAACCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTCC‑CCGACGGCTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGcc                                     <  8:49553/147‑1 (MQ=255)                                                                                                           gcgcctggcTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGccc                                    >  3:41643/1‑149 (MQ=255)                                                                                                           gcgcctggcTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGc                                      <  4:760/147‑1 (MQ=255)                                                                                                             gcctggcTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTg                                       <  1:366006/144‑1 (MQ=255)                                                                                                             gcctggcTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTg                                       >  2:366006/1‑144 (MQ=255)                                                                                                             gcctggcTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTg                                  <  6:15963/149‑1 (MQ=255)                                                                                                              cctggcTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTgg                                 <  3:118028/149‑1 (MQ=255)                                                                                                                    tGAAGCA‑CAGCGCCTGCTGGCCTTTGCCTAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATAGCTTTGCCCTGGCCCAcc                           >  5:356990/1‑149 (MQ=255)                                                                                                                        gCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTg                         >  3:113002/1‑147 (MQ=255)                                                                                                                        gCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCa                       <  5:120735/149‑1 (MQ=255)                                                                                                                         ca‑caGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCAg                      <  4:141004/149‑1 (MQ=255)                                                                                                                          a‑caGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCAgg                     <  4:360430/149‑1 (MQ=255)                                                                                                                          a‑caGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCAgg                     <  3:237599/149‑1 (MQ=255)                                                                                                                          a‑caGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCAgg                     <  3:232752/149‑1 (MQ=255)                                                                                                                                gcCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCAGGGCATc                <  2:363745/149‑1 (MQ=255)                                                                                                                                   tgctgGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCAg                      <  6:327560/140‑1 (MQ=255)                                                                                                                                   tgctgGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCAg                      >  5:327560/1‑140 (MQ=255)                                                                                                                                        gCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCAGGGCATCCCCGGTa         <  4:113002/148‑1 (MQ=255)                                                                                                                                        gCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCAGGGCATCCCCGGTAg        <  4:41643/149‑1 (MQ=255)                                                                                                                                         ccTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCAGGGCATCCCCGGTAGc       >  5:341980/1‑149 (MQ=255)                                                                                                                                               ccAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCAGGGCATCCCCGGTAGCCTGGc  <  3:40114/148‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                     CCCGCTGGCTAGAATTCGCAAACGACTGGCAAGGATCAAGCTGGCCCATGACCCCCCTCAACCTGTATGTGAACAGTTGGCTCAATGCCCCGGCCCACACCGCCTGGCGCCTGGCTGAAGCA‑CAGCGCCTGCTGGCCTTTGCCAAGGCCGCGAAACTAC‑CCGACGGTTTCGGTAACCTGGATGCCAATGGGCAACTCGCCCCCGGTGCCCGTG‑CCGAAACCATCAACACCGCCCGCATGACCCATTGCTTTGCCCTGGCCCACCTGCAGGGCATCCCCGGTAGCCTGGC  >  NC_002947/1156466‑1156754

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