Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	3,983,797	A→G	6.0%	Q256R (CAA→CGA)	ilvE →	branched‑chain‑amino acid aminotransferase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	3,983,797	0	A	G	6.0%	200.1 / 2.3	66	Q256R (CAA→CGA)	ilvE	branched‑chain‑amino acid aminotransferase
Reads supporting (aligned to +/- strand):  ref base A (30/32);  new base G (2/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.44e-01

CATCTACCTCGACCCGCTGACCCACACCAAGATCGAGGAAGTCGGTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGT  >  NC_002947/3983648‑3983937                                                                                                                                                      |                                                                                                                                             cATCTACCTCGACCCGCTGACCCACACCAAGATCGAGGAAGTCGGTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCa                                                                                                                                              >  1:812866/1‑150 (MQ=255)      aCCTCGACCCGCTGACCCACACCAAGATCGAGGAAGTCGGTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATcgc                                                                                                                                         >  2:1071675/1‑150 (MQ=255)        cTCGACCCGCTGACCCACACCAAGATCGAGGAAGTCGGCTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCgcgc                                                                                                                                       >  1:731403/1‑150 (MQ=255)         tCGACCCGCTGACCCACACCAAGATCGAGGAAGTCGGTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCAGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGcc                                                                                                                                      >  1:726094/1‑150 (MQ=255)           gACCCGCTGACCCACACCAAGCTCGAGGAAGTCGGTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTg                                                                                                                                    <  1:573029/150‑1 (MQ=255)             cccGCTGACCCACACCAAGATCGAGGAAGTCGGTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTggg                                                                                                                                  <  1:488122/150‑1 (MQ=255)              ccGCTGACCCACACCAAGATCGAGGAAGTCGGTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCAGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGc                                                                                                                                 <  2:726094/150‑1 (MQ=255)                   gACCCACACCAAGATCGAGGAAGCTGTCTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGAc                                                                                                                            <  1:36/150‑1 (MQ=255)                      ccACACCAAGATCGAGGAAGTCGGTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGt                                                                                                                         >  2:780614/1‑150 (MQ=255)                        acacCAAGATCGAGGAAGTCGGTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGa                                                                                                                       <  1:780614/150‑1 (MQ=255)                            cAAGATCGAGGAAGTCGGTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGa                                                                                                                   >  1:498325/1‑150 (MQ=255)                            cAAGATCGAGGAAGTCGGTTCGGCCAACTTATTCGGCATCACCGCAATCGCAGGATTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAAGCGCGCCGGGGCATGACCGTTATCGa                                                                                                                   >  2:470681/1‑150 (MQ=255)                              aGATCGAGGAAGTCGGCTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAg                                                                                                                 <  2:731403/150‑1 (MQ=255)                               gATCGAGGAAGTCGGTTCTCCCAACTTCTTCGTCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAgg                                                                                                                <  2:568163/150‑1 (MQ=255)                                   gAGGAAGTCGGTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGAt                                                                                                            >  1:436175/1‑150 (MQ=255)                                     ggAAGTCGGTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGGCCGTGATCGAAGGCGATGt                                                                                                          >  2:427824/1‑150 (MQ=255)                                             gTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCa                                                                                                  >  2:971937/1‑150 (MQ=255)                                             gTTCGACCAAACTCCTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATCGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCa                                                                                                  <  2:335820/150‑1 (MQ=255)                                                 ggCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATATCGAGATCAGCaa                                                                                              <  2:812866/150‑1 (MQ=255)                                                     aaCTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTg                                                                                          >  1:224341/1‑150 (MQ=255)                                                      aCTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCGATTGCGccagtt                                                                                                                                  <  2:190195/109‑5 (MQ=255)                                                      aCTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCGATTGCGccagtt                                                                                                                                  >  1:190195/1‑105 (MQ=255)                                                       cttcttCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGa                                                                                        <  1:776219/150‑1 (MQ=255)                                                        ttcttcGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGAc                                                                                       <  2:909730/150‑1 (MQ=255)                                                                    aCCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCGATTGCGcca                                                                                                                                     <  2:815611/92‑2 (MQ=255)                                                                    aCCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCGATTGCGcca                                                                                                                                     >  1:815611/1‑91 (MQ=255)                                                                     ccGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAg                                                                          >  2:487475/1‑150 (MQ=255)                                                                      cGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGc                                                                         >  2:898276/1‑150 (MQ=255)                                                                       gCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGcc                                                                        >  1:1003941/1‑150 (MQ=255)                                                                          aacaacGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGt                                                                     >  1:1146305/1‑150 (MQ=255)                                                                            caacGAGTTCGTAACCCCGCAACCGGCCTAGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGc                                                                   <  2:498325/150‑1 (MQ=255)                                                                             aacGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGa                                                                            >  2:810975/1‑140 (MQ=255)                                                                             aacGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGa                                                                            <  1:810975/140‑1 (MQ=255)                                                                               cGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCgg                                                                      >  2:539503/1‑144 (MQ=255)                                                                               cGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCgg                                                                      <  1:539503/144‑1 (MQ=255)                                                                                  gTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCgg                                                             <  2:224341/150‑1 (MQ=255)                                                                                           cccGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGc                                                                                                                                 >  1:817244/1‑73 (MQ=255)                                                                                           cccGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGc                                                                                                                                 <  2:817244/73‑1 (MQ=255)                                                                                               aaaTCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGtt                                                <  1:898276/150‑1 (MQ=255)                                                                                                    ggCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCAc                                           <  1:1071675/150‑1 (MQ=255)                                                                                                    ggCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCCGCCGCTATCAc                                           >  2:305678/1‑150 (MQ=255)                                                                                                           gTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCGTGCGGCACCGCGGCCGTTATCACCCCGATc                                    >  1:319814/1‑150 (MQ=255)                                                                                                               tGCCAGCCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATcggcg                                <  1:470681/150‑1 (MQ=255)                                                                                                                gCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGAt                                                                                                            >  2:462091/1‑73 (MQ=255)                                                                                                                gCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGAt                                                                                                            <  1:462091/73‑1 (MQ=255)                                                                                                                 ccAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCggcggc                              >  1:746570/1‑150 (MQ=255)                                                                                                                  cAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCa                                            >  2:202099/1‑135 (MQ=255)                                                                                                                  cAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCa                                            <  1:202099/135‑1 (MQ=255)                                                                                                                      cATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGa                         <  2:1003941/150‑1 (MQ=255)                                                                                                                        tCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGGCCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGt                       <  1:427824/150‑1 (MQ=255)                                                                                                                           cccGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTtcgt                                                                               <  1:517081/91‑1 (MQ=255)                                                                                                                           cccGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTtcgt                                                                               >  2:517081/1‑91 (MQ=255)                                                                                                                           cccGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACa                    >  2:855235/1‑150 (MQ=255)                                                                                                                            ccGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCgaga                                                                                                     <  1:825941/68‑1 (MQ=255)                                                                                                                            ccGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCgaga                                                                                                     >  2:825941/1‑68 (MQ=255)                                                                                                                              gCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACg                 <  1:971937/150‑1 (MQ=255)                                                                                                                                ctctCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGc                                                                                                >  1:781988/1‑69 (MQ=255)                                                                                                                                ctctCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGc                                                                                                <  2:781988/69‑1 (MQ=255)                                                                                                                                 tctcGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGaa                                                                           <  2:1035354/89‑1 (MQ=255)                                                                                                                                 tctcGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGaa                                                                           >  1:1035354/1‑89 (MQ=255)                                                                                                                                 tctcGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCa              >  2:33405/1‑150 (MQ=255)                                                                                                                                 tctcGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCa              >  1:293027/1‑150 (MQ=255)                                                                                                                                  ctcGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCaa             <  2:1146305/150‑1 (MQ=255)                                                                                                                                   tcGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAg            <  1:33405/150‑1 (MQ=255)                                                                                                                                             aaCTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACgtgt  <  1:855235/150‑1 (MQ=255)                                                                                                                                                  gcgcAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTg                                                                                          <  1:253707/57‑1 (MQ=255)                                                                                                                                                  gcgcAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTg                                                                                          >  2:253707/1‑57 (MQ=255)                                                                                                                                                      |                                                                                                                                             CATCTACCTCGACCCGCTGACCCACACCAAGATCGAGGAAGTCGGTTCGGCCAACTTCTTCGGCATCACCGCCAACAACGAGTTCGTCACCCCGAAATCGGCCTCGGTACTGCCAGGCATCACCCGCCTCTCGCTGATGGAACTGGCGCAATCGCGCCTGGGCATGACCGTGATCGAAGGCGATGTCGAGATCAGCAAGCTGGACCGCTTCGTCGAAGCCGGTGCCTGCGGCACCGCGGCCGTTATCACCCCGATCGGCGGCATCGAGTACAACGGCAAGCTGCACGTGT  >  NC_002947/3983648‑3983937

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑