BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	1,311,610	2 bp→AG	coding (522‑523/720 nt)	tonB →	membrane spanning protein in TonB‑ExbB‑ExbD transport complex

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,311,610	0	G	A	97.3%	1141.0 / 9.3	335	G174G (GGG→GGA)	tonB	membrane spanning protein in TonB‑ExbB‑ExbD transport complex
Reads supporting (aligned to +/- strand): ref base G (8/1); new base A (158/167); total (166/169)
Fisher's exact test for biased strand distribution p-value = 1.93e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
*	NC_000913	1,311,611	0	C	G	97.3%	1275.0 / 15.0	336	Q175E (CAG→GAG)	tonB	membrane spanning protein in TonB‑ExbB‑ExbD transport complex
Reads supporting (aligned to +/- strand): ref base C (8/1); new base G (158/168); total (167/169)
Fisher's exact test for biased strand distribution p-value = 1.90e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.48e-01

Only 100 of 336 total aligned reads displayed.

AAAATACGGCACCGGCACGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGG  >  NC_000913/1311465‑1311757

||

aaaaTACGGCACCGGCACGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGgt                                                                                                                                                <  3:312248/151‑1 (MQ=255)

aaaaTACGGCACCGGCACGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCAGTGAAGGAGAGgt                                                                                                                                                >  8:1583968/1‑151 (MQ=255)

        gcaccggcacGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGgttaaagtta                                                                                                                                        <  7:1016549/151‑1 (MQ=255)

         caccggcacGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGgttaaagtta                                                                                                                                        >  6:1010147/1‑150 (MQ=255)

          accggcacGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGgttaaagttaat                                                                                                                                      <  2:58264/151‑2 (MQ=255)

           ccggcacGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAAt                                                                                                                                     >  4:1308299/1‑151 (MQ=255)

             ggcacGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAAttt                                                                                                                                   <  3:1611043/151‑1 (MQ=255)

              gcacGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTg                                                                                                                                  <  6:1150758/151‑1 (MQ=255)

               cacGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGa                                                                                                                                 <  8:194271/151‑1 (MQ=255)

               aaagacTGACATGAATTACAGCAACTGCTGCAACCAGCAAGCCTGTTACCAGTGTGGCTTAAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCAGGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGa                                                                                                                                 <  3:1153308/146‑1 (MQ=255)

                actccTGAAATCAAGTACAGCAACGGCTTCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGAc                                                                                                                                <  2:706816/148‑1 (MQ=255)

                acGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGAc                                                                                                                                >  5:445856/1‑151 (MQ=255)

                acGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGAc                                                                                                                                >  1:1136210/1‑151 (MQ=255)

                acGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGAc                                                                                                                                <  7:428853/151‑1 (MQ=255)

                acGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGAc                                                                                                                                <  6:822363/151‑1 (MQ=255)

                  gCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGt                                                                                                                              >  4:1504561/1‑151 (MQ=255)

                  gCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGt                                                                                                                              >  7:523024/1‑151 (MQ=255)

                  gCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGt                                                                                                                              <  3:716231/151‑1 (MQ=255)

                  gCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGt                                                                                                                              >  3:402729/1‑151 (MQ=255)

                  gCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGt                                                                                                                              <  3:1498583/151‑1 (MQ=255)

                   ccTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGtt                                                                                                                             <  7:1649563/151‑1 (MQ=255)

                      gACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTAc                                                                                                                           <  5:112836/150‑1 (MQ=255)

                      gACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACg                                                                                                                          <  3:196995/151‑1 (MQ=255)

                       aCATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGc                                                                                                                         >  2:356064/1‑151 (MQ=255)

                        cATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGcc                                                                                                                        <  2:516262/151‑1 (MQ=255)

                        cATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGcc                                                                                                                        <  7:252533/151‑1 (MQ=255)

                           cAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCgg                                                                                                                      <  4:1157883/150‑1 (MQ=255)

                           cAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCgg                                                                                                                      <  7:763254/150‑1 (MQ=255)

                            aaGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATAGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGAt                                                                                                                    >  8:822900/1‑151 (MQ=255)

                              gTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATgg                                                                                                                  >  2:964893/1‑151 (MQ=255)

                                 cAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGAACACGCGCATTACGCCGTAATCAGCCGCAGTATCAGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATTGTcg                                                                                                               >  8:1042237/1‑151 (MQ=255)

                                     aaCGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTg                                                                                                           >  1:1308373/1‑151 (MQ=255)

                                     aaCGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTg                                                                                                           >  3:1342201/1‑151 (MQ=255)

                                     aaCGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTg                                                                                                           >  4:1591719/1‑151 (MQ=255)

                                        ggCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGAt                                                                                                        >  7:1363304/1‑151 (MQ=255)

                                        ggCTGCAACCAGCAACCCGGTAACCAGTCTGCCTTCACGACCACGCGCATTAAGCCGTAATCAGCCGCACTATCCGGCACGAGCACAGCCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGATACGCAGGATCGTCGCGTCGa                                                                                                         >  2:851198/1‑150 (MQ=255)

                                            gCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAAc                                                                                                     <  5:1199109/150‑1 (MQ=255)

                                            gCAACCACCAAGCCGCTTACCAGTGTGGCATCAGGACCACGCGCATTAAGCCGAAATCAGCCGCAGTATCCGGCACGAGCACAGGCAATGCGCATTGAAGGAGACCTAAAAGTTAAATTTGACGTTACGCCGGATGCTcgccaggaatacg                                                                                                    >  1:1474613/1‑141 (MQ=255)

                                              aaCCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTa                                                                                                  <  4:817094/151‑1 (MQ=255)

                                                 cAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACa                                                                                                >  4:185773/1‑150 (MQ=255)

                                                  aGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACaa                                                                                               >  7:1494377/1‑150 (MQ=255)

                                                     aaGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATcc                                                                                           >  5:57361/1‑151 (MQ=255)

                                                      aGCCGGTTACCAGTGTGGCTTCAGGAGCACGCGCATTAAGCCGTATTCAGCCGGAGTATCCGGCACGAGCACAGGCTTTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCct                                                                                          <  6:1690560/151‑1 (MQ=255)

                                                       gCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTAGCGTGGATAACGTACAAATCct                                                                                          >  6:1208123/1‑150 (MQ=255)

                                                          ggTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCa                                                                                      <  6:1346482/151‑1 (MQ=255)

                                                          ggTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCa                                                                                      <  3:1580975/151‑1 (MQ=255)

                                                               ccAGTGTGGCTTCAGGACCACGCGCATCAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCaa                                                                                 <  8:1474857/151‑1 (MQ=255)

                                                                cAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAg                                                                                <  8:1097892/151‑1 (MQ=255)

                                                                     tgGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGc                                                                           >  3:1111821/1‑151 (MQ=255)

                                                                     tgGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGc                                                                           >  2:1390469/1‑151 (MQ=255)

                                                                     gggCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTATCAGCCAAGCCTGc                                                                           <  3:661904/150‑1 (MQ=255)

                                                                      ggCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGc                                                                           <  3:524238/150‑1 (MQ=255)

                                                                       gCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGa                                                                         <  5:133242/151‑1 (MQ=255)

                                                                            aGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATg                                                                    >  3:1238781/1‑151 (MQ=255)

                                                                             ggACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGt                                                                   <  2:282133/151‑1 (MQ=255)

                                                                               aCCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGttt                                                                 >  5:751135/1‑151 (MQ=255)

                                                                                ccccgcgCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTg                                                                <  4:1549071/148‑1 (MQ=255)

                                                                                  aCGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTg                                                                >  4:1222950/1‑149 (MQ=255)

                                                                                  aCGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAg                                                              <  7:610018/151‑1 (MQ=255)

                                                                                  aCGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAg                                                              <  5:1450609/151‑1 (MQ=255)

                                                                                     cgcATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGt                                                           <  4:287197/151‑1 (MQ=255)

                                                                                     cgcATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGt                                                           <  3:929381/151‑1 (MQ=255)

                                                                                     cgcATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGt                                                           >  8:410858/1‑151 (MQ=255)

                                                                                        aTTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTg                                                          >  6:294242/1‑149 (MQ=255)

                                                                                           aaGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGt                                                      <  5:365738/150‑1 (MQ=255)

                                                                                           aaGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTg                                                     >  8:607242/1‑151 (MQ=255)

                                                                                           aaGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTATCAGCCAAGCCTGCGAACATGTTTGAGCGGGAGTTg                                                     >  2:422683/1‑151 (MQ=255)

                                                                                              ccGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGaa                                                   >  4:496884/1‑150 (MQ=255)

                                                                                              ccGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGaa                                                   >  2:1136877/1‑150 (MQ=255)

                                                                                                gTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGaaaaa                                                <  1:231749/151‑1 (MQ=255)

                                                                                                 tAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAat                                               >  8:357951/1‑151 (MQ=255)

                                                                                                 tAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAat                                               >  1:574623/1‑151 (MQ=255)

                                                                                                   aTCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAatgc                                             <  2:1601025/151‑1 (MQ=255)

                                                                                                         cGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAatgcgatgcg                                       <  8:1355168/151‑1 (MQ=255)

                                                                                                             gTATCGGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGGACAAATCCTCTCAGCCAAGCCTGCGAACATCTTGGAGCGTGAGGTGAAATATGCGATgcgc                                      >  6:525209/1‑148 (MQ=255)

                                                                                                               aTCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATg                                 >  6:1491903/1‑151 (MQ=255)

                                                                                                               aTCCGGCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGCTg                                 >  8:728677/1‑151 (MQ=255)

                                                                                                                    gCACGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGtt                            <  8:274571/151‑1 (MQ=255)

                                                                                                                    gCACGAGCACAGGCATTGCGCATTGAAGGAGAGCTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGtt                            <  4:244995/151‑1 (MQ=255)

                                                                                                                      aCGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTAt                          >  7:323445/1‑151 (MQ=255)

                                                                                                                       cGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTtgtg                         >  2:1011344/1‑151 (MQ=255)

                                                                                                                       cGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATg                         <  4:1026670/151‑1 (MQ=255)

                                                                                                                       cGAGCACAGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATg                         <  5:911638/151‑1 (MQ=255)

                                                                                                                           cacaGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGcc                     >  4:136152/1‑151 (MQ=255)

                                                                                                                            acaGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCg                    >  7:126554/1‑151 (MQ=255)

                                                                                                                            acaGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCg                    >  8:1685703/1‑151 (MQ=255)

                                                                                                                             caGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCgg                   >  5:1503881/1‑151 (MQ=255)

                                                                                                                              aGGCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCggg                  >  4:460390/1‑151 (MQ=255)

                                                                                                                               ggCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGt                 >  3:561521/1‑151 (MQ=255)

                                                                                                                                gCATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGt                 >  1:320109/1‑150 (MQ=255)

                                                                                                                                 cATTGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTa                >  4:980478/1‑150 (MQ=255)

                                                                                                                                   tttcgcATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGc             >  5:1065218/4‑151 (MQ=255)

                                                                                                                                    tGCGCATTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGcc            >  1:858788/1‑151 (MQ=255)

                                                                                                                                         aTTGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAACCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGc        >  8:1217322/1‑150 (MQ=255)

                                                                                                                                          ttGAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATCGCGTTATGAGCCGGGTAAGCCAgcc        >  1:1647108/1‑147 (MQ=255)

                                                                                                                                            gAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCGGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGt     <  4:677539/150‑1 (MQ=255)

                                                                                                                                            gAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGt     <  1:639817/150‑1 (MQ=255)

                                                                                                                                            gAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGt     <  2:1222223/150‑1 (MQ=255)

                                                                                                                                            gAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGt     <  5:191024/150‑1 (MQ=255)

                                                                                                                                            gAAGGAGAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTg    >  8:581430/1‑151 (MQ=255)

||

AAAATACGGCACCGGCACGCCTGACATCAAGTACAGCAACGGCTGCAACCAGCAAGCCGGTTACCAGTGTGGCTTCAGGACCACGCGCATTAAGCCGTAATCAGCCGCAGTATCCGGCACGAGCACAGGCATTGCGCATTGAAGGGCAGGTTAAAGTTAAATTTGACGTTACGCCGGATGGTCGCGTGGATAACGTACAAATCCTCTCAGCCAAGCCTGCGAACATGTTTGAGCGTGAGGTGAAAAATGCGATGCGCAGATGGCGTTATGAGCCGGGTAAGCCAGGCAGTGGG  >  NC_000913/1311465‑1311757

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