Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,621,118	0	A	C	44.1%	96.7 / 14.4	60	V7G (GTG→GGG)	eamA	cysteine and O‑acetyl‑L‑serine efflux system
Reads supporting (aligned to +/- strand):  ref base A (25/8);  new base C (0/26);  total (25/35)
Fisher's exact test for biased strand distribution p-value = 5.04e-10
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTT  >  NC_000913/1620970‑1621267                                                                                                                                                     |                                                                                                                                                      tACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACAcc                                                                                                                                                     >  1:75891/1‑151 (MQ=255)  aCTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACAcc                                                                                                                                                     <  1:235139/150‑1 (MQ=255)         gTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCttt                                                                                                                                             <  1:532750/151‑1 (MQ=255)         gTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCttt                                                                                                                                             >  2:2842/1‑151 (MQ=255)          tCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCtttt                                                                                                                                            >  1:365938/1‑151 (MQ=255)          tCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCtttt                                                                                                                                            >  1:256711/1‑151 (MQ=255)              gCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTcgcg                                                                                                                                        >  2:219525/1‑151 (MQ=255)                  cAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACAt                                                                                                                                    <  1:758693/151‑1 (MQ=255)                  cAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACAt                                                                                                                                    >  1:520987/1‑151 (MQ=255)                                  aaaaGCGACCAGCAAAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCAGACGACCCCGCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAAt                                                                                                                    <  2:199181/151‑1 (MQ=255)                                    aggCGACCAGCATAAAGCGCAACCCGGCCACCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTACCCCCCATCCGCCCCCTCCCAGTAGCGCCACCCCCCCATCTTTCCCCGACATTCTACCGCCTCTGAAttt                                                                                                                  <  2:440218/149‑1 (MQ=255)                                     aGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCATCTGAAttt                                                                                                                  >  1:812529/1‑150 (MQ=255)                                         aCCAGAATAAAGCGCAAACCGGCCAGCATCAGTGGTGGCATGTTATGAAGCCCCACTTTGATGACCAAAAAATTTACCCCCCATACGCCCCCTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCAt                                                                                                               <  2:164988/149‑1 (MQ=255)                                           cAGCAAAAAGCGCAACCCGGCCAGCATCAGCGGTGGCATTTTATGAAGCCCCACTTTGATGACCACAAAATTTACCCCCCATCCGCCCCCTCCCATTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCt                                                                                                             <  2:7285/149‑1 (MQ=255)                                            aGCATAAAGCGCAACCCGGCCAGCATCAGCGGTGGCATTTTATGAAGCCCCACTTGGATGACCACAAAATTTACCCCCCATACGACCCCTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCtttt                                                                                                          <  1:573895/151‑1 (MQ=255)                                            aGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCtttt                                                                                                          >  1:226674/1‑151 (MQ=255)                                               aTAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTGGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGt                                                                                                        >  2:692902/1‑150 (MQ=255)                                               aTAAAGAGAAAACCGCCCACCACCCGCGGTGGTATTTTATGAACCCCAATTTTAATGACCAAAAAATTTCGCCCCAATACGCCCCCTCCCAATACCCCCACCCCCCCACCTTTCCCCTCCTTTCTACCGCCTCTGAATTTCATCtttt                                                                                                          <  2:836860/148‑1 (MQ=255)                                                  aaGCGCAACCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGCCCCCTACCAGTAGCGCCACCCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGc                                                                                                    <  1:444408/151‑1 (MQ=255)                                                  aaGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGc                                                                                                    >  1:147803/1‑151 (MQ=255)                                                   aGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAg                                                                                                     >  2:424651/1‑149 (MQ=255)                                                    gcccAACCCGCCCAGCATAAGGGGTGGCATGTTATAAACCCCCCTTTTGTTGACCACAAAATTTACCCCCCAAACGACCTCTCCCATTAGCTCCAACACCCCATCTTTTCCCAACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGc                                                                                                    <  2:878581/146‑1 (MQ=255)                                                     cccAACCCGCCCACCACCAGCGGGGGCATGTAATGAGCCCCCATTTTGATGACCAAAAAATTTAGCCCCCAAACGCCCTCTCCCGGTAGCTCCAACCCCCCATCTTTTCCCGACTTGCTACCGCCTCTGAATTTCATCTTTTGTAAGc                                                                                                    <  1:403784/146‑1 (MQ=255)                                                        aacccGGCCAGCATCAGCGGTGGCATTTAATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATCCGCCCCCTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCaa                                                                                              <  2:620705/148‑1 (MQ=255)                                                        aaaCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCCCTTTGATGACCACAAATTTTAGCCCCCATACGCCCACTCCCAGTAGCGCCAACCCCCAATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAAt                                                                                                 <  1:40592/148‑1 (MQ=255)                                                         aaCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAAc                                                                                             >  2:271640/1‑151 (MQ=255)                                                         aaCCGGCCAGAATCAGCGGTGGCATGTTATGAAGCCCCATTTTGATGACCACAAAATTTAGCCCCCATCCGCCCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAAc                                                                                             <  1:18668/151‑1 (MQ=255)                                                          aCCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACt                                                                                            >  2:864646/1‑151 (MQ=255)                                                               ccAGCATCGGCGGGGGCATTTTATGAAGCCCCACTTTGATGCCCCAAAAATTTAGCCCCCATACGCCCCCTACCATTAGCGCCACCCCCCCATTTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCt                                                                                       <  1:208724/151‑1 (MQ=255)                                                                cAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTg                                                                                      >  2:580856/1‑151 (MQ=255)                                                                 agcacaaggggTGGCATGTTATGAAGCCCCCCTTTGATGCCCACAAAATTTAGCCCCCATACGCCCCCTACCATTAGCGCCACCCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCt                                                                                       <  2:66193/140‑1 (MQ=255)                                                                 aGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTg                                                                                      >  1:786036/1‑150 (MQ=255)                                                                    aTCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTCCCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAt                                                                                   <  1:30752/150‑1 (MQ=255)                                                                    aTCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAt                                                                                   <  1:436304/150‑1 (MQ=255)                                                                    aTCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCAATACGCCCACCACCAGTAGCGCCAACACCCCCTCTTTTCGCGACATTTTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAtt                                                                                  <  1:649637/151‑1 (MQ=255)                                                                      cAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCAACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTa                                                                                <  1:622974/151‑1 (MQ=255)                                                                      cAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAAttt                                                                                 <  1:747046/150‑1 (MQ=255)                                                                        gCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACt                                                                              >  1:49703/1‑151 (MQ=255)                                                                         cGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACtt                                                                             >  2:433118/1‑151 (MQ=255)                                                                           gTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACtttt                                                                           <  1:732996/151‑1 (MQ=255)                                                                                attttATGAGCCCCCACTTTGATGACCAAAAAATTTACCCCCCATACCCCCCCTACCATTAGCGCCAACCCCCCATCTTTCCACGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCttt                                                                       <  2:775816/147‑1 (MQ=255)                                                                                        agcccccGCTTTAAGGACCACAAATTTTAGCCCCCATACGCCCCCTTCCATTCGCGCCAACCCCCCATCTTTTCGCGCCATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGtt                                                                <  2:495857/146‑1 (MQ=255)                                                                                        aaGCCCCTTTTTGATTCCCACAAAATTTAGCCCCCAAACGCCCCCAACCAGTAGCCCCAACCCCCCATCTTTTCGCAACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGttt                                                               <  2:636480/150‑2 (MQ=255)                                                                                                     atgaCCACAAAATTTAGCCCCCATCCGCCCCCTCCCATTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGcc                                                 <  2:308122/151‑1 (MQ=255)                                                                                                      tgaCCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCg                                                >  1:626636/1‑151 (MQ=255)                                                                                                              aaaTTTAGCCCCCATACGCCCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacgac                                        <  1:859041/151‑1 (MQ=255)                                                                                                               atttgaccccccATACTACCCCCCCCAGTAGCCCCAACCCCCCCTCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTacgacg                                       <  2:685267/144‑1 (MQ=255)                                                                                                                     gCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCAtt                                  >  1:340972/1‑150 (MQ=255)                                                                                                                       ccccATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGCAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAAt                               >  1:546080/1‑151 (MQ=255)                                                                                                                         ccATACGCCCACTACCAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgc                             <  1:205926/151‑1 (MQ=255)                                                                                                                          cATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgcg                            >  1:176569/1‑151 (MQ=255)                                                                                                                           aTACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATgcgc                           >  1:610338/1‑151 (MQ=255)                                                                                                                                 cctctcccAGTACCGCAAACGCCCCTTCTTTTCCCGGTTTTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAAt                       <  2:708661/143‑1 (MQ=255)                                                                                                                                   aCTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGt                   >  2:293489/1‑151 (MQ=255)                                                                                                                                       ccAGTAGCGCCAACCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCt               <  1:788802/151‑1 (MQ=255)                                                                                                                                       ccAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCt               >  2:471279/1‑151 (MQ=255)                                                                                                                                            aGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATAc           <  2:69896/150‑1 (MQ=255)                                                                                                                                            aGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACt          >  1:113045/1‑151 (MQ=255)                                                                                                                                              cgcCACCCCCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACtt         <  1:662077/150‑1 (MQ=255)                                                                                                                                              ccccAACCCCCCATCTTTCCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACt          <  2:478740/147‑1 (MQ=255)                                                                                                                                                    cacCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGAttt  <  2:462388/151‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                                      TACTTTCGGTCGTGCGACAAAAAAGATAGCCGGAAAAGCGACCAGCATAAAGCGCAAACCGGCCAGCATCAGCGGTGGCATGTTATGAAGCCCCACTTTGATGACCACAAAATTTAGCCCCCATACGACCACTACCAGTAGCGCCAACACCCCATCTTTTCGCGACATTCTACCGCCTCTGAATTTCATCTTTTGTAAGCAATCAACTTAGCTGAATTTACTTTTCTTTAACAGTTGATTCGTTAGTCGCCGGTTACGACGGCATTAATGCGCAAATAAGTCGCTATACTTCGGATTT  >  NC_000913/1620970‑1621267

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