Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,364,396	0	A	T	43.0%	126.5 / 40.8	86	P530P (CCA→CCT)	yhcD	putative outer membrane fimbrial subunit usher protein
Reads supporting (aligned to +/- strand):  ref base A (36/13);  new base T (0/37);  total (36/50)
Fisher's exact test for biased strand distribution p-value = 1.66e-13
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ACGGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGA  >  NC_000913/3364248‑3364538                                                                                                                                                     |                                                                                                                                               acgGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCAtt                                                                                                                                              >  1:182819/1‑151 (MQ=255)   ggTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTg                                                                                                                                            >  1:309827/1‑151 (MQ=255)    gTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTgg                                                                                                                                           >  2:10961/1‑151 (MQ=255)       aCAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAg                                                                                                                                        <  2:741129/151‑1 (MQ=255)        cAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGa                                                                                                                                       <  1:305187/151‑1 (MQ=255)         aaTGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGa                                                                                                                                       <  2:1351908/150‑1 (MQ=255)         aaTGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGa                                                                                                                                       <  1:1131388/150‑1 (MQ=255)          aTGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGc                                                                                                                                          >  2:907444/1‑146 (MQ=255)           tGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATc                                                                                                                                     >  1:787596/1‑150 (MQ=255)             aCTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCaa                                                                                                                                   >  1:497880/1‑150 (MQ=255)               tCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAac                                                                                                                                 >  1:70632/1‑150 (MQ=255)                cAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAac                                                                                                                                 <  1:817134/149‑1 (MQ=255)                   gCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGc                                                                                                                            >  2:606908/1‑151 (MQ=255)                     aaCCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGcc                                                                                                                           >  2:1378560/1‑150 (MQ=255)                              ataaacggaaaaaaaGCGTAAAAACACGCCACGGTTTCCGTCAACTGGCAGCACCAACTTATTGCTTCAGAAAATAAAGAAGACGATGGTGATCTGTTCTACGTCAACATCATTATTCCTTTTGGCAGATAAACCACCGCCACACTGtat                                                                                                                  <  2:1287750/138‑1 (MQ=255)                               ttAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATAc                                                                                                                >  2:397837/1‑151 (MQ=255)                                tAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACt                                                                                                               >  1:69950/1‑151 (MQ=255)                                aaatgggAACAAGCCTTTCAAACCCCCCACGGTTCCCGTAACCTGGAAGAATCAATTTAGTGCTTAAAAAAATAAGAAAAACGATGGTGTTCTTTCCTACTTCAACATAAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATAc                                                                                                                <  2:1421162/145‑1 (MQ=255)                                 aacggaaaCAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAAAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTc                                                                                                              <  1:132872/145‑1 (MQ=255)                                 aaCTGGAACAAAGCGTTCAAACCCGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAAAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGATCAAACACCGCCACACTGTATACt                                                                                                               <  1:285545/150‑1 (MQ=255)                                 aaCTGAAACAAAGCGTTCAACCACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAAAAAATAATGAAGACGATGGTGATCTTTTCTACGTCAACATCAGTATTCCTTTTGGCAGATCAAACACCGCCACACTGTATACTc                                                                                                              <  1:1356536/151‑1 (MQ=255)                                     ggAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGcc                                                                                                           >  1:890282/1‑150 (MQ=255)                                          aaggcttaaaaaCCCGCCCGGGTTCCCGTCAACGGGCAGCATAAACTTAGTTCTTAAAAAAATAATGAAGACGATGGTGTTCTTTTCTCCGTAAACATCAGTATTCCTTTTGGCAGATAAAACACCGCCACACTGTATACTCGCCATGACg                                                                                                     <  2:30167/142‑1 (MQ=255)                                          aaaGCTTTAAAAAACGCCACGGTTTCCGTAAACTGGCAGCTTCATCTTAGTCCTTAAAAAAATAATAAAGACGAGGGTGATCTTTTCAACGCAAAAATCAATATTCTTTTTGGCATACAAAACCCCGCCACACTGTATACTCGCCATGAc                                                                                                      <  1:1253435/150‑1 (MQ=255)                                          aaaGCGTTAAAACACGCCACGGTTCCCGTAAACTGGCAGCATCAACTTAGTGCTTCAAAAAATAAGAAAGACGTTGGTGTTCTGTTCTACGTAAACATCATTATTCCTTTTGGCAGATAAACCACCGCCACACTGTATACTCGCCATGAc                                                                                                      <  1:555498/150‑1 (MQ=255)                                              cGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATc                                                                                                  >  1:1037502/1‑150 (MQ=255)                                               gTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCAt                                                                                                >  1:838390/1‑151 (MQ=255)                                                  cAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATaa                                                                                              >  2:251374/1‑150 (MQ=255)                                                  aaaaCACGCCACGTTTCCCGTCAACGGCCAGCACAATTTTAGTGCTTAAAAAAAAAAGGAAGACGATGGGGATCTTTTTTCGGTCAACATCATTATTCATTTTGGCAAATAAAACACCGCCACACTGTATACTCGCCATGACGATc                                                                                                  <  2:69231/145‑1 (MQ=255)                                                   aaccccGCCACGTTTCCCGTCAACTGGCAGCATCAACTTAGTCTTTCAGAAAAAAAGGAAGACGATGGTGATTTGTTTTACGTAAACATCATTTTTCCTTTTGGAAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATa                                                                                               <  2:964279/143‑1 (MQ=255)                                                    aaCACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAAc                                                                                           >  2:910660/1‑151 (MQ=255)                                                       acGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAAcc                                                                                          >  2:1097365/1‑149 (MQ=255)                                                          ccACGTTTTCCTTCAACGGGCAGCATCAACTTAGTGCTTCAAAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATCCCTTTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACt                                                                                      <  1:316534/150‑1 (MQ=255)                                                          ccACGTTTTCCGTAAACTGGCACCATCAACTTAGTGCTTAAAAAAATAATGAAGACGATGGTGTTCTTTTCTACGTCAACACCAGTATCCCTTTTGGCAGACCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACt                                                                                      <  2:263122/150‑1 (MQ=255)                                                           cACGGTTTCCGACAACGGGCAGCATCAACTTAGTGCTTCAAAAAATAAGGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTa                                                                                     <  1:1069903/150‑1 (MQ=255)                                                              ggTTTCCGTCAACTGGCGGCATCATCTTAGTGCTTCAAAAAAAAATAAAGACGATGGTGTTCTGTTCTCCGTCAACATCATTTTTACTTTTGGCAGACCAACCGCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGt                                                                                 <  1:1204278/151‑1 (MQ=255)                                                              ggTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGt                                                                                 >  2:1230716/1‑151 (MQ=255)                                                              ggTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGt                                                                                 >  1:217052/1‑151 (MQ=255)                                                                tttCCGTAAACGGGAACCATCAATTTAGTGCTTAAAAAAATAATGAAGACGAGGGTGATCTGTTTTACTTCAACATCAGTATTCCTTTTGGAAAATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATg                                                                                   <  2:42117/147‑1 (MQ=255)                                                                ttcccGTCAACTGGCGGCATAAATTTAGTGTTTCAAAAAAAAATGAAGACGAGGGTGATCTGTTTTACGTCAACATCAGTATTCCTTTTGGAAAAACAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATg                                                                                   <  2:497868/144‑1 (MQ=255)                                                                 ttCCGACAGCTGGCAGCATAACTTTAGTCCTTAAAAAAATAAGAAAGGCGATGGTGATCTGTTCTCCGTAAACTCCAGTATCCCTTTTGGCAAATAAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTAc                                                                               <  2:614828/150‑1 (MQ=255)                                                                   ccgcaaccgggCAGCATAAATTTAGTTCTTAAAAAAAAAAGAAAGACGATGGTAATCTGTTCTACGTAAACACCATTATTCATTTTGGAAAACCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACt                                                                              <  2:903248/140‑1 (MQ=255)                                                                   ccgcaaacgggaggcaCCACCTTAGTGCTTAAAAAAAAAATGAAGAGGATGGTGATCTGTTCTCCGTAAACATAAGTATTCCTTTTGGCAAACCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTg                                                                             <  2:587005/137‑1 (MQ=255)                                                                   ccGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGAAGATAAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACt                                                                              <  1:483412/149‑1 (MQ=255)                                                                   ccGTCAACGGGCAGCATCAACTTAGTGCTTAAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCATTATTCCTTTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACt                                                                              <  2:1043684/149‑1 (MQ=255)                                                                   ccGTAAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGATCAAACGCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTgg                                                                            <  1:599365/151‑1 (MQ=255)                                                                    cGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGgt                                                                           >  2:421258/1‑151 (MQ=255)                                                                        aaCTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTc                                                                        >  1:1213162/1‑150 (MQ=255)                                                                              cagcaTCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATggg                                                                   >  2:447142/1‑149 (MQ=255)                                                                               agcaTCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAgt                                                                >  2:610361/1‑151 (MQ=255)                                                                                        ttAGTGCTTAAAAAAAAAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCTTTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAgt                                                                <  2:955189/142‑1 (MQ=255)                                                                                           gTGCTTCAAAAAATAATGAAGACGATGGTAATCTGTTTTACGTCAACATCAGTATTCCTTTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCagatga                                                    <  2:326311/151‑1 (MQ=255)                                                                                              cTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCagatg                                                     >  2:337956/1‑147 (MQ=255)                                                                                               ttCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATAAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCagatgaga                                                  <  1:357762/149‑1 (MQ=255)                                                                                               ttCAAAAAATAATAAAGACGAGGGTGATCTGTTCTACGTAAACATCAGTATTCCTTTTGGCAAATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGGGTCGTCTCa                                                         <  2:996294/142‑1 (MQ=255)                                                                                                  aaaaaaTAATGAAGACGATGGTGATCTGTTCTACGTAAACACCATTATTCCTTTTGGCAGAACAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGt                                               <  2:941945/147‑1 (MQ=255)                                                                                                  aaaaaaTAAGGAAGACGAGGGTGATCTGTTCTCCGTCAACATCAGTATTCCTTTTGGCAGATAAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCagatgagatg                                                <  1:1258295/146‑1 (MQ=255)                                                                                                    aaaaTAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCct                                            >  2:814850/1‑150 (MQ=255)                                                                                                      aataatGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCctact                                         >  1:1360297/1‑151 (MQ=255)                                                                                                       ataatGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCctact                                         >  2:1388389/1‑150 (MQ=255)                                                                                                        taatGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTAt                                       >  2:1293056/1‑151 (MQ=255)                                                                                                            gAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTAtgtg                                    >  1:1276549/1‑150 (MQ=255)                                                                                                            gAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTAtgtg                                    >  2:1279418/1‑150 (MQ=255)                                                                                                                aCGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATg                                <  1:898245/150‑1 (MQ=255)                                                                                                                aCGATGGTGATCTGTTCTACGACAACATCAGTTTTCCTTTTGGAAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATg                                <  2:762485/150‑1 (MQ=255)                                                                                                                 cGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCt                              >  2:354783/1‑151 (MQ=255)                                                                                                                 cGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCt                              >  1:706573/1‑151 (MQ=255)                                                                                                                  gATGGTTACGTGTCCTGCGCCAACTGCATTATCCCTTTTGGAAGATCAAACTCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCt                              <  2:1073062/150‑1 (MQ=255)                                                                                                                  gATGGTGATCTTTTCTACGTCAACATCAGTATTCCTTTTGGAAAATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGc                               <  2:493997/149‑1 (MQ=255)                                                                                                                     gttGATCTTTTTTCCTTCAACATCAGTTTCCCTTTTGGAAGATCAAACCCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGa                            <  1:501913/147‑1 (MQ=255)                                                                                                                     ggTGTTTTGTTCTCCGTAAACATCATTATTCTTTTTGGAAAACAAAACCCCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGaa                           <  2:165905/150‑1 (MQ=255)                                                                                                                     ggTGTTTTGTTCTACGTAAACATATGTTTTCCTTTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGaa                           <  1:907465/150‑1 (MQ=255)                                                                                                                      gTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACg                         >  1:50487/1‑151 (MQ=255)                                                                                                                             gTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCac                  >  1:678154/1‑151 (MQ=255)                                                                                                                              ttCTCCGTCAACATCAGTATTCAATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCac                  <  1:1279446/150‑1 (MQ=255)                                                                                                                              tccTACGTAAACATCAGTATTCCATTTGGCAGAAAAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCa                   <  2:228374/147‑1 (MQ=255)                                                                                                                                  aCGTCAACATCAGTATTCCTTTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCacg                 <  2:182812/147‑1 (MQ=255)                                                                                                                                      cAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGt         >  2:620267/1‑151 (MQ=255)                                                                                                                                       aaCATCAGTTTTCCTTTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGt         <  2:217046/150‑1 (MQ=255)                                                                                                                                       aaCATCAGTATTCCTTTTGGCAAATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAAc           <  1:832595/148‑1 (MQ=255)                                                                                                                                       aaCATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTg        >  1:89232/1‑151 (MQ=255)                                                                                                                                       aaCACCAGTATTCCTTTTGGAAGACCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTg        <  2:677129/151‑1 (MQ=255)                                                                                                                                             aGTATTCCTTTTGGAAGACCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGa  <  2:193217/151‑1 (MQ=255)                                                                                                                                             aGTATTCCATTTGGCAAAAAAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGaaa     <  1:1230745/148‑1 (MQ=255)                                                                                                                                             aGTATCCCTTTTGGCAAAAAAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGaaa     <  2:1301644/148‑1 (MQ=255)                                                                                                                                               ttttCCTTTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGaaa     <  1:1315202/144‑1 (MQ=255)                                                                                                                                                     |                                                                                                                                               ACGGTGACAATGACTCAAGCAACCTTAATATTAACTGGAACAAAGCGTTCAAACACGCCACGGTTTCCGTCAACTGGCAGCATCAACTTAGTGCTTCAGAAAATAATGAAGACGATGGTGATCTGTTCTACGTCAACATCAGTATTCCATTTGGCAGATCAAACACCGCCACACTGTATACTCGCCATGACGATCATAAAACCCACTATGGTACTGGTGTCATGGGAGTCGTCTCAGATGAGATGTCCTACTATGTGAATGCTGAACGAGATCACGACGAACGTGAAACGA  >  NC_000913/3364248‑3364538

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