Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,302,658	0	T	G	43.6%	87.0 / 13.4	56	T141P (ACC→CCC)	mdtO	membrane translocase (MDR) of MdtNOP efflux pump, PET family
Reads supporting (aligned to +/- strand):  ref base T (7/24);  new base G (24/0);  total (32/24)
Fisher's exact test for biased strand distribution p-value = 1.11e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AGCGCCTGATGCATTTGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAACGGTTCGCCGTATGACCATTTGTAGATCAA  >  NC_000913/4302510‑4302802                                                                                                                                                     |                                                                                                                                                 aGCGCCTGATGCATTTGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGttt                                                                                                                                                >  2:1058308/1‑151 (MQ=255) aGCGCCTGATGCATTTGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGttt                                                                                                                                                >  2:196494/1‑151 (MQ=255)  gcgcCTGATGTATTTGCGAAATGGCACGACTGGGAAGCCACAGCGCGCCGATTAACGTCATCAGCAAGGTTAGGTAGAGGCCGACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGGTAGTCGAGCATGGCGGTGAAGGttt                                                                                                                                                <  1:141015/150‑1 (MQ=255)  gcgcCTGATGCATTTGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGgtt                                                                                                                                                >  2:520817/1‑150 (MQ=255)        gATGCATTTGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTTccc                                                                                                                                            >  1:172747/1‑148 (MQ=255)             aTTTGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTaaaaa                                                                                                                                     >  2:454991/1‑148 (MQ=255)              tttGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAGc                                                                                                                                   >  1:778047/1‑151 (MQ=255)                    aaaTGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAAt                                                                                                                             <  2:6524/151‑1 (MQ=255)                    aaaTGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAAt                                                                                                                             <  2:59551/151‑1 (MQ=255)                     aaTGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATg                                                                                                                            <  1:316079/151‑1 (MQ=255)                        ggCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGc                                                                                                                          <  1:322402/150‑1 (MQ=255)                        ggCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGc                                                                                                                          <  1:4297/150‑1 (MQ=255)                        ggCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGc                                                                                                                          <  2:970802/150‑1 (MQ=255)                          cACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGAc                                                                                                                       <  2:1121338/151‑1 (MQ=255)                            cGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGCAAATAGCGACAAAGgcggcg                                                                                                                      >  2:924472/1‑150 (MQ=255)                                tGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTTCCCGCAAAAAGCGACAATGGCGACgg                                                                                                                     >  2:501731/1‑147 (MQ=255)                                  ggAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGTCCGCAAAAAACGAAAATGGAGACGGGGaaa                                                                                                                >  1:877651/1‑150 (MQ=255)                                   gAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTCCCGGAAATAACGAAAAAGGCGACGGGGgaaa                                                                                                               >  1:777802/1‑150 (MQ=255)                                    aaaCCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGTTTGCCCGGAAATAGAGACAAAGGGGACGGCGaaaaa                                                                                                              >  1:1062324/1‑150 (MQ=255)                                         acaGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAAGGTTGCCCGGAAATAACGACAATGGCGACGGCGgaaaa                                                                                                              >  1:914758/1‑145 (MQ=255)                                             cACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGcccc                                                                                                    <  2:347379/151‑1 (MQ=255)                                              aCGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGTCCGTGAATAGGGACAAAGGCGACGGCGAAAAAAACCAGgccc                                                                                                    >  2:984237/1‑150 (MQ=255)                                              aCGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGTAAATAGCGACAATGGCGGCGGCGAAAAAGAACCGcccc                                                                                                    >  2:850499/1‑150 (MQ=255)                                               cGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCaa                                                                                                  <  1:532273/151‑1 (MQ=255)                                                   gATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAACGACAAATGAGAAGGCGGAAAAAACCAGcccccaacg                                                                                               >  2:1064008/1‑145 (MQ=255)                                                    aTTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAACGACAAAGGCGACGGCGGAAAAAAACAGCCCCAAGCg                                                                                               >  2:571403/1‑149 (MQ=255)                                                       aaCGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATgg                                                                                           <  1:599199/150‑1 (MQ=255)                                                         cGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTTCCCGCAAATAACGACAAAGGCGACGGCGAAAAAGACCAGACCCCAACGAAggg                                                                                          >  2:250736/1‑149 (MQ=255)                                                            catcaGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGATTGCCCGGAAAAAACGACAATGGCGCCGGAGGAAAAGACACGCCCCACGCGACggggg                                                                                        >  1:79112/1‑146 (MQ=255)                                                                  cAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAAGGTTGCCCGGAAATAGCGACAATGGCGACGGCGAAAAAAACCAGACCCAAACGATGGGGGCGCATCaa                                                                                >  1:149572/1‑150 (MQ=255)                                                                     ggTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAAGGTTGCCCGGAAATAACGACAATGGCGACCGCGAAAAAGACCAGCCCCAAGCGATGGGGGCGCCTCaaaaa                                                                             >  2:538336/1‑150 (MQ=255)                                                                        tGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGGAAAAAGACACAAAAGCGACGGCGGAAAAAACCAGCCCCCAACGAGGGGTGGGCATCaaaaaaa                                                                           >  1:620168/1‑147 (MQ=255)                                                                         ggATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAGCGACAAAGGCGGCGGCGGAAAAGACCAGCCCCCAGCGATGGGGGCGCATCAAAAACAAGc                                                                        >  2:681509/1‑151 (MQ=255)                                                                          gATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCCTAAATAGCGACAATGGCGCCGGCGAAAAAAAACCGCCCCCAACGATGGGTGCGCATCaaaaaaa                                                                           >  1:1082391/1‑145 (MQ=255)                                                                             agagGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGcc                                                                    <  2:595463/151‑1 (MQ=255)                                                                                   caacaaCGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGCAAAAAGAGAAAATGGCGACGGGGAAAAAGAACAGCCCCAAACGATGGGGGCGCCACAAAAAAAAGCAGccccac                                                                >  2:1132062/1‑146 (MQ=255)                                                                                       aacGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGGAAAAAGCGACAATGGCGACGGCGGAAAAAACCAGCCCCAAACGATGGGTGCGCATCCAAAAAAATCCGCCCATCAgg                                                             >  2:1091943/1‑148 (MQ=255)                                                                                                 cacaGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAAAAGCGACAATGGCGACGGCGAAAAAGAACAGCCCCCAGCGATGGGGGCGCATCAAAAACATGCAGCCCAACAGGATCGGTcc                                                     >  2:302574/1‑146 (MQ=255)                                                                                                  acaGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCgat                                               <  1:876426/151‑1 (MQ=255)                                                                                                   caGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCgatg                                              <  1:248485/151‑1 (MQ=255)                                                                                                   caGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCgat                                               <  1:108769/150‑1 (MQ=255)                                                                                                       gTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCa                                          <  2:1025555/151‑1 (MQ=255)                                                                                                       gTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCa                                          <  2:820820/151‑1 (MQ=255)                                                                                                       gTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAAAAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCCAACCATGGGGGGGCATCAAAAAAAAGCAGCCCCTCAGGATCGGGCCGGCGATGGTc                                           >  2:785527/1‑150 (MQ=255)                                                                                                        ttAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTTCCCCTCAAAAAAGACAAAGGCGACGGCGGAAAAGACCAGCCCCAAACGATGGGGGCGCATCAAAAAAATGCAGCCCCTCAGGAACGGGCCGGCGATGATCa                                          >  1:827307/1‑150 (MQ=255)                                                                                                          aaGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCCTAAATAGCGGCAAAGGCGACCGCGAAAAAGACCAGCCCCCAGCGATGGGGGGGCATCAAAAACATGCAGCCCCTCCGGAACGGGCCGGCGATGATCaa                                         >  1:176791/1‑149 (MQ=255)                                                                                                          aaGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATc                                       <  2:79117/151‑1 (MQ=255)                                                                                                          aaGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATc                                       <  1:1027307/151‑1 (MQ=255)                                                                                                           aGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCCTAAATAGCGACAATGGCGACGGCGAAAAAAACCCGCCCCCAGCGATGGGGGCGCATCAAAAAAATGCAGCCCCTCAGGATCGGTCCGGCGATGATCAATCg                                      >  2:668019/1‑151 (MQ=255)                                                                                                             cgcACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGa                                    <  1:115297/151‑1 (MQ=255)                                                                                                                  gACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTCCCGGAAAAAACGACAAAGGGGACGGCGAAAAAAACCAACCCCCAACGATGGGGGCGCATCAAAAAAATGCAGCCCATCAGGATCGGGCCGGCGGAGATCAAACGGGTCaa                                >  2:1060371/1‑150 (MQ=255)                                                                                                                  gACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAAc                               <  1:135038/151‑1 (MQ=255)                                                                                                                     cACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGTAAAAAGCGACAATGGCGGCTGCGGAAAAAACCAGCCCCCAGCGATGGGTGCGCATCAAAAAAATGCAACCCATCAAGAACGGTCCCGGGATGATTAATCGGATCAAccg                             >  2:814030/1‑148 (MQ=255)                                                                                                                      aCTTCCGGATAGTCGAGCATGGCGGAGAAGGGTTGCCCCGAAATAGCGACAATGGGGCAGGCGAAAAAGACCAAGCCCAAACGACGGGGGGGCCTCCAAAAAAAACAGCCCATCCGGGACGGTCCGGCGGAGAACCAACGGGTCAACCGtt                           >  1:501318/1‑151 (MQ=255)                                                                                                                         tCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAACGGTTc                          <  2:49559/149‑1 (MQ=255)                                                                                                                          ccGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGGGACGGCGGAAAAAACCAGCCCCAAGCGATGGGGGCGCATCAAAAAAATACAGCCCAACAGGATCGGGACCGCGAAGATCAATAGGATCAACAGGTCGcc                       >  2:530143/1‑151 (MQ=255)                                                                                                                            ggATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAACGGTTCGCCGt                     <  2:603706/151‑1 (MQ=255)                                                                                                                              aTAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAACGGTTCGCCGTAt                   <  1:73887/151‑1 (MQ=255)                                                                                                                                               ggAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAACGGTTCGCCGTATGACCATTTGTAGATCaa  <  1:770169/151‑1 (MQ=255)                                                                                                                                                gAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAACGGTTCGCCGTATGACCATTTGTAGATCaa  >  2:945598/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                 AGCGCCTGATGCATTTGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAACGGTTCGCCGTATGACCATTTGTAGATCAA  >  NC_000913/4302510‑4302802

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