BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	812,771	C→G	5.9%	L168V (CTG→GTG) ‡	bioD →	dethiobiotin synthetase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	812,771	0	C	G	5.9%	346.7 / 9.5	101	L168V (CTG→GTG) ‡	bioD	dethiobiotin synthetase
Reads supporting (aligned to +/- strand): ref base C (51/44); new base G (2/4); total (53/48)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.55e-01

Only 100 of 103 total aligned reads displayed.

CTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACC  >  NC_000913/812622‑812917

cTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGgac                                                                                                                                                    <  2:1639394/150‑1 (MQ=255)

    cggcTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGgactgac                                                                                                                                                <  2:1479930/150‑1 (MQ=255)

       cTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGActct                                                                                                                                             <  1:561772/150‑1 (MQ=255)

          gTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGc                                                                                                                                          >  1:1192387/1‑150 (MQ=255)

           tttACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCg                                                                                                                                         <  1:1571191/150‑1 (MQ=255)

            ttACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCgg                                                                                                                                        >  2:1090706/1‑150 (MQ=255)

             tACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCggg                                                                                                                                       >  2:460530/1‑150 (MQ=255)

              aCGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGt                                                                                                                                      >  2:1349029/1‑150 (MQ=255)

               cgccgcTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGtt                                                                                                                                     >  1:1101710/1‑150 (MQ=255)

                gccgcTTTCTGACACTTTCACTTTTGCAGATTGGGTAACTCAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCAGAGTCAGTCCGGCGTGTTg                                                                                                                                    <  2:1087255/150‑1 (MQ=255)

                gccgcTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGGAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACCCGCCAGAGTCAGTCCGGCGTGTTg                                                                                                                                    <  1:78020/150‑1 (MQ=255)

                gccgcTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTg                                                                                                                                    <  1:8434/150‑1 (MQ=255)

                gccgcTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTg                                                                                                                                    <  2:1192387/150‑1 (MQ=255)

                  cgcTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTggg                                                                                                                                  >  1:201103/1‑150 (MQ=255)

                  cgcTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTggg                                                                                                                                  <  1:1062120/150‑1 (MQ=255)

                     tttCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGtag                                                                                                                               <  2:1494130/150‑3 (MQ=255)

                       tCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCg                                                                                                                             <  2:1207584/150‑1 (MQ=255)

                         tGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGaa                                                                                                                           >  2:395578/1‑150 (MQ=255)

                          gACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAAc                                                                                                                          <  1:1090706/150‑1 (MQ=255)

                           acacTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACg                                                                                                                         <  1:1549275/150‑1 (MQ=255)

                                ttcactttTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGtt                                                                                                                    >  1:349630/1‑150 (MQ=255)

                                  cactttTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTAc                                                                                                                  >  2:1629419/1‑150 (MQ=255)

                                  cactttTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTAc                                                                                                                  <  1:194927/150‑1 (MQ=255)

                                   actttTGCAGATTGGGTAACACAGGAACAACTGCCGGTGCTACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACg                                                                                                                 <  1:1659416/150‑1 (MQ=255)

                                   actttTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACg                                                                                                                 >  1:1436919/1‑150 (MQ=255)

                                    ctttTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGc                                                                                                                <  1:1134300/150‑1 (MQ=255)

                                     ttttGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGcc                                                                                                               <  2:1177004/150‑1 (MQ=255)

                                         gCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCg                                                                                                           >  1:1536836/1‑150 (MQ=255)

                                           aGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCAGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCggg                                                                                                         >  2:1200931/1‑150 (MQ=255)

                                            gATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGa                                                                                                        >  2:1266916/1‑150 (MQ=255)

                                                gggTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAAc                                                                                                    >  2:1261705/1‑150 (MQ=255)

                                                gcgTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAAc                                                                                                    >  1:1041069/3‑150 (MQ=255)

                                                  gTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGt                                                                                                  <  1:202306/150‑1 (MQ=255)

                                                  gTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGt                                                                                                  >  2:1328604/1‑150 (MQ=255)

                                                   tAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTc                                                                                                 <  2:1436919/150‑1 (MQ=255)

                                                   tAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTc                                                                                                 >  2:1558408/1‑150 (MQ=255)

                                                    aaCACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCa                                                                                                >  2:558534/1‑150 (MQ=255)

                                                     acacAGGAACAACTGCCGGTGATACTGGTAGTTGTTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCAc                                                                                               <  1:422438/150‑1 (MQ=255)

                                                     acacAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCAc                                                                                               <  2:1297826/150‑1 (MQ=255)

                                                     acacAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCAc                                                                                               <  1:921651/150‑1 (MQ=255)

                                                     acacAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCAc                                                                                               <  1:915711/150‑1 (MQ=255)

                                                     acacAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTCTATTAATCACGAAATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCAc                                                                                               >  2:1651722/1‑150 (MQ=255)

                                                      cacaGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACg                                                                                              <  2:187287/150‑1 (MQ=255)

                                                      cacaGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACg                                                                                              <  2:1019692/150‑1 (MQ=255)

                                                       acaGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGc                                                                                             >  2:1536371/1‑150 (MQ=255)

                                                        caGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCt                                                                                            <  2:235090/150‑1 (MQ=255)

                                                           gAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGaa                                                                                         <  1:1261705/150‑1 (MQ=255)

                                                            aacaacTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAat                                                                                        <  1:646643/150‑1 (MQ=255)

                                                                acTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGGCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAatatat                                                                                    <  2:457368/150‑1 (MQ=255)

                                                                   gCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGAt                                                                                                                       >  2:352012/1‑112 (MQ=255)

                                                                   gCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGAt                                                                                                                       <  1:352012/112‑1 (MQ=255)

                                                                     cGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAGCACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCa                                                                               <  2:762916/150‑1 (MQ=255)

                                                                           tACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCa                                                                         >  2:299131/1‑150 (MQ=255)

                                                                             cTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCAcc                                                                       <  1:393023/150‑1 (MQ=255)

                                                                              tGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCAccc                                                                      <  1:599303/150‑1 (MQ=255)

                                                                               ggTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCg                                                                     >  1:865298/1‑150 (MQ=255)

                                                                                   gTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATg                                                                 <  1:1349029/150‑1 (MQ=255)

                                                                                    ttGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGCCCACGCTCACCCGCATGa                                                                >  2:1389307/1‑150 (MQ=255)

                                                                                    ttGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGa                                                                >  2:30653/1‑150 (MQ=255)

                                                                                        tgtgAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTcc                                                            >  1:297573/1‑150 (MQ=255)

                                                                                           gAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCcgc                                                         >  1:772952/1‑150 (MQ=255)

                                                                                           gAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCcgc                                                         >  1:622000/1‑150 (MQ=255)

                                                                                              aCTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGcgcc                                                      <  1:395578/150‑1 (MQ=255)

                                                                                                  ggCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCgctg                                                  >  1:1374638/1‑150 (MQ=255)

                                                                                                   gCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCgctgc                                                 >  2:1024742/1‑150 (MQ=255)

                                                                                                    cTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAatat                                                                                      <  2:590449/112‑1 (MQ=255)

                                                                                                    cTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAatat                                                                                      >  1:590449/1‑112 (MQ=255)

                                                                                                           aaTCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGTGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGgagag                                         >  1:652999/1‑150 (MQ=255)

                                                                                                           aaTCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGgagag                                         <  2:349630/150‑1 (MQ=255)

                                                                                                               aCGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATcc                                     >  1:1438000/1‑150 (MQ=255)

                                                                                                                   gATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTg                                 <  2:806910/150‑1 (MQ=255)

                                                                                                                    aTGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTgg                                >  1:1407072/1‑150 (MQ=255)

                                                                                                                     tGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGc                               <  1:558534/150‑1 (MQ=255)

                                                                                                                      gTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCt                              <  1:1024742/150‑1 (MQ=255)

                                                                                                                        tGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTg                            >  1:1366255/1‑150 (MQ=255)

                                                                                                                          aCTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCa                          <  2:772952/150‑1 (MQ=255)

                                                                                                                            tGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTCCGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCGGCTGGGAGAGATCCCCTGGCTTGCAGa                        <  2:409490/150‑1 (MQ=255)

                                                                                                                            tGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTAACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGa                        <  2:1536836/150‑1 (MQ=255)

                                                                                                                            ggCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGa                        <  2:415769/149‑1 (MQ=255)

                                                                                                                             gCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGaa                       >  1:1615472/1‑150 (MQ=255)

                                                                                                                             gCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGaa                       >  1:272781/1‑150 (MQ=255)

                                                                                                                              cacaGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGaaa                      >  2:232426/1‑150 (MQ=255)

                                                                                                                              cacaGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGaaa                      >  2:1652533/1‑150 (MQ=255)

                                                                                                                              cacaGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGaaa                      >  2:1488372/1‑150 (MQ=255)

                                                                                                                                caGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAAt                    >  1:1151636/1‑150 (MQ=255)

                                                                                                                                  ggTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATcc                  >  1:192898/1‑150 (MQ=255)

                                                                                                                                  ggTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATcc                  >  1:1215448/1‑150 (MQ=255)

                                                                                                                                  ggTAATACAACACGCCGGACTGACTCTGCCGGGTAGGCTGGCGAACGATGTTACGACTCCCGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCAGCTGCTCGCAGACATCCCCTGGCATGCAGAAAATcc                  >  1:429715/1‑150 (MQ=255)

                                                                                                                                     aaTACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGa               >  2:171881/1‑150 (MQ=255)

                                                                                                                                      acccaaccCGCCAGAGTCAGTCCGGCGTGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATc                   <  2:507502/142‑1 (MQ=255)

                                                                                                                                      acccaaccCGCCAGAGTCAGTCCGGCGTGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATc                   >  1:507502/4‑145 (MQ=255)

                                                                                                                                      acccaaccCGCCAGAGTCAGTCCGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGaa              >  2:210535/4‑150 (MQ=255)

                                                                                                                                         caacaCGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAAt           >  1:132186/1‑150 (MQ=255)

                                                                                                                                            cacGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCg        <  1:1286381/150‑1 (MQ=255)

                                                                                                                                            cacGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCg        <  1:117438/150‑1 (MQ=255)

                                                                                                                                             acGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCgg       <  1:1558408/150‑1 (MQ=255)

                                                                                                                                             acGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCgg       >  2:656295/1‑150 (MQ=255)

                                                                                                                                              cGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGc      >  2:1554578/1‑150 (MQ=255)

                                                                                                                                              cGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGc      <  1:1536371/150‑1 (MQ=255)

                                                                                                                                                  ggACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAAcc  >  1:817344/1‑150 (MQ=255)

CTGGCGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACC  >  NC_000913/812622‑812917

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