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

CTGATGCATTTGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAACGGTTCGCCGTATGACCATTTGTAGATCAAAAACA  >  NC_000913/4302515‑4302807                                                                                                                                                |                                                                                                                                                      cTGATGCATTTGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGccc                                                                                                                                                 >  1:1232776/1‑150 (MQ=255)  tGATGCATTTGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTCCCg                                                                                                                                                >  1:921419/1‑150 (MQ=255)   gATGCATTTGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTa                                                                                                                                              <  2:1097577/151‑1 (MQ=255)   gATGCATTTGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTa                                                                                                                                              <  2:928850/151‑1 (MQ=255)   gATGCATTTGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGccc                                                                                                                                                 >  1:1133119/1‑148 (MQ=255)     tGCATTTGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGGaa                                                                                                                                             >  1:626825/1‑150 (MQ=255)      gCATTTGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAAt                                                                                                                                           <  2:1266365/151‑1 (MQ=255)       cATTTGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTaaa                                                                                                                                            >  1:22428/1‑149 (MQ=255)        aTTTGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAAGGTTGcccccaaa                                                                                                                                            >  1:828856/1‑143 (MQ=255)        aTTTGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAAGGTTGCCCGTcaa                                                                                                                                            >  2:298225/1‑148 (MQ=255)        aTTTGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGCAAATa                                                                                                                                          >  1:887957/1‑150 (MQ=255)         tttGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGTAGGGTTTcccctcaa                                                                                                                                            >  2:66255/1‑142 (MQ=255)             cGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACa                                                                                                                                    <  2:505758/151‑1 (MQ=255)             cGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACa                                                                                                                                    <  1:611398/151‑1 (MQ=255)              gAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGGCCCTAAAAAga                                                                                                                                        >  1:78253/1‑145 (MQ=255)               aaaTGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGTAAATAGCGACaa                                                                                                                                   >  1:877045/1‑150 (MQ=255)                   ggCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGc                                                                                                                               <  2:1222107/150‑1 (MQ=255)                   ggCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGc                                                                                                                               <  1:1071570/150‑1 (MQ=255)                   ggCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCg                                                                                                                              <  2:815080/151‑1 (MQ=255)                    gCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAGCGACAATgg                                                                                                                                >  2:732354/1‑148 (MQ=255)                        gACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAAGGTTTCCCGTAAATAACGACAATGGCGACgg                                                                                                                          >  2:2396/1‑150 (MQ=255)                        gACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGACGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGc                                                                                                                         <  1:1139293/151‑1 (MQ=255)                           tGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAGCGACAATGGCGACGGCg                                                                                                                        >  1:1248448/1‑149 (MQ=255)                            gggAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGaaa                                                                                                                     <  2:610075/151‑1 (MQ=255)                            gggAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGaa                                                                                                                      <  2:1175856/150‑1 (MQ=255)                             ggAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGGAAATAGCGACAAAGGGGACGGCGaaa                                                                                                                     >  2:295373/1‑150 (MQ=255)                               aaaCCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGaaaaa                                                                                                                   >  2:189015/1‑150 (MQ=255)                                  ccACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGAcc                                                                                                               <  2:1046018/151‑1 (MQ=255)                                    acaGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAg                                                                                                             <  2:1256943/151‑1 (MQ=255)                                     caGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAGCGACAAAAGCGACGGCGaaaaaaa                                                                                                                 >  1:1041890/1‑144 (MQ=255)                                         aCGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGTAAATAGCGACAATGGCGACGGCGGAAAAAACCAGcccc                                                                                                         >  2:1131504/1‑150 (MQ=255)                                           gCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAg                                                                                                      <  1:696013/151‑1 (MQ=255)                                               aTTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGAt                                                                                                  <  2:685612/151‑1 (MQ=255)                                                 tAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATg                                                                                                 <  2:87598/150‑1 (MQ=255)                                                  aaCGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATgg                                                                                                <  1:841535/150‑1 (MQ=255)                                                               aGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGTAAATAGCGACAAAGGCGACGGCGAAAAAGACCAGCCCCCAGCGATGGGTGCGCATCaaaa                                                                                   >  2:2937/1‑150 (MQ=255)                                                                 gTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACa                                                                                <  1:1187106/151‑1 (MQ=255)                                                                 gTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACa                                                                                <  2:926857/151‑1 (MQ=255)                                                                            gCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGTTTGCCCCAAAAAAACGACAAATGCGGCGGCGAAAAAAACCAGCCCCCCACGAAGGGTGGGCATTAAAAAAAAGCAGccc                                                                        >  1:943968/1‑148 (MQ=255)                                                                            gCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATc                                                                     <  2:986898/151‑1 (MQ=255)                                                                              caacaaCGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAg                                                                   <  2:395798/151‑1 (MQ=255)                                                                                acaacGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGa                                                                 <  1:228483/151‑1 (MQ=255)                                                                                acaacGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCCTAAATAGCGACAATTGCGGCGGCGAAAAAAAACAGCCCCAAGCGATGGGGGGGCATCAAAAAAAAACAGCCCCTc                                                                     >  2:1254803/1‑147 (MQ=255)                                                                                   acGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCg                                                              <  2:541162/151‑1 (MQ=255)                                                                                    cGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAGCGGCAATGGGGGCCGCGAAAAAAACCAGCCCCCAACGATGGGTGGGCATCAAAAACATACAGCCCATCAGGATccg                                                             >  2:470935/1‑149 (MQ=255)                                                                                       tACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGTAAAAAACGACAATGGCGACGGCGAAAAAAAACAGCCCCCAGCGAAGGGGGGGCCTCAAAAACATGCAGCCCAACAAGATCGGGcc                                                          >  1:1087186/1‑151 (MQ=255)                                                                                        accccacCGCGTTAAGCGCACCACCACTTCCGGATAGGCGGGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCgg                                                             <  1:120506/144‑1 (MQ=255)                                                                                          accacAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGGAAATAGCGACAATGGCGACGGCGAAAAAAACCAGCCCCAAGCGATGGGGGCGCATCAAAAACAAGCAGCCCATCAGGATCGGTCCgg                                                        >  2:1120307/1‑150 (MQ=255)                                                                                           ccacAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAACCCCCAGCGATGGGGGCGCATCAAAAAAAAGCAGCCCATCAGGAACGGTCCgg                                                        >  2:340583/1‑149 (MQ=255)                                                                                            cacaGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTTCCCCCAAATAGCGACAATGGCGCCGGCGAAAAAAACCAGCCCCAACCGCTGGGGGCGCATCAAAAACATGCCCCCCATCAAGATCGGGcc                                                          >  1:380536/1‑146 (MQ=255)                                                                                            cacaGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGCAAATAACGACAATGGCGACGGCGAAAAAAACCAGCCCCAAGCGATGGGTGGGCATCAAAAAAAACCAGCCCATCAGGATCGGTCCGGCg                                                      >  1:327908/1‑150 (MQ=255)                                                                                            cacaGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGGAAAAAACGACAAAGGCGACGGCGAAAAAGACCAGCCCCAAACGATGGGGGCGCATCAAAAACAAGCAGCCCATCAGGATCGGTCCGGCgg                                                     >  1:1087487/1‑150 (MQ=255)                                                                                            cacaGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAAGGTTTCCCGGTAATAAACAAAATTGCGACGGCGCAAAAAACCAGCCCCACGCGAAGGGGGCGCAACAAAAACCATCAGCCCATCAGGATCGGGCCcgcg                                                      >  2:1190421/1‑150 (MQ=255)                                                                                             acaGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGCAAATAGCGACAATGGCGACGGCGAAAAAAACCCGCCCCAAGCGATGGGTGCGCATCAAAAACAAGCAGCCCATCAGGATCGGGCCgggg                                                      >  1:283842/1‑147 (MQ=255)                                                                                              caGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCgatg                                                   <  2:940729/151‑1 (MQ=255)                                                                                                   ttAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAACGACAATGGCGACGGCGAAAAAAACCAGCCCCCAGCGATGGGGGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGGTGATc                                                >  2:937912/1‑149 (MQ=255)                                                                                                     aaGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATc                                            <  1:816341/151‑1 (MQ=255)                                                                                                      aGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCg                                           <  2:1162276/151‑1 (MQ=255)                                                                                                      aGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGGTTGCCCGTAAATAACGACAATGGCGACGGCGAAAAAAAACAGCCCCCAGCGAAGGGTGCGCATCAAAAAAATGCCGCCCATCAAGATCGGTCCGGCgaag                                                   >  2:1011587/1‑141 (MQ=255)                                                                                                         gcACGACCACTTCCGGATAGTCGAGCATGGCGGGGGAGGGTTGCCCGGAAAAAACGGCAATGGCGACGGCGAAAAAAACCAGCCCCAAGCGATGGGTGCGCATCAAAAACAAGCAGCCCCTCAGGATCGGTCCGGCGATGATCaa                                              >  1:1159048/1‑145 (MQ=255)                                                                                                          cACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAAGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAAACCAGCCCCCAGCGATGGGGGCGCCTCAAAAAAATACAGCCCATCAGGATCGGTCCGGCGGCGATCAATCGGAt                                        >  2:709795/1‑150 (MQ=255)                                                                                                             gACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAAc                                    <  2:524557/151‑1 (MQ=255)                                                                                                             gACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAAc                                    <  2:996116/151‑1 (MQ=255)                                                                                                                           aGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAACGGTTCGCCGTATg                       <  1:1107868/150‑1 (MQ=255)                                                                                                                                     ggCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAACGGTTCGCCGTATGACCATTTGTa             <  2:650029/150‑1 (MQ=255)                                                                                                                                     ggCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAACGGTTCGCCGTATGACCATTTGTa             <  2:138643/150‑1 (MQ=255)                                                                                                                                          ggAAGGTTTGCCCGTAAATATGGACAATGGGAACGGCGAAAAAGACCAGCCCCAAGCGATGGGTACGGATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAACGGTTCGCGGTATGACCATTTGTAGATCaa       <  1:1138049/151‑1 (MQ=255)                                                                                                                                          ggAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAACGGTTCGCCGTATGACCATTTGTAGATCaa       <  2:351604/151‑1 (MQ=255)                                                                                                                                               gTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAACGGTTCGCCGTATGACCATTTGTAGATCAAAAACa  <  2:877049/151‑1 (MQ=255)                                                                                                                                                |                                                                                                                                                      CTGATGCATTTGCGAAATGGCACGACTGGGAAACCACAGCACGCCGATTAACGTCATCAGCAAGGTTGGATAGAGGCCAACAACGATACACCACAGCGTTAAGCGCACGACCACTTCCGGATAGTCGAGCATGGCGGGGAAGGTTTGCCCGTAAATAGCGACAATGGCGACGGCGAAAAAGACCAGCCCCAAGCGATGGGTGCGCATCAAAAACATGCAGCCCATCAGGATCGGTCCGGCGATGATCAATCGGATCAACGGTTCGCCGTATGACCATTTGTAGATCAAAAACA  >  NC_000913/4302515‑4302807

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