Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_002947	4,091,449	A→C	6.2%	N20H (AAC→CAC)	garD →	galactarate dehydratase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_002947	4,091,449	0	A	C	6.2%	189.9 / 2.4	61	N20H (AAC→CAC)	garD	galactarate dehydratase
Reads supporting (aligned to +/- strand):  ref base A (30/27);  new base C (2/2);  total (32/29)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.49e-01

TCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAGGGCGAGCCCGTGCGCCGCTACGGGC  >  NC_002947/4091303‑4091593                                                                                                                                                   |                                                                                                                                                 tCCCGTCACTGAAAAAACACACCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCTAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAAc                                                                                                                                                >  3:18331/1‑149 (MQ=255)           gAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTTATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGgtcg                                                                                                                                      >  3:119556/1‑149 (MQ=255)               aaaCACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGgtcgtcgt                                                                                                                                  <  1:122337/149‑1 (MQ=255)                aaCACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACCACGACGTTgtcgtcgtc                                                                                                                                 <  3:77560/149‑1 (MQ=255)                  cacGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCaa                                                                                                                               <  3:233823/149‑1 (MQ=255)                   acGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGgtcgtcgtc                                                                                                                                 >  7:170196/1‑146 (MQ=255)                   acGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGgtcgtcgtc                                                                                                                                 <  8:170196/146‑1 (MQ=255)                      cccGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAacgac                                                                                                                           <  3:59160/149‑1 (MQ=255)                      cccGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGGCGTGGTCGTCATCAacgac                                                                                                                           >  5:145601/1‑149 (MQ=255)                       ccGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCCCGACGACGACAACGTCGTGGTCGTCGTCAacgacg                                                                                                                          >  4:141756/1‑149 (MQ=255)                        cgcgCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGACCCGGACAAGTAGGGCCCAAACATGCAGTGGATCGAACCTGCCGACTCGCCCCGCTCCGGCCGCCTGCCCGACGACGACCACGTCGTGGTCGTCGGCAacgacgc                                                                                                                         >  3:194137/1‑148 (MQ=255)                        cgcgCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTTATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGAcgg                                                                                                                         <  4:119556/149‑1 (MQ=255)                              gggTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCCACGACGGCGGCtt                                                                                                                   >  8:162795/1‑149 (MQ=255)                                 tCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATACCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCCTGGTCGTCGTCAACGACGCCGGCTTggg                                                                                                                >  8:121849/1‑149 (MQ=255)                                    cccGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCg                                                                                                              >  7:191459/1‑148 (MQ=255)                                    cccGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCg                                                                                                              >  5:209585/1‑148 (MQ=255)                                                       cAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGAc                                                                                          >  3:218680/1‑149 (MQ=255)                                                        aaaaGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACCACGACGTTgtcgtc                                                                                                                                    >  6:201358/1‑106 (MQ=255)                                                        aaaaGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACCACGACGTTgtcgtc                                                                                                                                    <  5:201358/106‑1 (MQ=255)                                                        aaaaGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGc                                                                                              <  5:50543/144‑1 (MQ=255)                                                        aaaaGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGc                                                                                              >  6:50543/1‑144 (MQ=255)                                                         aaaGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACgg                                                                                        >  3:228417/1‑149 (MQ=255)                                                          aaGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGc                                                                                       <  6:209585/149‑1 (MQ=255)                                                                        gACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCg                                                                          >  5:277138/1‑148 (MQ=255)                                                                         aCAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCg                                                                          <  4:115181/147‑1 (MQ=255)                                                                          cAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAg                                                                       <  4:228417/149‑1 (MQ=255)                                                                           aaGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAgg                                                                      >  3:95085/1‑149 (MQ=255)                                                                           aaGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAgg                                                                      >  5:170943/1‑149 (MQ=255)                                                                            aGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGc                                                                     >  4:6461/1‑149 (MQ=255)                                                                              tAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAg                                                                       <  7:222493/145‑1 (MQ=255)                                                                              tAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAg                                                                       >  8:222493/1‑145 (MQ=255)                                                                              tAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGt                                                                   >  4:195931/1‑149 (MQ=255)                                                                              tAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGt                                                                   <  3:34018/149‑1 (MQ=255)                                                                                 ggCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGcc                                                                <  6:277138/149‑1 (MQ=255)                                                                                      aaaCATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCaaa                                                           <  1:107097/149‑1 (MQ=255)                                                                                              aGTTGATCGAACATTCTGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAg                                                   >  1:156812/1‑149 (MQ=255)                                                                                                ttGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGt                                                 >  4:58572/1‑149 (MQ=255)                                                                                                    tCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGc                                              >  4:222864/1‑148 (MQ=255)                                                                                                       aaCATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCAcc                                          <  6:170943/149‑1 (MQ=255)                                                                                                        aCATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCg                                         >  8:57887/1‑149 (MQ=255)                                                                                                          aTTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCt                                                                                     <  1:197104/103‑1 (MQ=255)                                                                                                          aTTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCt                                                                                     >  2:197104/1‑103 (MQ=255)                                                                                                          aTTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCgtg                                       <  7:57887/149‑1 (MQ=255)                                                                                                           ttCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCgtgt                                      <  3:58572/149‑1 (MQ=255)                                                                                                               gACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAg                                                   >  5:267729/1‑132 (MQ=255)                                                                                                               gACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAg                                                   <  6:267729/132‑1 (MQ=255)                                                                                                                    gCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGcc                              >  7:77376/1‑148 (MQ=255)                                                                                                                     ccccGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCaa                            <  3:6461/149‑1 (MQ=255)                                                                                                                       ccGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAgg                          <  3:222864/149‑1 (MQ=255)                                                                                                                         gCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAGGGc                        <  3:195931/149‑1 (MQ=255)                                                                                                                         gCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAGGGc                        <  4:218680/149‑1 (MQ=255)                                                                                                                             cGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAGGGCGAg                     <  4:95085/148‑1 (MQ=255)                                                                                                                                   ccTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAGGGCGAGCCCGTgc              >  4:189074/1‑149 (MQ=255)                                                                                                                                   ccTGCACGACGACGACAACGGCGTGGTCGTCGTCAACGACGACGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAGGGCGAGCCCGTgc              <  4:194137/149‑1 (MQ=255)                                                                                                                                      gCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCAt                                  <  5:122163/126‑1 (MQ=255)                                                                                                                                      gCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCAt                                  >  6:122163/1‑126 (MQ=255)                                                                                                                                       cACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTg                                                >  3:246178/1‑111 (MQ=255)                                                                                                                                       cACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTg                                                <  4:246178/111‑1 (MQ=255)                                                                                                                                              acgacAACGTCGTGGTCGTCATCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAGGGCGAGCCCGTGCGCCGCTACggg   <  6:145601/149‑1 (MQ=255)                                                                                                                                                gacAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAGGGCGAGCCCGTGCGCCGCTACGGGc  >  7:101909/1‑148 (MQ=255)                                                                                                                                                gacAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAGGGCGAGCCCGTGCGCCGCTACGGGc  >  1:207513/1‑148 (MQ=255)                                                                                                                                                   |                                                                                                                                                 TCCCGTCACTGAAAAAACACGCCCGCGCCGGGTCGCCCGGCATGGGCGGATGTACAAAAGCCTGAGAACCTGACAAGTAGGGCCCAAACATGCAGTTGATCGAACATTCCGACTCGCCCCGCTACGTCCGCCTGCACGACGACGACAACGTCGTGGTCGTCGTCAACGACGGCGGCTTGGGCGAAGGCGCCCGCTTCGCCGACGGCCTGACCCTGATCGAAGGCGTGCCGCAAAGCCATAAGGTGGCCACCGTGTCCATCGCCAAGGGCGAGCCCGTGCGCCGCTACGGGC  >  NC_002947/4091303‑4091593

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