BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	423,989	A→C	5.0%	N491T (AAC→ACC)	malZ →	maltodextrin glucosidase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	423,989	0	A	C	5.0%	633.7 / 19.1	175	N491T (AAC→ACC)	malZ	maltodextrin glucosidase
Reads supporting (aligned to +/- strand): ref base A (82/84); new base C (5/4); total (87/88)
Fisher's exact test for biased strand distribution p-value = 7.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.57e-01

Only 100 of 184 total aligned reads displayed.

TCAGCTCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAA  >  NC_000913/423843‑424136

tCAGCTCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACg                                                                                                                                                   >  2:1660623/1‑149 (MQ=255)

tCAGCTCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACg                                                                                                                                                   >  4:648659/1‑149 (MQ=255)

tCAGCTCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATAGCAAAAACg                                                                                                                                                   >  3:1161822/1‑149 (MQ=255)

 cAGCTCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGa                                                                                                                                                  <  4:730008/149‑1 (MQ=255)

 cAGCTCGACAGCCACGATACTGCGCGATTTAAAACGCTGCGCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGa                                                                                                                                                  <  1:1660623/149‑1 (MQ=255)

  aGCTCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGAt                                                                                                                                                 >  1:536557/1‑149 (MQ=255)

   gCTCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATc                                                                                                                                                <  1:472462/149‑1 (MQ=255)

    cTCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATcc                                                                                                                                               >  2:1246443/1‑149 (MQ=255)

          aGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTg                                                                                                                                         >  1:683559/1‑149 (MQ=255)

            ccACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAACGGATCGTTTTTGcc                                                                                                                                       <  2:896962/149‑1 (MQ=255)

            ccACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGcc                                                                                                                                       >  3:568963/1‑149 (MQ=255)

            ccACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGcc                                                                                                                                       >  1:362933/1‑149 (MQ=255)

            ccACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGcc                                                                                                                                       >  1:73319/1‑149 (MQ=255)

             cACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAACGGATCGTTTTTGcca                                                                                                                                      <  1:1417949/149‑2 (MQ=255)

              aCGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGt                                                                                                                                     <  4:1270072/149‑1 (MQ=255)

              aCGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGt                                                                                                                                     >  4:160480/1‑149 (MQ=255)

                gATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTaa                                                                                                                                   >  2:1840803/1‑149 (MQ=255)

                 aTACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTaaa                                                                                                                                  <  3:319974/149‑1 (MQ=255)

                      gcgcGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGtt                                                                                                                             <  1:617037/149‑1 (MQ=255)

                       cgcgATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTc                                                                                                                            >  2:1167227/1‑149 (MQ=255)

                        gcgATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGc                                                                                                                                        >  4:1440081/1‑136 (MQ=255)

                          gATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTcccc                                                                                                                         >  2:314319/1‑149 (MQ=255)

                          gATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTcccc                                                                                                                         <  1:1167227/149‑1 (MQ=255)

                          gATTTAAAACGCTGCCCAGTCGGGATATTGTGCGCTTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTATACCGTGCATTTATAACGGTGATGAAGTAGGACTGGATGGCAAATACAATCCGTTTTGCCGTAAACCGTTcccc                                                                                                                         >  4:1173591/1‑149 (MQ=255)

                            tttAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTg                                                                                                                       >  3:601589/1‑149 (MQ=255)

                            tttAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTg                                                                                                                       >  4:189894/1‑149 (MQ=255)

                                  aCGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGt                                                                                                                 <  3:189894/149‑1 (MQ=255)

                                   cGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTg                                                                                                                <  2:367468/149‑1 (MQ=255)

                                       gctCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGGGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGaaa                                                                                                            >  1:812860/1‑149 (MQ=255)

                                         tcggtcggGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAg                                                                                                          >  2:1286401/1‑149 (MQ=255)

                                           ggtcggGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCt                                                                                                                        >  4:1301801/1‑133 (MQ=255)

                                           ggtcggGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGGAGGTGGAAAAGct                                                                                                        <  4:1223345/149‑2 (MQ=255)

                                           ggtcggGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCa                                                                                                        <  1:1246443/149‑1 (MQ=255)

                                           ggtcggGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCa                                                                                                        <  2:1711097/149‑1 (MQ=255)

                                           ggtcggGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCa                                                                                                        <  2:362933/149‑1 (MQ=255)

                                           ggtcggGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATAGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCa                                                                                                        <  4:1161822/149‑1 (MQ=255)

                                            gtcggGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAg                                                                                                       >  2:1615178/1‑149 (MQ=255)

                                             tcggGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTccc                                                                                                                          <  4:1415103/129‑1 (MQ=255)

                                              cggGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGa                                                                                                     <  4:601589/149‑1 (MQ=255)

                                               gggATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGAt                                                                                                    >  1:1477332/1‑149 (MQ=255)

                                                 gATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATAc                                                                                                  >  4:542315/1‑149 (MQ=255)

                                                   tatTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAACGGATCGTTTTTGcca                                                                                                                                      >  1:756806/1‑110 (MQ=255)

                                                   tatTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGt                                                                                                                                             <  3:873260/104‑1 (MQ=255)

                                                   tatTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACgg                                                                                                >  4:634701/1‑149 (MQ=255)

                                                       gcgcgcCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGtt                                                                                            <  3:648659/149‑1 (MQ=255)

                                                           gcCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGa                                                                                                                                                  >  2:1770235/1‑91 (MQ=255)

                                                           gcCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTc                                                                                        <  1:1286401/149‑1 (MQ=255)

                                                           gcCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTc                                                                                        >  3:1434932/1‑149 (MQ=255)

                                                            ccTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTcg                                                                                       <  4:779264/149‑1 (MQ=255)

                                                                    tGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATa                                                                                                   >  1:1744384/1‑129 (MQ=255)

                                                                     ggCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTAcc                                                                              <  3:634701/149‑1 (MQ=255)

                                                                        ggtggtCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGt                                                                                                                 >  2:1148528/1‑111 (MQ=255)

                                                                           ggtCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGaa                                                                        <  4:568963/149‑1 (MQ=255)

                                                                            gtCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAAt                                                                       <  1:1840803/149‑1 (MQ=255)

                                                                                ggCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGAtt                                                                   <  1:314319/149‑1 (MQ=255)

                                                                                ggCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGAtt                                                                   <  3:542315/149‑1 (MQ=255)

                                                                                 gCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCg                                                                          >  3:1395395/1‑141 (MQ=255)

                                                                                     ttCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAAt                                                                       >  4:879785/1‑140 (MQ=255)

                                                                                       cACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGc                                                            <  2:248242/149‑1 (MQ=255)

                                                                                        aCCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCg                                                           >  3:1482520/1‑149 (MQ=255)

                                                                                          ctggcctggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTa                                                         >  1:1051968/1‑149 (MQ=255)

                                                                                           tggcctggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGt                                                                                                                 >  3:1464871/1‑92 (MQ=255)

                                                                                             gcctggTGTACCGTTCATTTATTACGGTGATTAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACTGTTCCCCTGGCAGGTGGTAAAGCAGGATACGGCGCTATTCGCGCTGTACCAGCGAACGATTGCGCTGCGTAAGa                                                      <  3:1173591/149‑1 (MQ=255)

                                                                                                tggTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGaaaa                                                   >  4:507005/1‑149 (MQ=255)

                                                                                                   tgtACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGt                                                                                             <  2:1455917/104‑1 (MQ=255)

                                                                                                   tgtACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTc                                                <  2:683559/149‑1 (MQ=255)

                                                                                                    gtACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCa                                               <  3:507005/149‑1 (MQ=255)

                                                                                                    gtACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCa                                               <  3:160480/149‑1 (MQ=255)

                                                                                                     tACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAg                                              <  4:1482520/149‑1 (MQ=255)

                                                                                                     tACCGTGCATTTATTACGGTGACGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAATAAAAGTCAg                                              <  1:1615178/149‑1 (MQ=255)

                                                                                                         gTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGgcgc                                          >  3:694592/1‑149 (MQ=255)

                                                                                                             aTTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCa                                                                                                        <  2:756205/83‑1 (MQ=255)

                                                                                                             aTTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTa                                        <  4:520171/147‑1 (MQ=255)

                                                                                                             aTTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTa                                        >  3:520171/1‑147 (MQ=255)

                                                                                                              tttATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTc                                                                                                                            <  4:553906/62‑1 (MQ=255)

                                                                                                              tttATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTc                                                                                                                            >  3:553906/1‑62 (MQ=255)

                                                                                                              tttATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTc                                                                                        >  2:1033299/1‑98 (MQ=255)

                                                                                                              tttATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTc                                                                                        <  1:1033299/98‑1 (MQ=255)

                                                                                                              tttATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGt                                     >  1:254759/1‑149 (MQ=255)

                                                                                                                tattaCGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCa                                   <  4:1161903/149‑1 (MQ=255)

                                                                                                                  ttaCGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTgc                                                                 <  4:61658/117‑1 (MQ=255)

                                                                                                                  ttaCGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTgc                                                                 >  3:61658/1‑117 (MQ=255)

                                                                                                                  ttaCGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAAACGTTCCCCTGGCAGGTGGaaa                                                                                                            <  2:54889/74‑1 (MQ=255)

                                                                                                                      ggTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcgg                             <  1:632841/149‑1 (MQ=255)

                                                                                                                      ggTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATggcgg                             >  2:460508/1‑149 (MQ=255)

                                                                                                                      ggTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATg                                 >  3:1054486/1‑145 (MQ=255)

                                                                                                                        tgatgaAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCt                           <  3:759418/149‑1 (MQ=255)

                                                                                                                                tAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTg                   >  3:1365482/1‑149 (MQ=255)

                                                                                                                                 aGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATgg                                >  4:1747562/1‑135 (MQ=255)

                                                                                                                                 aGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATgg                                <  3:1747562/135‑1 (MQ=255)

                                                                                                                                 aGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGc                  <  3:1270628/149‑1 (MQ=255)

                                                                                                                                  ggACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCt                 >  3:976290/1‑149 (MQ=255)

                                                                                                                                  ggACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCt                 <  2:766242/149‑1 (MQ=255)

                                                                                                                                   gACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTg                >  2:560105/1‑149 (MQ=255)

                                                                                                                                      tGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTAt             <  4:786419/149‑1 (MQ=255)

                                                                                                                                        gATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGaaaa                                                   >  3:1169144/1‑109 (MQ=255)

                                                                                                                                           ggCAAAACGGATCGTTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCgg         >  2:648451/1‑148 (MQ=255)

                                                                                                                                            gCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGaa       <  1:460508/149‑1 (MQ=255)

                                                                                                                                               aaaaCGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTa                                        <  2:504230/113‑1 (MQ=255)

                                                                                                                                                 aaCGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATaa  <  1:1815485/149‑1 (MQ=255)

TCAGCTCGACAGCCACGATACTGCGCGATTTAAAACGCTGCTCGGTCGGGATATTGCGCGCCTGCCGCTGGCGGTGGTCTGGCTGTTCACCTGGCCTGGTGTACCGTGCATTTATTACGGTGATGAAGTAGGACTGGATGGCAAAAACGATCCGTTTTGCCGTAAACCGTTCCCCTGGCAGGTGGAAAAGCAGGATACGGCGTTATTCGCGCTGTACCAGCGAATGATTGCGCTGCGTAAGAAAAGTCAGGCGCTACGTCATGGCGGCTGTCAGGTGCTGTATGCGGAAGATAA  >  NC_000913/423843‑424136

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 27 ≤` `ATCG/ATCG` `< 34 ≤` `ATCG/ATCG` `< 39 ≤` `ATCG/ATCG` `< 40 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`