BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	2,838,231	A→T	6.6%	intergenic (‑126/+23)	hydN ← / ← ascG	formate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit/asc operon transcriptional repressor; prpBC operon repressor

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	2,838,231	0	A	T	6.6%	517.6 / 21.5	151	intergenic (‑126/+23)	hydN/ascG	formate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit/asc operon transcriptional repressor; prpBC operon repressor
Reads supporting (aligned to +/- strand): ref base A (66/75); new base T (5/5); total (71/80)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.70e-02

Only 100 of 153 total aligned reads displayed.

ATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCACCGCTACGCCGCGCTCGTGTAACGCTTTCATCG  >  NC_000913/2837971‑2838501

                                                     ccTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTcc                                                                                                                                                                                                             >  2:628542/1‑275 (MQ=255)

                                                     ccTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAAAGCGCCGATACATTCCCCAATCATCCCAGTTTTTAATTTACTGTTATGTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGCTGTCATCGACACTATTGACGTTGACATGTGATAACAATGTTTCTCTCGAAGGAGCAATGAGCGAGTCGCGGCGGATCAGttttgt                                                                                                                                                                                                             >  4:514258/1‑273 (MQ=255)

                                                                                             aCCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgagtga                                                                                                                                                                                                                                  >  2:357787/1‑214 (MQ=255)

                                                                                                                tCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGa                                                                                                                                                                                                                                                                                <  5:55469/149‑1 (MQ=255)

                                                                                                                tCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGa                                                                                                                                                                                                                                                                                <  8:143765/149‑1 (MQ=255)

                                                                                                                 cAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGAc                                                                                                                                                                                                                                                                               <  8:381408/149‑1 (MQ=255)

                                                                                                                   gCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGa                                                                                                                                                                                                                                                                             <  5:397861/149‑1 (MQ=255)

                                                                                                                    cAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGAt                                                                                                                                                                                                                                                                            <  8:550721/149‑1 (MQ=255)

                                                                                                                    cAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGAt                                                                                                                                                                                                                                                                            <  8:419119/149‑1 (MQ=255)

                                                                                                                    cAATCATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTATTGCGCCGATACATTCACAAATCATGGCAGTTTTCAAATTACTGTTCTTTCAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACTAt                                                                                                                                                                                                                                                                            <  8:672112/149‑1 (MQ=255)

                                                                                                                     aaTGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATg                                                                                                                                                                                                                                                                           >  7:391730/1‑149 (MQ=255)

                                                                                                                     aaTGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATg                                                                                                                                                                                                                                                                           >  8:1995/1‑149 (MQ=255)

                                                                                                                       tgatgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGATGAc                                                                                                                                                                                                                                                                         <  7:754304/149‑1 (MQ=255)

                                                                                                                        gatgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTGATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACa                                                                                                                                                                                                                                                                        <  7:670521/149‑1 (MQ=255)

                                                                                                                         atgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAt                                                                                                                                                                                                                                                                       <  6:71337/149‑1 (MQ=255)

                                                                                                                          tgaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGATGACAtc                                                                                                                                                                                                                                                                      <  7:884442/149‑2 (MQ=255)

                                                                                                                           gaAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAtgt                                                                                                                                                                                                                                                                     <  5:263814/149‑1 (MQ=255)

                                                                                                                            aagcGGTTCATCAAATTCTCTAGCCATGACAGCTAATGCGCCGATACATTCACAAATTATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAtgtg                                                                                                                                                                                                                                                                    >  6:900351/4‑149 (MQ=255)

                                                                                                                            aaaCGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTGTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAtgtg                                                                                                                                                                                                                                                                    >  8:748703/1‑149 (MQ=255)

                                                                                                                             aaCGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtga                                                                                                                                                                                                                                                                   <  8:332095/149‑1 (MQ=255)

                                                                                                                                ggTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAtgtg                                                                                                                                                                                                                                                                    >  6:638083/1‑145 (MQ=255)

                                                                                                                                ggTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtgatga                                                                                                                                                                                                                                                                <  7:782912/149‑1 (MQ=255)

                                                                                                                                   tcatcaAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACaa                                                                                                                                                                                                                                                             <  6:240580/149‑1 (MQ=255)

                                                                                                                                     atcaAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATg                                                                                                                                                                                                                                                           >  5:1019339/1‑149 (MQ=255)

                                                                                                                                      tcaAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtga                                                                                                                                                                                                                                                                   >  6:795356/1‑140 (MQ=255)

                                                                                                                                      tcaAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtga                                                                                                                                                                                                                                                                   <  5:795356/140‑1 (MQ=255)

                                                                                                                                      tcaAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGt                                                                                                                                                                                                                                                          >  5:705863/1‑149 (MQ=255)

                                                                                                                                        aaaTTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGttt                                                                                                                                                                                                                                                        <  5:809872/149‑1 (MQ=255)

                                                                                                                                          aTTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTAt                                                                                                                                                                                                                                                      <  5:35732/149‑1 (MQ=255)

                                                                                                                                          aTTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTAt                                                                                                                                                                                                                                                      >  5:738989/1‑149 (MQ=255)

                                                                                                                                                cAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGaa                                                                                                                                                                                                                                                <  6:1065298/149‑1 (MQ=255)

                                                                                                                                                cAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGaa                                                                                                                                                                                                                                                <  6:13454/149‑1 (MQ=255)

                                                                                                                                                    aaTGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAg                                                                                                                                                                                                                                            <  5:505116/149‑1 (MQ=255)

                                                                                                                                                     aTGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTTTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGc                                                                                                                                                                                                                                           <  6:526839/149‑1 (MQ=255)

                                                                                                                                                     aTGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCAGCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGc                                                                                                                                                                                                                                           <  6:401786/149‑1 (MQ=255)

                                                                                                                                                       gACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCaa                                                                                                                                                                                                                                         <  6:904045/149‑1 (MQ=255)

                                                                                                                                                          aGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtga                                                                                                                                                                                                                                      <  5:774962/149‑1 (MQ=255)

                                                                                                                                                          aGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtga                                                                                                                                                                                                                                      <  5:850737/149‑1 (MQ=255)

                                                                                                                                                           gTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAtgtg                                                                                                                                                                                                                                                                    <  5:944782/118‑1 (MQ=255)

                                                                                                                                                           gTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgag                                                                                                                                                                                                                                     >  7:355737/1‑149 (MQ=255)

                                                                                                                                                           gTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgag                                                                                                                                                                                                                                     >  5:889956/1‑149 (MQ=255)

                                                                                                                                                            ttAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtgat                                                                                                                                                                                                                                                                  <  7:972706/119‑1 (MQ=255)

                                                                                                                                                                     cGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtgatga                                                                                                                                                                                                                                                                >  8:1009000/1‑112 (MQ=255)

                                                                                                                                                                      taTACATTCACAAATTATGCCAGTTTTGAATTTACTGTTATTTAAGGAAATTAATTTCTGGAATGGAGGAAAAACGATATAATCGACACTAGTGACGATGACATGTGATGAGAATGGTTATCGCGAAGGGGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGaa                                                                                            <  4:628542/274‑1 (MQ=255)

                                                                                                                                                                          cATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGAt                                                                                                                                                                                                                      <  6:293678/149‑1 (MQ=255)

                                                                                                                                                                           aTTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATc                                                                                                                                                                                                                     >  5:926279/1‑149 (MQ=255)

                                                                                                                                                                            ttCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCa                                                                                                                                                                                                                    >  5:386042/1‑149 (MQ=255)

                                                                                                                                                                              cacaAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGt                                                                                                                                                                                                                  <  8:721087/149‑1 (MQ=255)

                                                                                                                                                                                 aaaTCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGtttt                                                                                                                                                                                                               <  8:474610/149‑1 (MQ=255)

                                                                                                                                                                                  aaTCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTc                                                                                                                                                                                                              >  5:1026970/1‑149 (MQ=255)

                                                                                                                                                                                       tGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGgaga                                                                                                     >  2:232199/1‑249 (MQ=255)

                                                                                                                                                                                          cAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAg                                                                                                                                                                                                      >  8:209473/1‑149 (MQ=255)

                                                                                                                                                                                           aGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAgcgc                                                                                                 <  4:232199/249‑1 (MQ=255)

                                                                                                                                                                                             tttttAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATcg                                                                                                                                                                                                                                                    >  4:58609/1‑100 (MQ=255)

                                                                                                                                                                                                 tAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGaa                                                                                                                                                                                                                                                <  2:58609/100‑1 (MQ=255)

                                                                                                                                                                                                     ttACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGaa                                                             <  4:121750/275‑1 (MQ=255)

                                                                                                                                                                                                       aCTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTg                                                                                                                                                                                         <  7:233169/149‑1 (MQ=255)

                                                                                                                                                                                                         tGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAATAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGttt                                                                                                                                                                                                                                                        >  7:910568/1‑84 (MQ=255)

                                                                                                                                                                                                         tGTTATTTAAGGAAATTAATTTCTCTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGttt                                                                                                                                                                                                                                                        <  8:910568/84‑1 (MQ=255)

                                                                                                                                                                                                          gTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTgaga                                                                                                                                                                                      >  5:116833/1‑149 (MQ=255)

                                                                                                                                                                                                          gTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCCGCGGTgaga                                                                                                                                                                                      >  5:202918/1‑149 (MQ=255)

                                                                                                                                                                                                           ttATTTAAGGAAATGAATTTATGTAATGCAGGAAAAACGATGTCATCGACACTAGGGACGATGAGATGTGATGACAAAGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGaa                                                                                            <  1:148298/238‑1 (MQ=255)

                                                                                                                                                                                                                tAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATccc                                                                                                                                                                                <  5:518032/149‑1 (MQ=255)

                                                                                                                                                                                                                 aaGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATcccc                                                                                                                                                                               <  7:618565/149‑1 (MQ=255)

                                                                                                                                                                                                                  aGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCg                                                                                                                                                                              >  8:692677/1‑149 (MQ=255)

                                                                                                                                                                                                                     aaaTTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCa                                                                                                                                                                           >  7:732188/1‑149 (MQ=255)

                                                                                                                                                                                                                        ttaattTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCg                                                                                                                                                                        >  5:57670/1‑149 (MQ=255)

                                                                                                                                                                                                                          aattTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAg                                                                                                                                                                      <  8:381849/149‑1 (MQ=255)

                                                                                                                                                                                                                                tGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATTACTTACACCTTCTCTGGCa                                   >  2:355524/1‑274 (MQ=255)

                                                                                                                                                                                                                                  tAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATaaaaa                                                                                                                                                              <  8:355737/149‑1 (MQ=255)

                                                                                                                                                                                                                                      gCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAg                                                                                                                                                          >  6:376692/1‑149 (MQ=255)

                                                                                                                                                                                                                                         ggAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGGTTTCGGCGGTGAGAAATCACCGCCATCGAGCATAAAAATCAGCCg                                                                                                                                                       >  7:84378/1‑149 (MQ=255)

                                                                                                                                                                                                                                          gAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGt                                                                                                                                                      >  8:844118/1‑149 (MQ=255)

                                                                                                                                                                                                                                          gAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGt                                                                                                                                                      >  6:757175/1‑149 (MQ=255)

                                                                                                                                                                                                                                          gAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGt                                                                                                                                                      >  8:1039835/1‑149 (MQ=255)

                                                                                                                                                                                                                                            aaaaCGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGt                                                                                                                                                                                          <  5:606782/111‑1 (MQ=255)

                                                                                                                                                                                                                                            aaaaCGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCTCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGcc                                                                                                >  4:266511/1‑201 (MQ=255)

                                                                                                                                                                                                                                              tcacaTGTCATCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCaa                                                                                                                                                  >  8:153621/5‑149 (MQ=255)

                                                                                                                                                                                                                                              aaCGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCCCTGGCACCGCTACGCCgcgc                     >  2:21126/1‑274 (MQ=255)

                                                                                                                                                                                                                                                cGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGaa                                                                                            <  2:266511/201‑1 (MQ=255)

                                                                                                                                                                                                                                                  aTGTCATCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCaataat                                                                                                                                              >  8:816363/1‑149 (MQ=255)

                                                                                                                                                                                                                                                   tGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGa                                                                                                                                                                       >  6:522328/1‑123 (MQ=255)

                                                                                                                                                                                                                                                    gTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCACCGCTACGCCGCGCTCgtgt               >  4:427502/1‑274 (MQ=255)

                                                                                                                                                                                                                                                     tCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTc                                                                                                                                           <  6:421472/149‑1 (MQ=255)

                                                                                                                                                                                                                                                     tCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTc                                                                                                                                           <  5:171794/149‑1 (MQ=255)

                                                                                                                                                                                                                                                     tCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTc                                                                                                                                           <  7:1062416/149‑1 (MQ=255)

                                                                                                                                                                                                                                                       aTCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATc                                                                                                                                                                                  <  8:1027212/108‑1 (MQ=255)

                                                                                                                                                                                                                                                       aTCGTCACTAGTGTCGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATc                                                                                                                                                                                  >  7:1027212/1‑108 (MQ=255)

                                                                                                                                                                                                                                                       aTCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGCGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCt                                                                                                                                         >  8:270721/1‑149 (MQ=255)

                                                                                                                                                                                                                                                       aTCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCt                                                                                                                                         >  5:597400/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        tCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTg                                                                                                                                        >  5:333793/1‑149 (MQ=255)

                                                                                                                                                                                                                                                        tCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCACCGCTAcgcc                         >  4:46570/1‑260 (MQ=255)

                                                                                                                                                                                                                                                           acacTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAat                                                                                                                                     <  6:949213/149‑1 (MQ=255)

                                                                                                                                                                                                                                                               tAGTGACGATGACATGTGGTGACAATGTTTATCGCGAAGGGGCAGTGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAatcatc                                                                                                                                 >  8:624830/1‑149 (MQ=255)

                                                                                                                                                                                                                                                               tAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGGGGATCAGTATTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTTAatcaac                                                                                                                                 >  5:557082/1‑147 (MQ=255)

                                                                                                                                                                                                                                                               tAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGcc                                                                                                >  1:94429/1‑182 (MQ=255)

                                                                                                                                                                                                                                                                 gTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCACCGCTACGCCGCGCTCGTGTAACGCTTTCATCg  <  2:216414/274‑1 (MQ=255)

                                                                                                                                                                                                                                                                  tGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCa                                                                                                                              <  5:757175/149‑1 (MQ=255)

                                                                                                                                                                                                                                                                  tGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCa                                                                                                                              >  6:950496/1‑149 (MQ=255)

                                                                                                                                                                                                                                                                   gACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGaa                                                                                            <  3:94429/182‑1 (MQ=255)

ATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCACCGCTACGCCGCGCTCGTGTAACGCTTTCATCG  >  NC_000913/2837971‑2838501

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