Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,693,417	G→C	6.3%	L283V (CTG→GTG)	uidB ←	glucuronide transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,693,417	0	G	C	6.3%	231.0 / 5.4	87	L283V (CTG→GTG)	uidB	glucuronide transporter
Reads supporting (aligned to +/- strand):  ref base G (38/43);  new base C (4/2);  total (42/45)
Fisher's exact test for biased strand distribution p-value = 4.23e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.82e-01

CCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAGCGCACCGATGCACA  >  NC_000913/1693269‑1693566                                                                                                                                                     |                                                                                                                                                      ccAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGc                                                                                                                                                      >  2:312307/1‑150 (MQ=255) ccAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGc                                                                                                                                                      <  2:219030/150‑1 (MQ=255) ccAATTGAAGCGAAGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGc                                                                                                                                                      <  2:31153/150‑1 (MQ=255)       gAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGCCGGAAACCCAGAAGAACCGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTTCCGCTCCTCGCGACCATCCCCGGCACCAGCGGTGcc                                                                                                                                                >  2:542436/1‑150 (MQ=255)       gAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGcc                                                                                                                                                >  1:541761/1‑150 (MQ=255)       gAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTAGTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGcc                                                                                                                                                >  1:1205540/1‑150 (MQ=255)          gCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCGCTGGTGCCggg                                                                                                                                             <  2:463401/150‑3 (MQ=255)          gCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCGCTGGTGCCggg                                                                                                                                             <  1:626660/150‑3 (MQ=255)                gCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGAGCTTCGCGACCAGCCCCGGAACCACCGGTCCCGATGccccg                                                                                                                                       >  2:587479/1‑147 (MQ=255)                gCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGAACAACGCGACCATCCCCGGCACAGCCGGTGCCGATGCccca                                                                                                                                       >  2:592828/1‑150 (MQ=255)                 ccAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAg                                                                                                                                      >  1:1030774/1‑150 (MQ=255)                   aaCGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTa                                                                                                                                    >  2:1127873/1‑150 (MQ=255)                   aaCGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTa                                                                                                                                    >  2:1069374/1‑150 (MQ=255)                    aCGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTac                                                                                                                                   >  1:444688/1‑150 (MQ=255)                       caaccaacGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTaccaa                                                                                                                                >  2:701069/1‑150 (MQ=255)                        aaccaacGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCCCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGCCCATCCCCCGCACCAGCGGTGCCGAAGCCAAAGTaccaac                                                                                                                               >  2:819200/1‑150 (MQ=255)                             acGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGt                                                                                                                          >  2:698874/1‑150 (MQ=255)                               gCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGttt                                                                                                                        <  2:434229/150‑1 (MQ=255)                                ccACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGtttt                                                                                                                       <  2:15363/150‑1 (MQ=255)                                   ccGGCTGTGACAAGAGGTAAACCCAGAATAACAGCAGATACCCGCCGGTCCCCAGCAACGCGCCAATCAGGAAGGTTTTCTTTTTCCCGACCCGCGCGACCATCCCCGGCCCCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCa                                                                                                                    <  1:1154741/150‑1 (MQ=255)                                   ccGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCACCCAGGTTTTGcc                                                                                                                    >  1:861813/1‑149 (MQ=255)                                    cGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCAc                                                                                                                   <  2:919514/150‑1 (MQ=255)                                    cGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCAc                                                                                                                   <  1:701069/150‑1 (MQ=255)                                     ggCAGTGACCAGACGGAAACGCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCAcc                                                                                                                  <  1:10627/150‑1 (MQ=255)                                            aCCAGACGGAAGCCCAGAAGTACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCCGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTCCCCTCCAGGTTTTGCCCCAGtccgg                                                                                                           >  2:172654/1‑146 (MQ=255)                                            aCCAGACGGAAACCCATAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCTCCAATCAGCAACGTATTCTTTTTACCGATCCTCGCGACCATAACCGGCACCAGCGGTGCCGATGCCACAGAAACAACCACCTTAAACACCAGTACCa                                                                                                           >  2:1107833/1‑150 (MQ=255)                                            aCCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCa                                                                                                           >  1:957312/1‑150 (MQ=255)                                             ccAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCACCCAGGTTTTGCACCAGTACCAg                                                                                                          >  2:442397/1‑150 (MQ=255)                                             ccAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAg                                                                                                          <  2:1068044/150‑1 (MQ=255)                                              cAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGc                                                                                                         <  1:542116/150‑1 (MQ=255)                                                   ggAAACCCAGAAGAACAGCAGATACCCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGt                                                                                                    <  2:154540/150‑1 (MQ=255)                                                    taaaCCCCGAAGATCAGCAGATAACCGCTGGTTCCCAGCAAAGCGCCAATCTGGAAGGTATTCTTTTTACCGATCCGCGCGACCATCCCCGGCACCAGCGGTGCGGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTg                                                                                                   <  1:272167/149‑1 (MQ=255)                                                       aCCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAAc                                                                                                <  2:98108/150‑1 (MQ=255)                                                           agaagaACAGCAGTTAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGCCCATCCCCGGCACCAGCGGTGCCGCTGCCACAGTACCAACCAGGTTTTGCCCCAGTACGAGCACAGTGCACAGcc                                                                                            >  2:338703/1‑150 (MQ=255)                                                           agaagaACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGcc                                                                                            <  2:482430/150‑1 (MQ=255)                                                            taagaaCAGCAGATAACCGCAGGTTCCCAGCACAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGccc                                                                                           <  2:1107599/149‑1 (MQ=255)                                                              agaaCAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTAGTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCgg                                                                                         >  2:445307/1‑150 (MQ=255)                                                                 aCAGCATATCACCGCATGTTCCCAGCAAAGCGCCAACCAGGAAGGAATTCTTTTAACCGTGCCTCGCGACCATCCCCGGCACCCGCGGGGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTTAACAGCCCGGTAt                                                                                      <  2:419799/150‑1 (MQ=255)                                                                    gcagATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACCGCCCGGTATCAt                                                                                   >  1:557164/1‑150 (MQ=255)                                                                       gATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGCCCATCCCCGGCACCCGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAGCAGCCCGGTATCATTTa                                                                                >  2:817249/1‑150 (MQ=255)                                                                        aTAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTaa                                                                               <  1:1270721/150‑1 (MQ=255)                                                                         tacccGCAGCTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTAACCGATCCGCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCCACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAac                                                                              <  1:1122939/147‑1 (MQ=255)                                                                         tAACCGCAGGTTCCAAGAAAAGCGCCAATCAGGAAGGAATTCTTTTTACCGATCCTCGCTACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAac                                                                              <  1:997551/150‑1 (MQ=255)                                                                         aacccGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAac                                                                              <  1:393284/147‑1 (MQ=255)                                                                          aaCCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAaca                                                                             <  2:742623/150‑1 (MQ=255)                                                                           aCCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCCGCACAGTGAACAGCCCGGTATCATTTAacac                                                                            >  2:598317/1‑150 (MQ=255)                                                                           aCCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAaccc                                                                            >  1:183858/1‑148 (MQ=255)                                                                           aCCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAacac                                                                            <  1:940486/150‑1 (MQ=255)                                                                           aCCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGAGAACAGCCCGGTATCCTTAAacac                                                                            >  2:1481/1‑150 (MQ=255)                                                                            ccGCAGGTTCCCTGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAAcaca                                                                           <  1:339591/150‑1 (MQ=255)                                                                             cGCAGGTACACAGCAAACCTCAAATCAGGAAGGTATTCTTTTTATCGATCCTCTCTACCATCCCCGGCAACTGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACAt                                                                          <  1:784740/150‑1 (MQ=255)                                                                                aGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAgc                                                                       <  2:1180451/150‑1 (MQ=255)                                                                                      ccAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGt                                                                 <  2:1040578/150‑1 (MQ=255)                                                                                          cAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGGCGATGACACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTACCACATAGCGCACGTAgca                                                             >  2:211491/1‑148 (MQ=255)                                                                                          cAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGCACAGCCCGGTATCAATTAACACATAGCGCACGTAGaa                                                             >  2:151462/1‑150 (MQ=255)                                                                                           ataGCGCCATTCAGGTAGGTATTCTTTTTACGTTTACTCGCGACCATCCCCGTCACAAGGGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGaac                                                            <  1:1086645/148‑1 (MQ=255)                                                                                           aaaGCGCCAATCCGGAAGGTTTAATTTTTACCTAACCTCGCTACCATCCCCGGCACCAGCGGTGCCTATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGaac                                                            <  1:1091906/150‑1 (MQ=255)                                                                                            aaGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGaaca                                                           <  2:1314374/150‑1 (MQ=255)                                                                                               cgcCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAacg                                                        >  2:138667/1‑150 (MQ=255)                                                                                                  cAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAacgacg                                                     <  2:327342/150‑1 (MQ=255)                                                                                                       aGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGgcg                                                <  1:576562/150‑1 (MQ=255)                                                                                                        ggAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTACCACAAAGCGCACGTAGAACAACGCCGAGgcgc                                               >  2:818899/1‑150 (MQ=255)                                                                                                        ggAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGgcgc                                               <  1:1069374/150‑1 (MQ=255)                                                                                                            ggTATTCTTTTAACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGAc                                           <  1:1013003/150‑1 (MQ=255)                                                                                                              tATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCg                                         <  1:484917/150‑1 (MQ=255)                                                                                                                ttCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCa                                       <  2:673778/150‑1 (MQ=255)                                                                                                                  cTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCaaa                                     >  2:956215/1‑150 (MQ=255)                                                                                                                    ttttACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCAGATGCCACAGTACCAACCAGGATTTGCACCAGTACCAGAACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAACgg                                   >  1:670104/1‑150 (MQ=255)                                                                                                                              ccTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAg                         <  1:445307/150‑1 (MQ=255)                                                                                                                               cTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGc                        <  1:398603/150‑1 (MQ=255)                                                                                                                                    gACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGcacaca                   <  2:444688/150‑1 (MQ=255)                                                                                                                                        aTCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGTCGCTGACCGCAAAGGTCGAAATCAGCACACACAgc               <  1:592828/150‑1 (MQ=255)                                                                                                                                        aTCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAgc               <  1:587479/150‑1 (MQ=255)                                                                                                                                          ccccGGCACCCGCGGTGCCGCTGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGGATCATTTGACACCTAGCGCACGTAGGACTACGCCGCGGCGCTGACCGCAAAGGTCGAACTCAGCACCCACAgcgc             >  1:668257/1‑150 (MQ=255)                                                                                                                                          catcGGCACCACTGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAAAACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAATGTCGAAATCAGCACACACAccgc             >  2:504480/4‑150 (MQ=255)                                                                                                                                           cccGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCTCCAGTACCAGCACAGTGAACAGCCCGGTATCAATTAACAAATAGAGCGCGTAGAACAACGACGAGGCGCTGACCGCATAGGTCGAAATCAGCACACACAGCGCa            >  2:614639/1‑150 (MQ=255)                                                                                                                                           cccGGCACCAGCGGTGCCGATGCCACAGTACAAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAGCGCa            >  2:697569/1‑150 (MQ=255)                                                                                                                                           atcGGCACCGCTGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCGGCCCAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGGAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACTCAGCGCa            >  2:1025131/3‑150 (MQ=255)                                                                                                                                            ccGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACCCAgcgccc           >  2:567821/1‑148 (MQ=255)                                                                                                                                               gCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTTCCAGCACAGTGAACAGCCCGGTATCATTTTACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACAGACAGCGCACCGa        >  2:923011/1‑150 (MQ=255)                                                                                                                                                cACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCCCACACAGCGCACCgcg       >  1:579659/1‑148 (MQ=255)                                                                                                                                                  ccAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGCACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAGCGCACCGATGc     >  1:1032635/1‑150 (MQ=255)                                                                                                                                                  ccAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAGCGCACCGATGc     <  1:493429/150‑1 (MQ=255)                                                                                                                                                  ccAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAGCGCACCGATGc     <  1:745472/150‑1 (MQ=255)                                                                                                                                                   cAGCGGTGCCGATGCCACAATACCAACCAGGTTTTGCACCAGTACCAGCACAGTGCACCGCCCGGGTGCATTTAACACATAGCGCACGTAGAACGACGACGCGGCGCTGCCCGCAAAGGGCGCAATCCGCACACCCCGCGCACCGCTGcc    >  2:707044/1‑149 (MQ=255)                                                                                                                                                    aGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAGCGCACCGATGcac   >  2:311246/1‑150 (MQ=255)                                                                                                                                                    aGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAGCGCACCGATGcac   <  1:1127873/150‑1 (MQ=255)                                                                                                                                                     gCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACAAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAGCGCACCGATGcaca  >  2:76366/1‑150 (MQ=255)                                                                                                                                                     |                                                                                                                                                      CCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAGCGCACCGATGCACA  >  NC_000913/1693269‑1693566

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