BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,356,445	A→T	6.1%	intergenic (‑34/+25)	dtpC ← / ← cadA	dipeptide and tripeptide permease/lysine decarboxylase, acid‑inducible

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,356,445	0	A	T	6.1%	919.7 / 34.0	262	intergenic (‑34/+25)	dtpC/cadA	dipeptide and tripeptide permease/lysine decarboxylase, acid‑inducible
Reads supporting (aligned to +/- strand): ref base A (117/128); new base T (7/9); total (124/138)
Fisher's exact test for biased strand distribution p-value = 8.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.56e-01

Only 100 of 265 total aligned reads displayed.

TCACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGA  >  NC_000913/4356177‑4356692

gcACTGCAGTGGGGGTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGGGCTGAACAGGCTGATGGAATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCAc                                                                                                                                                                                                                                                   <  1:94500/274‑1 (MQ=255)

                                   agccagccGCCGAGAATAGGGGTAACGTAAACCAGAGAAACATATGCGCTGAACAGGCTGATGGCATGGTTATAATCAAAACCAAGCTGATGGGTGAGATATAGAATGAGTAAGGCACGCAGGCCGTAAAAGCTGAAGTACTCCCAGATTTGGATCGCCGCGGTATAGTATATCGCGCGCGGCTGTGGGGGTGTTTTCATGTGTTCTCCTAATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTttct                                                                                                                                                                                                                 <  2:2149/274‑1 (MQ=255)

                                       ggcgaggaaagaaaGGGGTAAAGGGNNCCGGANAAGCAAAAGCGATGGAAAGGTCGATGGATTGATTATCATCAAAACAAAGCTGATGGGGTAGTGAAGTGATGAGGAAGGAACGCATGCCGTAAATACTGAATTACTCCCAGATTTGGTTCGCCACGATATAGTATATCGCGCGCGGGTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCaa                                                                                                                                                                                                         <  3:15947/265‑1 (MQ=255)

                                             caggaAAAGTGGTAACGAAAACAAGAGAGGCAAATAGGCTTAAAAGGCTTAGGGAATGGTTATCAGAAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATAGAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTC‑TTCAATAc                                                                                                                                                                                                      <  4:119317/271‑1 (MQ=255)

                                                                           cATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCAc                                                                                                                                                                                                                                                   <  2:488265/200‑1 (MQ=255)

                                                                             tatGCGCTGAACGAGCTGATGGCATGGTTAGCATCAAAACCAAGCTGTTGGGTGAGATAGAGAATGAGTAAGGAACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCGCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGc                                                                                                                                                                  <  1:1254/279‑1 (MQ=255)

                                                                             tatGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATCGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAGGTGGCAAGCAACTTCCCTTTTACGAGCTAATTATTTTTTGCTTTCTTTTTTCAATTCATTAACGGTTTCTCGGCCATTCGCCCAACGGt                                                                                                                                                                      >  4:468111/1‑275 (MQ=255)

                                                                             tatGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTcttctt                                                                                                                                                                                                             >  2:238924/1‑236 (MQ=255)

                                                                                                       gTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAg                                                                                                                                            <  2:595056/275‑1 (MQ=255)

                                                                                                                                gTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACtt                                                                                                                                                                                                                                                 <  6:1447745/149‑1 (MQ=255)

                                                                                                                                gTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACtt                                                                                                                                                                                                                                                 >  8:1648244/1‑149 (MQ=255)

                                                                                                                                  gagaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTCTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTcc                                                                                                                                                                                                                                               <  7:1101141/149‑1 (MQ=255)

                                                                                                                                    gaTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCt                                                                                                                                                                                                                                             <  6:573092/149‑1 (MQ=255)

                                                                                                                                     aTAGAGAATGAGTAAGGCACGCATGCCGTAAGAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTGTCCTTATGAGCTAAAAAGGCAAGTGGCAAGCCACTTCCCtt                                                                                                                                                                                                                                            >  7:1408254/1‑149 (MQ=255)

                                                                                                                                     aTAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCtt                                                                                                                                                                                                                                            <  5:573139/149‑1 (MQ=255)

                                                                                                                                        gagaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTa                                                                                                                                                                                                                                         <  5:1262355/149‑1 (MQ=255)

                                                                                                                                         agaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCtttttt                                                                                                                                                                                                                                        <  5:441212/149‑5 (MQ=255)

                                                                                                                                         agaATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTAc                                                                                                                                                                                                                                        <  5:600691/149‑1 (MQ=255)

                                                                                                                                           aaTGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGa                                                                                                                                                                                                                                      >  6:1673373/1‑149 (MQ=255)

                                                                                                                                                   aGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTccc                                                                                                                                                                                                                                              <  8:1137029/133‑1 (MQ=255)

                                                                                                                                                    ggCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAattat                                                                                                                                                                                                                             >  5:862731/1‑149 (MQ=255)

                                                                                                                                                     gCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAACAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCTGCCATCATCCtt                                                                                                                                                                           >  3:28657/1‑197 (MQ=255)

                                                                                                                                                     gCACGCATGCCGTAAAAACTGAAGTACTCCCAGATGTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTg                                                                                                                                                                           <  1:28657/198‑1 (MQ=255)

                                                                                                                                                         gCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAtttttt                                                                                                                                                                                                                        <  5:316082/149‑1 (MQ=255)

                                                                                                                                                            tGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAtttttt                                                                                                                                                                                                                        >  7:1128264/1‑146 (MQ=255)

                                                                                                                                                            tGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAtttt                                                                                                                                                                                                                          >  5:1372783/1‑144 (MQ=255)

                                                                                                                                                              ccGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCttt                                                                                                                                                                                                                   <  8:55525/149‑1 (MQ=255)

                                                                                                                                                               cGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTttc                                                                                                                                                                                                                  >  7:1589954/1‑149 (MQ=255)

                                                                                                                                                                  aaaaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTttcttc                                                                                                                                                                                                               <  7:542873/149‑1 (MQ=255)

                                                                                                                                                                   aaaaCTGAAGTACTCCCAGGTTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCtttttt                                                                                                                                                                                                                                        <  5:1070376/123‑5 (MQ=39)

                                                                                                                                                                    aaaCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTcttctt                                                                                                                                                                                                             <  5:258430/149‑1 (MQ=255)

                                                                                                                                                                      aCTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTc                                                                                                                                                                                                           <  7:353232/149‑1 (MQ=255)

                                                                                                                                                                           aGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAc                                                                                                                                                                                                      <  5:381384/149‑1 (MQ=255)

                                                                                                                                                                                tCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTaa                                                                                                                                                                                                 >  5:1185046/1‑149 (MQ=255)

                                                                                                                                                                                tCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTaa                                                                                                                                                                                                 >  7:908985/1‑149 (MQ=255)

                                                                                                                                                                                tCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTaa                                                                                                                                                                                                 >  7:1209702/1‑149 (MQ=255)

                                                                                                                                                                                 cccAAATTTGGATCGCCACGATATGGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAAc                                                                                                                                                                                                >  5:1461480/1‑149 (MQ=255)

                                                                                                                                                                                  ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCttttt                                                                                                                                                                                                                                         >  7:291724/1‑104 (MQ=255)

                                                                                                                                                                                  ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCTTGCCACTTCCCttttt                                                                                                                                                                                                                                         <  8:291724/107‑4 (MQ=255)

                                                                                                                                                                                  ccAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACg                                                                                                                                                                                               >  7:164317/1‑149 (MQ=255)

                                                                                                                                                                                   cAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCtt                                                                                                                                                                                                                    >  6:258591/1‑127 (MQ=255)

                                                                                                                                                                                    aGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCttt                                                                                                                                                                                                            >  5:1309097/1‑134 (MQ=255)

                                                                                                                                                                                    aGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCttt                                                                                                                                                                                                            <  6:1309097/134‑1 (MQ=255)

                                                                                                                                                                                      aTTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTAAACGGtat                                                                                                                                                                                           >  5:1547794/1‑149 (MQ=255)

                                                                                                                                                                                         tGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGc                                                                                                                                                                                        <  7:1346833/149‑1 (MQ=255)

                                                                                                                                                                                            aTCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGc                                                                                                                                                                                     >  5:1605216/1‑149 (MQ=255)

                                                                                                                                                                                            aTCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGc                                                                                                                                                                                     <  6:862731/149‑1 (MQ=255)

                                                                                                                                                                                               gCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGAGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCtt                                                                                                                                                                                                                    >  6:1047762/1‑115 (MQ=255)

                                                                                                                                                                                                   cGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGc                                                                                                                                                                              <  8:908985/149‑1 (MQ=255)

                                                                                                                                                                                                      tataGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTg                                                                                                                                                                           <  6:747692/149‑1 (MQ=255)

                                                                                                                                                                                                        taGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGAc                                                                                                                                                                         >  5:944677/1‑149 (MQ=255)

                                                                                                                                                                                                           tataTCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATAcc                                                                                                                                                                                                     >  4:761552/1‑118 (MQ=255)

                                                                                                                                                                                                           tataTCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGt                                                                                                                                                                      >  7:1568354/1‑149 (MQ=255)

                                                                                                                                                                                                                cgcgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCgg                                                                                                                                                                                      >  7:920483/1‑128 (MQ=255)

                                                                                                                                                                                                                  cgcgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGAtt                                                                                                                       >  4:439597/1‑189 (MQ=255)

                                                                                                                                                                                                                    cgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAtttttt                                                                                                                                                                                                                        <  8:747728/90‑1 (MQ=255)

                                                                                                                                                                                                                    cgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAtttttt                                                                                                                                                                                                                        >  7:747728/1‑90 (MQ=255)

                                                                                                                                                                                                                    cgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGt                                                                                                                                                             <  7:210546/149‑1 (MQ=255)

                                                                                                                                                                                                                    cgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCaa                                                                                                                                              <  4:62444/164‑1 (MQ=255)

                                                                                                                                                                                                                    cgcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAg                                                                                                                                            <  4:518095/166‑1 (MQ=255)

                                                                                                                                                                                                                     gcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTAt                                                                                                                                                                    >  1:139110/1‑141 (MQ=255)

                                                                                                                                                                                                                     gcgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTAt                                                                                                                                                                    <  3:139110/141‑1 (MQ=255)

                                                                                                                                                                                                                      cgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCaa                                                                                                                                                                                                         <  4:93461/103‑1 (MQ=255)

                                                                                                                                                                                                                      cgGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGa                                                                                                                                                           <  7:448080/149‑1 (MQ=255)

                                                                                                                                                                                                                       ggCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGaa                                                                                                                                                          <  6:709369/149‑1 (MQ=255)

                                                                                                                                                                                                                       ggCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATGTTTTGCTTTCTTCTTTCAATACCTTAACAGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACg                            <  2:422415/275‑1 (MQ=255)

                                                                                                                                                                                                                        gCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAat                                                                                                                                                         >  5:674797/1‑149 (MQ=255)

                                                                                                                                                                                                                         cTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAgcgc                                                                                                                            <  2:779476/177‑1 (MQ=255)

                                                                                                                                                                                                                          tgtgAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAatat                                                                                                                                                       >  8:1125691/1‑149 (MQ=255)

                                                                                                                                                                                                                            tgAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCg                                                                                                                                                     <  5:904724/149‑1 (MQ=255)

                                                                                                                                                                                                                            tgAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCg                                                                                                                                                     <  6:479401/149‑1 (MQ=255)

                                                                                                                                                                                                                              aGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGt                                                                                                                                                   >  6:1594076/1‑149 (MQ=255)

                                                                                                                                                                                                                              aGGGTGTTTCCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGt                                                                                                                                                   >  8:217960/1‑149 (MQ=255)

                                                                                                                                                                                                                                 gtgtTTTCATGTGTTCTCCTTATGAGCAAAATAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGTGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTc                                                                                                                                                <  5:500429/149‑1 (MQ=255)

                                                                                                                                                                                                                                    ttttCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGc                                                                                                                                                                              >  8:668734/1‑116 (MQ=255)

                                                                                                                                                                                                                                    ttttCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGc                                                                                                                                                                              <  7:668734/116‑1 (MQ=255)

                                                                                                                                                                                                                                    ttttCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCaaa                                                                                                                                             <  8:122486/149‑1 (MQ=255)

                                                                                                                                                                                                                                     tttCATGTGTTTGCCTTATGACGAAAAAAGGGAAGGGGCAAGCCACTTCCCTTGTAAGGGCTAAGTATTTTTTTCTTTCTTCTATCAAATCCATAACGGTATAGCGGGCAGCAGCCTGACGGTATGCACCGTGAAGATCGGTTTCAAAGCCCGGATAGTGGGCGCCGATTTCACACAGCATCGGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATAtt               <  4:30578/274‑1 (MQ=255)

                                                                                                                                                                                                                                        cATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGccc                                                                                                                                         <  6:1185046/149‑1 (MQ=255)

                                                                                                                                                                                                                                              ttCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAg                                                                                                                                   >  5:843149/1‑149 (MQ=255)

                                                                                                                                                                                                                                              ttCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATa     >  2:440272/1‑275 (MQ=255)

                                                                                                                                                                                                                                              ttCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGATAATACGa  >  1:182995/1‑278 (MQ=255)

                                                                                                                                                                                                                                                  ccTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTttct                                                                                                                                                                                                                 <  2:732711/67‑1 (MQ=255)

                                                                                                                                                                                                                                                   cTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCgg                                                         >  4:156388/1‑218 (MQ=255)

                                                                                                                                                                                                                                                    ttATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTAtttttt                                                                                                                                                                                                                        >  6:33557/1‑58 (MQ=255)

                                                                                                                                                                                                                                                    ttATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTCGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAgcg                                                                                                                             >  5:183094/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       tGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCAtt                                                     <  2:156390/218‑1 (MQ=255)

                                                                                                                                                                                                                                                          gCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACg                                                                                                                                                                        >  8:1471974/1‑100 (MQ=255)

                                                                                                                                                                                                                                                          gCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCcg                                     <  4:227615/231‑1 (MQ=255)

                                                                                                                                                                                                                                                            aaaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGAc                                                                                                                                                                         >  5:1279370/1‑97 (MQ=255)

                                                                                                                                                                                                                                                             aaaaaGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTca                                                                                                                    >  7:256091/1‑149 (MQ=255)

                                                                                                                                                                                                                                                              aaaaGGGAAGTGCCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTcac                                                                                                                   <  7:61130/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                aaGGGAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTcacac                                                                                                                 >  8:1569907/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                  gggAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAg                                                                                                               >  6:519731/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                    gAAGTGGCTTGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAat                                                                                                                                                         >  6:388526/1‑105 (MQ=255)

                                                                                                                                                                                                                                                                    gAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCa                                                                                                             >  6:924822/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                       gTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCt                                                                                                          <  6:828966/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                         ggCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATa                                                                                                                                    <  6:850117/121‑1 (MQ=255)

                                                                                                                                                                                                                                                                          gCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGt                                                                                                                                  >  8:33336/1‑122 (MQ=255)

                                                                                                                                                                                                                                                                          gCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCa                                                                                                       >  6:1087549/1‑149 (MQ=255)

TCACTGCAGTGCGGTTGCCGAGCAGGCGGTCGGCAAGCCAGCCGCCGAGAATAGGGGTAACGTAAACCAGAGAAGCATATGCGCTGAACAGGCTGATGGCATGGTTATCATCAAAACCAAGCTGATGGGTGAGATAGAGAATGAGTAAGGCACGCATGCCGTAAAAACTGAAGTACTCCCAGATTTGGATCGCCACGATATAGTATATCGCGCGCGGCTGTGAGGGTGTTTTCATGTGTTCTCCTTATGAGCAAAAAAGGGAAGTGGCAAGCCACTTCCCTTGTACGAGCTAATTATTTTTTGCTTTCTTCTTTCAATACCTTAACGGTATAGCGGCCATCAGCCTGACGGTATGCACCGTGAATATCGGTTTCAAAGCCCGGATAGTGAGCGCCGATTTCACACAGCATCTGCAGGAACTCCAGAACCGGACGGCTTTCTTCGGTGATCATTTCACCCGGCATTACCAGAGGAACTCCCGGCGGGTACGGAAGGATCATATTGGCGTTAATACGA  >  NC_000913/4356177‑4356692

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