BRESEQ :: Evidence

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

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	812,773	0	G	C	5.9%	347.4 / 9.8	101	L168L (CTG→CTC) ‡	bioD	dethiobiotin synthetase
Reads supporting (aligned to +/- strand): ref base G (51/44); new base C (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 = 8.35e-01

Only 100 of 104 total aligned reads displayed.

CGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAA  >  NC_000913/812626‑812921

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:1062120/150‑1 (MQ=255)

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

                cTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACCCGCCAGAGTCAGTCCGGCGTGTTgtatt                                                                                                                                    <  2:1395987/150‑5 (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                                                                                                                      <  1:194927/150‑1 (MQ=255)

                              cactttTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTAc                                                                                                                      >  2:1629419/1‑150 (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                                                                                                   <  1:915711/150‑1 (MQ=255)

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

                                                 acacAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCAc                                                                                                   <  1:921651/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                                                                                          >  1:590449/1‑112 (MQ=255)

                                                                                                cTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAatat                                                                                          <  2:590449/112‑1 (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:272781/1‑150 (MQ=255)

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

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

                                                                                                                          cacaGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGaaa                          >  2:1652533/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                       >  1:507502/4‑145 (MQ=255)

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

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

                                                                                                                                    acaacaCGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGaaaa                >  1:382626/1‑150 (MQ=255)

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

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

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

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

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

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

                                                                                                                                                  tgactCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGaa  <  1:299131/150‑1 (MQ=255)

                                                                                                                                                  tgactCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGaa  <  1:1200931/150‑1 (MQ=255)

CGGCTGGTTTACGCCGCTTTCTGACACTTTCACTTTTGCAGATTGGGTAACACAGGAACAACTGCCGGTGATACTGGTAGTTGGTGTGAAACTCGGCTGTATTAATCACGCGATGTTGACTGCACAGGTAATACAACACGCCGGACTGACTCTGGCGGGTTGGGTGGCGAACGATGTTACGCCTCCGGGAAAACGTCACGCTGAATATATGACCACGCTCACCCGCATGATTCCCGCGCCGCTGCTGGGAGAGATCCCCTGGCTTGCAGAAAATCCAGAAAATGCGGCAACCGGAA  >  NC_000913/812626‑812921

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