BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,693,416	A→G	5.6%	L283P (CTG→CCG) ‡	uidB ←	glucuronide transporter

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,693,416	0	A	G	5.6%	373.9 / 11.3	146	L283P (CTG→CCG) ‡	uidB	glucuronide transporter
Reads supporting (aligned to +/- strand): ref base A (62/69); new base G (2/6); total (66/80)
Fisher's exact test for biased strand distribution p-value = 2.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.36e-01

Only 100 of 147 total aligned reads displayed.

GACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAGCGCACCGATGCA  >  NC_000913/1693267‑1693564

  ccAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGc                                                                                                                                                    <  1:385887/150‑1 (MQ=255)

  ccAATTGAAGCGAAGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGc                                                                                                                                                    <  2:1105357/150‑1 (MQ=255)

    aaTTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCgg                                                                                                                                                  <  1:1534379/150‑1 (MQ=255)

       tGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGc                                                                                                                                               >  1:927977/1‑150 (MQ=255)

        gAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGcc                                                                                                                                              >  2:622050/1‑150 (MQ=255)

        gAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGcc                                                                                                                                              <  1:1593512/150‑1 (MQ=255)

           gCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCGCTGGTGCCggg                                                                                                                                           <  2:1776694/150‑3 (MQ=255)

           gCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCGCTGGTGCCggg                                                                                                                                           <  1:593502/150‑3 (MQ=255)

           gCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCGCGGGTGCCggg                                                                                                                                           <  1:1890261/150‑3 (MQ=255)

           gCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGAt                                                                                                                                           <  2:1207634/150‑1 (MQ=255)

           gCGATGGCCAACGCAACCAACGCCACCGGCAGGGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCACTCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGAt                                                                                                                                           <  2:1737687/150‑1 (MQ=255)

            cGATGGCCAAGGCAACCTACGCCACCGGCAGTGACCAGACGGAAACCCAGTTGAACAGCTCATAACCGCCGGTTCCCACCAAAGCGCCAATCAGGAAGGTATTCTTTTTAGCGAACCTCGCGACCATCCCCCGCACCAGCGGTGCCGATg                                                                                                                                          <  2:1171667/150‑1 (MQ=255)

            cGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCGCTGGTGCCgggg                                                                                                                                          <  1:1976801/150‑4 (MQ=255)

            cGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCGCTGGTGCCgggg                                                                                                                                          <  1:1984545/150‑4 (MQ=255)

              aTGGCCAACGCAACCAACTCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCTGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTATACCGATCCTCGCGACCATCCCCGGCACCCGCGGTGCCGATGcc                                                                                                                                        <  1:318797/150‑1 (MQ=255)

                ggCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAACCCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCcac                                                                                                                                      >  2:860858/1‑150 (MQ=255)

                ggCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGCCCATCCCCGGCACCAGCGGTGCCGCTGCcac                                                                                                                                      >  1:1498132/1‑150 (MQ=255)

                 gCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCcaca                                                                                                                                     >  1:1896440/1‑150 (MQ=255)

                    aaCGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGCTCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTa                                                                                                                                  >  1:154724/1‑150 (MQ=255)

                     aCGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAATGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTac                                                                                                                                 <  1:1420705/150‑1 (MQ=255)

                       gCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTacca                                                                                                                               >  2:1868086/1‑150 (MQ=255)

                            caacGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAg                                                                                                                          <  2:1543959/150‑1 (MQ=255)

                            caacGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCACATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAg                                                                                                                          <  2:1952533/150‑1 (MQ=255)

                            caacGCCACCGGCAGTGACCAGACGCAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAg                                                                                                                          <  2:630083/150‑1 (MQ=255)

                             aacGCCACCGGAAGTGACCATACGCAAACCCAGAAGAACAGCTGATAACCGCAGGTTCCCAGAAAAGCGCAAATCAGTAAGGTATTCTTTTTACCTATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAgg                                                                                                                         <  2:1668199/150‑1 (MQ=255)

                              acGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCACCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGt                                                                                                                        <  1:1353446/150‑1 (MQ=255)

                                gCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGttt                                                                                                                      >  1:1923119/1‑150 (MQ=255)

                                gCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGttt                                                                                                                      >  2:751480/1‑150 (MQ=255)

                                gCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGttt                                                                                                                      >  2:687920/1‑150 (MQ=255)

                                gCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGttt                                                                                                                      >  2:572407/1‑150 (MQ=255)

                                 ccACCGGCAGTGACCAGACGGAAACCCAGCAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCGTCAGTACCAACCAGGtttt                                                                                                                     >  2:1027831/1‑150 (MQ=255)

                                 ccACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGtttt                                                                                                                     <  1:1418678/150‑1 (MQ=255)

                                 ccACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGtttt                                                                                                                     <  2:1207225/150‑1 (MQ=255)

                                  cACCGGCAGTGCCCAGAAGGAAACCCACCAGACCAGCCCATAACACCAGGTTCCCTCCAAAGCCCCAATCAGCACGGTATTCTTTTTACCGATCGGCGGTCCCTTCCCCCGCACCCCCGGGGCCGAAGCCACAGTACCAACCTCGTTTTg                                                                                                                    <  2:837063/150‑1 (MQ=255)

                                  cACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTg                                                                                                                    >  2:477614/1‑150 (MQ=255)

                                  cACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGAAACCCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCCGCGGTGCCGATGCCACAGTACCAACCAGGTTTTg                                                                                                                    <  2:639029/150‑1 (MQ=255)

                                   aCCGGCAGTGACCAGACGGAAACCCAGAATAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGTAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGc                                                                                                                   <  1:911568/150‑1 (MQ=255)

                                   aCCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGc                                                                                                                   <  2:1292066/150‑1 (MQ=255)

                                   aCCGGCAGTGACAAGACTGAAACCCAGAAGAACAGCAGATAACCGCAGTTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCTACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGc                                                                                                                   <  1:67069/150‑1 (MQ=255)

                                    ccGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCCCAGTACCAACCAGGTTTTGCa                                                                                                                  >  1:1767381/1‑150 (MQ=255)

                                    ccGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCa                                                                                                                  <  2:1268316/150‑1 (MQ=255)

                                        cAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTTTTCTTTATACCGATCCGCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAg                                                                                                              <  1:139333/150‑1 (MQ=255)

                                          gTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTa                                                                                                            >  1:460724/1‑150 (MQ=255)

                                            gACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTAcc                                                                                                          <  2:808534/150‑1 (MQ=255)

                                               cAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGATGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGc                                                                                                       <  2:502684/150‑1 (MQ=255)

                                               cAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGc                                                                                                       >  1:670317/1‑150 (MQ=255)

                                               cAGACGGAAACCCAGAAGAACAGCAGATAACCGCACGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGCCCAACCCCGGCACCAGCGGTGCCGATGCCCCAGTACCAACCAGGTTTTGCCCCAGTCCCAGc                                                                                                       >  1:978549/1‑150 (MQ=255)

                                                 ccagtaaaCCCAGAATAACATCAGATAACCGCAGGTACCCAGCAACCCTCCAATCAGGAAGGAATTCTTTTTACCGAGCCTCGCGACCATCACCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGcac                                                                                                     <  2:1262808/145‑1 (MQ=255)

                                                   cGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAg                                                                                                   >  1:1164375/1‑150 (MQ=255)

                                                   cGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAg                                                                                                   >  2:1666119/1‑150 (MQ=255)

                                                     gAAACCCAGAAGAACAGCATATACCCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTg                                                                                                 <  1:1842126/150‑1 (MQ=255)

                                                     gAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTg                                                                                                 <  1:730825/150‑1 (MQ=255)

                                                       aaCCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCCGCACAGTGaa                                                                                               >  1:782880/1‑150 (MQ=255)

                                                        aCCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAAc                                                                                              <  1:1787983/150‑1 (MQ=255)

                                                           cAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGc                                                                                           <  2:426011/150‑1 (MQ=255)

                                                           cAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACAAGCACAGTGAACAGc                                                                                           >  1:567219/1‑150 (MQ=255)

                                                             gaagaaCAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTAACGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCAACAGTACCAGCACAGTGAACAGccc                                                                                         <  2:1098030/150‑1 (MQ=255)

                                                                  aCAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTAt                                                                                    >  2:1529586/1‑150 (MQ=255)

                                                                     gcagATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCAt                                                                                 >  2:1356096/1‑150 (MQ=255)

                                                                     gcagATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCAt                                                                                 <  1:622050/150‑1 (MQ=255)

                                                                      cagATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCACCCAGGTTTTGCACCAGTACCAGCCCAGTGAACAGCCCGGTATCAtt                                                                                >  2:1597345/1‑150 (MQ=255)

                                                                      cagATAACCGCAGGTTCCCAGCAAAGCGCCAAGCAGGAAGGTATTCTTTTTACCGATCCTCGCGCCCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCAtt                                                                                <  2:252999/150‑1 (MQ=255)

                                                                       agATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCACGCGACCATCCCCGGCACCCGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCGCCAGTACCAGCACAGTGAACAGCCCGGTATCAgtt                                                                               >  2:630402/1‑150 (MQ=255)

                                                                          tAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAac                                                                            <  2:724326/150‑1 (MQ=255)

                                                                          aaaCCGCAGGTTCCAAGCAAAGCGCCAATCTGTATGTTTTTCTTTTTACGGATCCTCGCCACCATCCCCGTCACCAGGGGTGCGGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAac                                                                            <  1:1760082/149‑1 (MQ=255)

                                                                           aaCCGCAGGTTCCCAGCAAAGCGCCAATCAGGAATGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAaca                                                                           <  1:863644/150‑1 (MQ=255)

                                                                            aCCGCAGGTTCCCAGCAAAGCGCCAATCAGGATGGTATTCTTTTTACCGATCCTCGCGACAATCCCAGGCACCAGCGGTGCCGATGCCACAGTACAAACCAGGTTTTGCACCAGTACCAGCACAGTGAACCGCCCGGTATCATTTAacac                                                                          >  2:962101/1‑150 (MQ=255)

                                                                            aCCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAacac                                                                          <  2:971823/150‑1 (MQ=255)

                                                                            aCCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAacaa                                                                          >  2:1355537/1‑149 (MQ=255)

                                                                             ccGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAAcaca                                                                         >  1:561592/1‑150 (MQ=255)

                                                                              cGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTGGCGACCATCGACGGGAGCAGCGGTGCCGATGCCACATTACCAACCATTTTTTGCACCATTACCAGCACAGTGAACAGCCCGGTATCATTTAACACAt                                                                        <  1:1131818/150‑1 (MQ=255)

                                                                              cGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACAt                                                                        <  1:1733793/150‑1 (MQ=255)

                                                                                    ttCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTATACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCa                                                                  <  2:896221/150‑1 (MQ=255)

                                                                                          gCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGa                                                            <  2:949337/150‑1 (MQ=255)

                                                                                            aaaGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACAAATAGCGCACGTAGaac                                                          >  1:1925708/1‑150 (MQ=255)

                                                                                             aaGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGaaca                                                         >  1:1654849/1‑150 (MQ=255)

                                                                                               gcgcCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTGTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGaccaac                                                       >  1:975856/1‑150 (MQ=255)

                                                                                                cgcCAATCAGGAAGGTATTCTTTTTACCGATCCACGCGACCATCCCCGGAACCAGCGGTGACGATGCCACAGTACCAACCAGGTTTTGAACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAacg                                                      <  2:691887/150‑1 (MQ=255)

                                                                                                  ccAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGATCAacgac                                                    >  1:498099/1‑150 (MQ=255)

                                                                                                  ccAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAacgac                                                    >  2:865943/1‑150 (MQ=255)

                                                                                                      tCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAgg                                                <  1:992468/150‑1 (MQ=255)

                                                                                                       cAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGgc                                               <  1:647146/150‑1 (MQ=255)

                                                                                                           aaGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGACACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTg                                           >  2:672872/1‑150 (MQ=255)

                                                                                                                   cTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCCCCAGTACCAGCACAGTGCACCGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCaaa                                   >  2:342727/1‑150 (MQ=255)

                                                                                                                         aCCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGCACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGa                             >  1:1122477/1‑150 (MQ=255)

                                                                                                                         aCCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAAACAGGTTTTGCACCAGTACCAGCCCAGTGAACAGCCCGGTATCATATAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGa                             <  1:1282993/150‑1 (MQ=255)

                                                                                                                          ccGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGaa                            <  2:1686834/150‑1 (MQ=255)

                                                                                                                            gATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACACCGACGAGGCGCTGACCGCAAAGGTCGAAAt                          >  2:1993520/1‑150 (MQ=255)

                                                                                                                            gATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGGGAACAGCCCGGTATCATTTAACACATAGCGCACGTCGAACATAGACGAGGCGCTGACCGCAATGGTCGaagt                          >  2:342344/1‑148 (MQ=255)

                                                                                                                             aTCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATc                         <  2:1662653/150‑1 (MQ=255)

                                                                                                                               ccTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAg                       <  2:1840063/150‑1 (MQ=255)

                                                                                                                                  cgcgACCATCCCCGGCACCAGAGGTGCCAATGCCACAGTACCAACCAGGTTTTGAACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGAAGAACAACGACGAAGAGCTGCCCGCAAAGGTCGAAATCAGcac                    >  1:1281595/1‑150 (MQ=255)

                                                                                                                                         aTCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAgc             >  1:1522176/1‑150 (MQ=255)

                                                                                                                                           ccccGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAgcgc           <  2:558367/150‑1 (MQ=255)

                                                                                                                                            cccGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACCGCCCGGTATCATTTAACACATAGCGCACGGTGAACAACGCAGCGGCGCTGACCGCAAAGGTCGAAATCAGCAGACACAGAGCa          >  2:1257784/1‑150 (MQ=255)

                                                                                                                                            cccGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAGCGCa          >  2:1376187/1‑150 (MQ=255)

                                                                                                                                            atcGGCACCGCTGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAgcgcc          >  2:1989397/3‑149 (MQ=255)

                                                                                                                                                  aCCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAGCGCACCGagg    >  2:1515241/1‑148 (MQ=255)

                                                                                                                                                    cAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAGCGCACCGATGca  <  1:718455/150‑1 (MQ=255)

                                                                                                                                                    cAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAGCGCACCGATGca  <  2:433228/150‑1 (MQ=255)

GACCAATTGAAGCGATGGCCAACGCAACCAACGCCACCGGCAGTGACCAGACGGAAACCCAGAAGAACAGCAGATAACCGCAGGTTCCCAGCAAAGCGCCAATCAGGAAGGTATTCTTTTTACCGATCCTCGCGACCATCCCCGGCACCAGCGGTGCCGATGCCACAGTACCAACCAGGTTTTGCACCAGTACCAGCACAGTGAACAGCCCGGTATCATTTAACACATAGCGCACGTAGAACAACGACGAGGCGCTGACCGCAAAGGTCGAAATCAGCACACACAGCGCACCGATGCA  >  NC_000913/1693267‑1693564

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