BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	385,724	T→A	5.4%	W165R (TGG→AGG)	tauA →	taurine ABC transporter periplasmic binding protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	385,724	0	T	A	5.4%	832.0 / 25.5	235	W165R (TGG→AGG)	tauA	taurine ABC transporter periplasmic binding protein
Reads supporting (aligned to +/- strand): ref base T (101/121); new base A (6/7); total (107/128)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.04e-01

Only 100 of 237 total aligned reads displayed.

GCAACCTCGGTTCCAGCCCGTTAGCGGTTGCAGCCAGCCAACAGGTGCCGATTGAAGTCTTCTTGCTGGCGTCAAAACTGGGTAACTCCGAAGCGCTGGTGGTAAAGAAAACTATCAGCAAACCGGAAGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGC  >  NC_000913/385458‑385967

                                                    tGAAGTGTTCTTGATTGCGGCAAAACCGGGTAACTCCGAAGCGCTGGTGGTAAAGAAAACTATCAGCAAACCGGAAGATCTGATGGGCAAACGCATCGCCGTACCGTTTATCTCAACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGaa                                                                                                                                                                                          <  4:83324/274‑1 (MQ=255)

                                                           ttctgGCTGGCGTCAAAACTGGGTAACTCCCAAGCGGTGGTGGTAAAGAAAACTATCAGCAAACCGGAAGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGa                                                                                                                                                                                     <  4:260220/272‑1 (MQ=255)

                                                                                                                            ggAAGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGc                                                                                                                                                                                                                                               <  6:719113/149‑1 (MQ=255)

                                                                                                                              aaGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTCCAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCAGGCAGCga                                                                                                                                                                                                                                             <  5:936617/149‑2 (MQ=255)

                                                                                                                               aGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCAGGCAGCgat                                                                                                                                                                                                                                            <  6:1338893/149‑3 (MQ=255)

                                                                                                                                 aTCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGga                                                                                                                                                                                                                                          <  8:1018579/149‑1 (MQ=255)

                                                                                                                                 aTCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCAGGCAGCgataa                                                                                                                                                                                                                                          <  5:972300/149‑5 (MQ=255)

                                                                                                                                 aTCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCAGGCAGCgataa                                                                                                                                                                                                                                          <  7:1413635/149‑5 (MQ=255)

                                                                                                                                    tgatctgCAAACGCATCGCCGTACCTTTTATCTCCACCACCCACTACAGCCTGATGGCGGCACTGAACCAGTGGGGCATCAGCCCCGGTCAAGTTGAGATTTTGCACGTGCGGCCGCCCGAGATTATCGCTGCCTGTCAGAGGGGACat                                                                                                                                                                                                                                       <  8:740217/143‑1 (MQ=255)

                                                                                                                                    tGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGat                                                                                                                                                                                                                                       >  6:651568/1‑149 (MQ=255)

                                                                                                                                     gATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGATCCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACAGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTAGGGCGGTGGGGCGGGCCAACGCTGGACGTCGGGGTGGTGCGCaagg                                                                                                        >  1:10484/1‑273 (MQ=255)

                                                                                                                                       ttGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATAtt                                                                                                                                                                                                                                    <  6:1228942/149‑1 (MQ=255)

                                                                                                                                        tGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGata                                                                                                                                                                                                                                      >  5:718890/1‑146 (MQ=255)

                                                                                                                                         ggCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGa                                                                                                                                                                                                                                  <  5:1436722/149‑1 (MQ=255)

                                                                                                                                           cAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATg                                                                                                                                                                                                                                <  7:303051/149‑1 (MQ=255)

                                                                                                                                            aaaCGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATgg                                                                                                                                                                                                                               >  6:686197/1‑149 (MQ=255)

                                                                                                                                            aaaCGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATgg                                                                                                                                                                                                                               >  5:664523/1‑149 (MQ=255)

                                                                                                                                              aCGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTg                                                                                                                                                                                                                             <  7:65152/149‑1 (MQ=255)

                                                                                                                                              aCGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTg                                                                                                                                                                                                                             <  5:498932/149‑1 (MQ=255)

                                                                                                                                               cGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGc                                                                                                                                                                                                                            <  6:305210/149‑1 (MQ=255)

                                                                                                                                               cGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGc                                                                                                                                                                                                                            <  7:513175/149‑1 (MQ=255)

                                                                                                                                                 cATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCtt                                                                                                                                                                                                                          <  8:478850/149‑1 (MQ=255)

                                                                                                                                                      ccGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTc                                                                                                                                                                                                                     >  8:558201/1‑149 (MQ=255)

                                                                                                                                                        gTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTg                                                                                                                                                                                                                   <  8:373254/149‑1 (MQ=255)

                                                                                                                                                         tACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTgg                                                                                                                                                                                                                  >  8:808992/1‑149 (MQ=255)

                                                                                                                                                            cGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGAt                                                                                                                                                                                                                                 >  7:756430/1‑131 (MQ=255)

                                                                                                                                                             gTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCAc                                                                                                                                                                                                              >  7:222995/1‑149 (MQ=255)

                                                                                                                                                             gTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCAc                                                                                                                                                                                                              >  5:1150712/1‑149 (MQ=255)

                                                                                                                                                             gTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCCGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCAc                                                                                                                                                                                                              >  6:398467/1‑149 (MQ=255)

                                                                                                                                                                    tcCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGt                                                                                                                                                                                                       <  7:1436523/149‑1 (MQ=255)

                                                                                                                                                                    tcCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGt                                                                                                                                                                                                       <  7:1203187/149‑1 (MQ=255)

                                                                                                                                                                    tcCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGt                                                                                                                                                                                                       >  8:1411843/1‑149 (MQ=255)

                                                                                                                                                                          acccacTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGc                                                                                                                                                                                                 >  5:382313/1‑149 (MQ=255)

                                                                                                                                                                            ccccgaCAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGccc                                                                                                                                                                                               <  5:1452536/144‑1 (MQ=255)

                                                                                                                                                                             cacTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCt                                                                                                                                                                                              <  7:536828/149‑1 (MQ=255)

                                                                                                                                                                             cacTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCt                                                                                                                                                                                              <  6:1413536/149‑1 (MQ=255)

                                                                                                                                                                                    gCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGaaaaa                                                                                                                                                                                       <  5:45354/149‑1 (MQ=255)

                                                                                                                                                                                     gggggggggggaaCTGAAACGGNNGGGCANTAACCCCGGGCACGTGGAGAGTGGGAGCCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGGCTGGGCACAGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAgcg                                                        <  3:9865/263‑1 (MQ=255)

                                                                                                                                                                                      ctgctgGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGa                                                                                                                                                                                     <  6:1288324/149‑1 (MQ=255)

                                                                                                                                                                                       tgctgGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGAc                                                                                                                                                                                    <  7:455437/149‑1 (MQ=255)

                                                                                                                                                                                       tgctgGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGAc                                                                                                                                                                                    <  8:1119430/149‑1 (MQ=255)

                                                                                                                                                                                        gctgGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACg                                                                                                                                                                                   <  7:1206426/149‑1 (MQ=255)

                                                                                                                                                                                        gctgGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACg                                                                                                                                                                                   <  7:891792/149‑1 (MQ=255)

                                                                                                                                                                                             cggcACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAg                                                                                                                                                                              <  5:322887/149‑1 (MQ=255)

                                                                                                                                                                                                cACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGgtg                                                                                                                                                                           <  6:205421/149‑1 (MQ=255)

                                                                                                                                                                                                cACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGgtg                                                                                                                                                                           <  6:565084/149‑1 (MQ=255)

                                                                                                                                                                                                 aCTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTg                                                                                                                                                                                                                   <  5:510675/108‑1 (MQ=255)

                                                                                                                                                                                                 aCTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGa                                                                                                                                                                                           >  8:674302/1‑132 (MQ=255)

                                                                                                                                                                                                 aCTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGa                                                                                                                                                                                           <  7:674302/132‑1 (MQ=255)

                                                                                                                                                                                                 aCTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGgtgt                                                                                                                                                                          >  7:113940/1‑149 (MQ=255)

                                                                                                                                                                                                  cTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGtt                                                                                                                                                                         <  8:222995/149‑1 (MQ=255)

                                                                                                                                                                                                   tGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTg                                                                                                                                                                        >  7:1420199/1‑149 (MQ=255)

                                                                                                                                                                                                    gAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGa                                                                                                                                                                       >  6:745486/1‑149 (MQ=255)

                                                                                                                                                                                                    gAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGa                                                                                                                                                                       <  7:375677/149‑1 (MQ=255)

                                                                                                                                                                                                    gAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGa                                                                                                                                                                       >  5:640370/1‑149 (MQ=255)

                                                                                                                                                                                                          cTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGAtt                                                                                                                                                                 <  6:441595/149‑1 (MQ=255)

                                                                                                                                                                                                           tGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTc                                                                                                                                                                >  6:786914/1‑149 (MQ=255)

                                                                                                                                                                                                           tGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTc                                                                                                                                                                >  5:1319830/1‑149 (MQ=255)

                                                                                                                                                                                                            ggggCATTAAACCCGGCCAAGTGGAGATTGTGAACCTGCAGCCGCCAGCGATTATAGCTGCCTGACAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCt                                                                                                                                                               <  8:870322/149‑1 (MQ=255)

                                                                                                                                                                                                             gggCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACcggcg                                                                                                                                                                                                         <  5:79312/106‑1 (MQ=255)

                                                                                                                                                                                                             gggCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGgt                                                                                                                                                                            <  6:507358/135‑1 (MQ=255)

                                                                                                                                                                                                             gggCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGCTTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGAAGTCTGGGTGGTTCGCAAAGATTTTGCCGAGTAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGc                                              >  2:265265/1‑261 (MQ=255)

                                                                                                                                                                                                               gCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGaa                                                                                                                                                            <  7:164974/149‑1 (MQ=255)

                                                                                                                                                                                                                 aTTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACa                                                                                                                                                          <  8:686055/149‑1 (MQ=255)

                                                                                                                                                                                                                        ccGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTAt                                                                                                                                                                                                                        >  8:169550/1‑80 (MQ=255)

                                                                                                                                                                                                                        ccGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCggg                                                                                                                                                   >  7:73657/1‑149 (MQ=255)

                                                                                                                                                                                                                          gggCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGCAGACTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTGCGCTAAAAGCGCCCTCGATGCTCAGCAACCGTACATGGCTAACCCAGa                   >  1:74203/1‑275 (MQ=255)

                                                                                                                                                                                                                            gCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGt                                                                                                                                               <  7:238682/149‑1 (MQ=255)

                                                                                                                                                                                                                               aGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTggg                                                                                                                                            <  8:1344199/149‑1 (MQ=255)

                                                                                                                                                                                                                                 tGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGAc                                                                                                                                                                      >  6:169699/1‑121 (MQ=255)

                                                                                                                                                                                                                                 tGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGAc                                                                                                                                                                      <  5:169699/121‑1 (MQ=255)

                                                                                                                                                                                                                                     gaTTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGgcgcgc                                                                                                                                      >  7:95513/1‑149 (MQ=255)

                                                                                                                                                                                                                                     gaTTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGgcgcgc                                                                                                                                      <  6:978791/149‑1 (MQ=255)

                                                                                                                                                                                                                                        tgtgAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATAtt                                                                                                                                                                                                                                    >  2:375940/1‑52 (MQ=255)

                                                                                                                                                                                                                                         gtgAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGG‑AAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACg                                                                                                                                 <  6:1406307/149‑1 (MQ=255)

                                                                                                                                                                                                                                           gAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGc                                                                                                                                >  6:802382/1‑149 (MQ=255)

                                                                                                                                                                                                                                           gAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGc  >  2:41454/1‑275 (MQ=255)

                                                                                                                                                                                                                                            aaCCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATg                                                                                                                                                                                                                                <  4:375939/52‑1 (MQ=255)

                                                                                                                                                                                                                                             aCCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTg                                                                                                                              <  5:293822/149‑1 (MQ=255)

                                                                                                                                                                                                                                             aCCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATCGATGGTGCATATGTCTGCGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTg                                                                                                                              <  8:1185435/149‑1 (MQ=255)

                                                                                                                                                                                                                                                tGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGAc                                                                                                                           <  7:808992/149‑1 (MQ=255)

                                                                                                                                                                                                                                                         ccgcGATCATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGggtg                                                                                                                  >  5:1083996/1‑149 (MQ=255)

                                                                                                                                                                                                                                                          cgcgATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACgg                                                                                                                                                                                  <  6:1062893/84‑1 (MQ=255)

                                                                                                                                                                                                                                                          cgcgATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACgg                                                                                                                                                                                  >  5:1062893/1‑84 (MQ=255)

                                                                                                                                                                                                                                                          cgcgATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGggtgg                                                                                                                 >  6:965995/1‑149 (MQ=255)

                                                                                                                                                                                                                                                          cgcgATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGggtgg                                                                                                                 >  6:1145957/1‑149 (MQ=255)

                                                                                                                                                                                                                                                             gATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTgc                                                                                                              <  7:56950/149‑1 (MQ=255)

                                                                                                                                                                                                                                                              aTTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCaa                                                                                                                                                                               <  6:857549/83‑1 (MQ=255)

                                                                                                                                                                                                                                                              aTTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCaa                                                                                                                                                                               >  5:857549/1‑83 (MQ=255)

                                                                                                                                                                                                                                                               ttATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTgcgc                                                                                                            >  8:1279753/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                ccccGCTGCCAGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCa                                                                                                           >  5:574552/4‑149 (MQ=255)

                                                                                                                                                                                                                                                                 cccGCTGCCAGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGggtggt                                                                                                                <  7:1372757/141‑1 (MQ=255)

                                                                                                                                                                                                                                                                 aTCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCaa                                                                                                          >  7:381420/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                  tCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGAtt                                                                                                                                                                 <  5:1282857/93‑1 (MQ=255)

                                                                                                                                                                                                                                                                  tCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCaaa                                                                                                         <  5:254546/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                     cTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGAt                                                                                                      <  5:786914/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                       gCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGaa                                                                                                                                                            >  8:440854/1‑93 (MQ=255)

                                                                                                                                                                                                                                                                         cTGGCATCGGGGAGATATTGATGTTGCTTATGTCTGGGCACCGGCGGTTAACGCCTTGGAAAAAGACGGCAAGGTGGTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTTGTGCGCAAAGATTTTg                                                                                                  <  6:1112875/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                         cTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTg                                                                                                  <  8:1420199/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                         cTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTg                                                                                                  >  8:1362411/1‑149 (MQ=255)

GCAACCTCGGTTCCAGCCCGTTAGCGGTTGCAGCCAGCCAACAGGTGCCGATTGAAGTCTTCTTGCTGGCGTCAAAACTGGGTAACTCCGAAGCGCTGGTGGTAAAGAAAACTATCAGCAAACCGGAAGATCTGATTGGCAAACGCATCGCCGTACCGTTTATCTCCACCACCCACTACAGCCTGCTGGCGGCACTGAAACACTGGGGCATTAAACCCGGGCAAGTGGAGATTGTGAACCTGCAGCCGCCCGCGATTATCGCTGCCTGGCAGCGGGGAGATATTGATGGTGCTTATGTCTGGGCACCGGCGGTTAACGCCCTGGAAAAAGACGGCAAGGTGTTGACCGATTCTGAACAGGTCGGGCAGTGGGGCGCGCCAACGCTGGACGTCTGGGTGGTGCGCAAAGATTTTGCCGAGAAACATCCTGAGGTCGTGAAAGCGTTCGCTAAAAGCGCCATCGATGCTCAGCAACCGTACATTGCTAACCCAGACGTGTGGCTGAAACAGC  >  NC_000913/385458‑385967

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