BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,334,703	A→G	5.5%	intergenic (+22/‑206)	ispB → / → sfsB	octaprenyl diphosphate synthase/malPQ operon transcriptional activator

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,334,703	0	A	G	5.5%	684.1 / 16.0	204	intergenic (+22/‑206)	ispB/sfsB	octaprenyl diphosphate synthase/malPQ operon transcriptional activator
Reads supporting (aligned to +/- strand): ref base A (89/100); new base G (5/6); total (96/108)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.62e-01

Only 100 of 206 total aligned reads displayed.

TGCATCATGGCACACCAGAACAGGCACAGATGATCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAA  >  NC_000913/3334440‑3334967

aataacgaacacgacCGGAAGAGACACAGAGGGTCCGTAGCGCAATCGGACAGGGTAACGGTCGCCATCGTCGGGANCCGGTGCGGGAAGCAATGAACGCTTGTGGTTCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGGG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGa                                                                                                                                                                                                                                                                <  3:163307/261‑1 (MQ=255)

                     gggCCCAGATGACACGTACGGCCATCGAACAGGGTAACGGCCGCCATCTTCTGGAACCGGTGCTGGAAGCAAGGAAGGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGATAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCAGCCCGCGCAGC‑GCGCGGGATGGTTCCAGGAAGTTCCATaaa                                                                                                                                                                                                                                          <  4:483120/274‑1 (MQ=255)

                     aggcgaaggtgtcccGTACCGCCAGCGGACAGGGTAACGTGCGCCATCTTCTGGAACCGGTTGTGGAAGGAGTGAACGCTTGTGGAGCTCTTGAAGTGAAGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCTGCGTTACAGGTG‑GTCCCGGACACCCCATGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATaaa                                                                                                                                                                                                                                          <  4:29066/262‑1 (MQ=255)

                       gCACAGATGATCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGACGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgc‑gcgc                                                                                                                                                                                                                                                                    >  1:163308/1‑247 (MQ=255)

                                                                        tGGAACCGGTTCTGCAAGCAATGAACGCTGGTGGATCTCTTGAATTGACGCGTCAGCGTGCCGAGGAAGAAGCATACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGAACCCATCGGCCTCGCGAACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTtat                                                                                                                                                                                       <  4:167921/275‑1 (MQ=255)

                                                                        gtgAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTGGAATGGAAGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTtat                                                                                                                                                                                       <  2:586664/273‑1 (MQ=255)

                                                                        gggAAACGGGTTTGGAAGCAATGAAGGGTTGGGGAGATCTTGAAGCGACGCGTCAGGGTGCCGAGGAAGAAGCAGACAAAGCCAGCGCAGCGTGACAGGTG‑GTCCCGGATAACCCTTGGCGATAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGGTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCGATAAAAACACTTATTCAGCTCTATATCACCTGAAAAATAGGGATGTGCATCTTtat                                                                                                                                                                                       <  4:329757/274‑1 (MQ=255)

                                                                             acgggtagggAAGAAATTGACGTGTGTGGGTGTCTTGAATGGACGCGTCTGCGTGGCGAGGAATAAGTAGGCAAAGGCATCGCAGCGTTAGCGGTGACTCCTGG‑CACGCCATGGCGGGGAGCAGTCATCGGACTCGCGCAAATCGCTGGTCAACGCGATAGTTAATCCCCTCCCCTAATCCCGCGGAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTatatat                                                                                                                                                                                   <  4:112053/266‑1 (MQ=255)

                                                                                          cAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTcaca                                                                                                                                                                     <  2:368721/275‑1 (MQ=255)

                                                                                           aaTGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTc                                                                                                                                                                                                                                                          >  2:168010/1‑189 (MQ=255)

                                                                                                            ctctTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCa                                                                                                                                                                                                                                              >  2:646518/1‑184 (MQ=255)

                                                                                                                tGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAATCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATaaa                                                                                                                                                                                                                                          <  4:646516/184‑1 (MQ=255)

                                                                                                                     ggACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAg                                                                                                                                                                                                                                                                          <  5:820454/148‑1 (MQ=255)

                                                                                                                            tCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCT‑GCGCggg                                                                                                                                                                                                                                                                 >  8:626353/1‑149 (MQ=255)

                                                                                                                             cAGCGTGCCGAGGAGGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCT‑GCGCggg                                                                                                                                                                                                                                                                 <  6:148968/148‑1 (MQ=255)

                                                                                                                                cGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCT‑GCGCGGGATga                                                                                                                                                                                                                                                             <  8:78118/149‑2 (MQ=255)

                                                                                                                                 gTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATaaaaa                                                                                                                                                                                                                                        <  2:599983/169‑1 (MQ=255)

                                                                                                                                      gAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCa                                                                                                                                                                                                                                                        <  6:383705/148‑1 (MQ=255)

                                                                                                                                      gAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAg                                                                                                                                                                                                                                                       <  7:872842/149‑1 (MQ=255)

                                                                                                                                         gaagaaGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTaa                                                                                                                                                                                                                                                    >  7:792259/1‑149 (MQ=255)

                                                                                                                                         gaagaaGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTaa                                                                                                                                                                                                                                                    <  8:170060/149‑1 (MQ=255)

                                                                                                                                                  gACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATa                                                                                                                                                                                                                                            >  6:889742/1‑148 (MQ=255)

                                                                                                                                                   aCAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTcc                                                                                                                                                                                                                                                         <  5:647019/134‑1 (MQ=255)

                                                                                                                                                      aaGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGc                                                                                                                                                                                               <  4:588855/189‑1 (MQ=255)

                                                                                                                                                         ccATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCA‑‑‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTa                                                                                                                                                                                                                                  >  7:184883/1‑149 (MQ=255)

                                                                                                                                                           aTCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTa                                                                                                                                                                                                                                  >  8:490233/1‑149 (MQ=255)

                                                                                                                                                            tCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTAt                                                                                                                                                                                                                                 >  8:440617/1‑149 (MQ=255)

                                                                                                                                                             cGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCta                                                                                                                                                                                                                        >  1:33377/1‑157 (MQ=255)

                                                                                                                                                                 gCGTTACAGGTG‑CTCCCGNACACCCNTTGNCGAGAAGNACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAac                                                                                                                                                                                                                    <  3:33376/157‑1 (MQ=255)

                                                                                                                                                                 gCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGAATCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACt                                                                                                                                                                                                                                    >  4:84977/1‑141 (MQ=255)

                                                                                                                                                                  cGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGc                                                                                                                                                                                                                            >  6:490932/1‑148 (MQ=255)

                                                                                                                                                                   gTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgc‑gcgc                                                                                                                                                                                                                                                                    >  2:656931/1‑107 (MQ=255)

                                                                                                                                                                    ttACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGctc                                                                                                                                                                                                                          <  6:64788/148‑1 (MQ=255)

                                                                                                                                                                     tACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTAtt                                                                                                                                                                                                                                <  2:84978/141‑1 (MQ=255)

                                                                                                                                                                         ggTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATaa                                                                                                                                                                                                                     <  6:34227/148‑1 (MQ=255)

                                                                                                                                                                              cTCCCGGACACCCCTTGGCGAGAAGAACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGa                                                                                   <  2:492689/274‑1 (MQ=255)

                                                                                                                                                                               tCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTg                                                                                                                                                                                                               >  8:42784/1‑149 (MQ=255)

                                                                                                                                                                               tCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGC‑GCTGCGCGGGa                                                                                                                                                                                                                                                                <  7:481627/100‑1 (MQ=255)

                                                                                                                                                                                 ccGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATg                                                                                                                                                                                                                                                              >  5:676817/1‑100 (MQ=255)

                                                                                                                                                                                       acCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTc                                                                                                                                                                                                                                                          <  8:582848/98‑1 (MQ=255)

                                                                                                                                                                                       acCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTc                                                                                                                                                                                                                                                          >  7:582848/1‑98 (MQ=255)

                                                                                                                                                                                        ccccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCGCT‑GCGCGGGATg                                                                                                                                                                                                                                                              >  6:690784/1‑93 (MQ=255)

                                                                                                                                                                                          ccTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                                                                                                                                    <  8:195722/149‑1 (MQ=255)

                                                                                                                                                                                           cTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                                                                                                                                    <  6:557603/148‑1 (MQ=255)

                                                                                                                                                                                           cTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAt                                                                                                                                                                                                   <  8:83828/149‑1 (MQ=255)

                                                                                                                                                                                           cTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCt                                                                                                                                                                                           >  1:1501/1‑157 (MQ=255)

                                                                                                                                                                                             tGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgt                                                                                                                                                                                                 <  7:523569/149‑1 (MQ=255)

                                                                                                                                                                                             tGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtg                                                                                                                                                                                                  <  5:462126/148‑1 (MQ=255)

                                                                                                                                                                                               gCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgc‑gcg                                                                                                                                                                                                                                                                     <  2:118825/79‑1 (MQ=255)

                                                                                                                                                                                               gCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGc                                                                                                                                                                                               <  7:13452/149‑1 (MQ=255)

                                                                                                                                                                                                 gagaAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAg                                                                                                                                                                                                                                                                          <  2:154758/73‑1 (MQ=255)

                                                                                                                                                                                                   gaAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                                <  2:667254/144‑1 (MQ=255)

                                                                                                                                                                                                     aGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAgc‑gcgc                                                                                                                                                                                                                                                                    >  2:226709/1‑74 (MQ=255)

                                                                                                                                                                                                     aGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCAACCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCtt                                                                                                                                                                                          <  6:651944/148‑1 (MQ=255)

                                                                                                                                                                                                        aCTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAg                                                                                                                                                                                                       <  7:416409/132‑1 (MQ=255)

                                                                                                                                                                                                        aCTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGAGATGTGCATCTTtata                                                                                                                                                                                      >  7:573657/1‑149 (MQ=255)

                                                                                                                                                                                                         cTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGa                                                                                                                                                                                                                                                                <  4:226708/74‑1 (MQ=255)

                                                                                                                                                                                                         cTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAACTTAAGGa                                                                       >  4:655580/1‑259 (MQ=255)

                                                                                                                                                                                                          tCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                                >  6:618439/1‑137 (MQ=255)

                                                                                                                                                                                                            aTCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATAtt                                                                                                                                                                                  >  8:531183/1‑149 (MQ=255)

                                                                                                                                                                                                             tCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                                <  4:546490/134‑1 (MQ=255)

                                                                                                                                                                                                             tCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                                <  4:731081/134‑1 (MQ=255)

                                                                                                                                                                                                             tCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATAtt                                                                                                                                                                                  >  6:168541/1‑148 (MQ=255)

                                                                                                                                                                                                             tCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTc                                                                                                                                                                                 <  8:798121/149‑1 (MQ=255)

                                                                                                                                                                                                               ggCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTg                                                                                                                                                                               <  8:87641/149‑1 (MQ=255)

                                                                                                                                                                                                                 ccTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGTCTGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTTCATCTTTTNATATTCTGAATATTCACAATCTTAACCGGAACTTTTTAGAGTAATAGGCCATCTGGAGAATAGTGATTCTCGAANGAAGCAGTGTCTTGAATGAAAGCAAATAAANNAAGAAGGAAGGAAAAAGGGtaaaacattt                                               >  3:88229/1‑268 (MQ=255)

                                                                                                                                                                                                                 ccTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGa                                                                                                                                                                              <  5:719694/148‑1 (MQ=255)

                                                                                                                                                                                                                  cTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGaa                                                                                                                                                                             <  5:889742/148‑1 (MQ=255)

                                                                                                                                                                                                                     gcgcACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAatat                                                                                                                                                                          <  5:157714/148‑1 (MQ=255)

                                                                                                                                                                                                                      cgcACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAg                                                                                                                                                                                                                                                                          >  2:29526/1‑52 (MQ=255)

                                                                                                                                                                                                                          caTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                                <  2:743244/121‑1 (MQ=255)

                                                                                                                                                                                                                           aTCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGCTGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTCATCCTCGACAGAAGAAGTTTTCGGAATGAAAGTGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTATTTTCATTGAATgg                                     >  4:18998/1‑275 (MQ=255)

                                                                                                                                                                                                                            tCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAACTTCCATAAAAACACTTATTGAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCGCtctc                                                                                                                                                                  <  8:573657/149‑1 (MQ=255)

                                                                                                                                                                                                                             cGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGAGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCa                                   <  1:88230/275‑1 (MQ=255)

                                                                                                                                                                                                                               cTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCTGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTc                                                                                                                                                                        >  5:661034/1‑140 (MQ=255)

                                                                                                                                                                                                                                     aaCGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAg                                                                                                                                                           <  7:88329/147‑1 (MQ=255)

                                                                                                                                                                                                                                      aCGCGATCGTTAATCCCCTCCCCCCATCCCGCGCGCT‑GCGCggg                                                                                                                                                                                                                                                                 <  6:284221/44‑1 (MQ=18)

                                                                                                                                                                                                                                      aCGCGATCGTTAATCCCCTCCCCCCATCCCGCGCGCT‑GCGCggg                                                                                                                                                                                                                                                                 >  5:284221/1‑44 (MQ=18)

                                                                                                                                                                                                                                        gcgATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTa                                                                                                                                                              <  6:739974/141‑1 (MQ=255)

                                                                                                                                                                                                                                            tCGTTAATACCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACttttt                                                                                                                                                  >  7:635756/1‑149 (MQ=255)

                                                                                                                                                                                                                                             cGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTa                                                                                                                                                 <  8:403580/149‑1 (MQ=255)

                                                                                                                                                                                                                                             cGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTa                                                                                                                                                 <  8:635756/149‑1 (MQ=255)

                                                                                                                                                                                                                                               ttAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTaga                                                                                                                                               <  7:820312/149‑1 (MQ=255)

                                                                                                                                                                                                                                               ttAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTag                                                                                                                                                >  5:34046/1‑148 (MQ=255)

                                                                                                                                                                                                                                                tAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTagag                                                                                                                                              <  7:312982/149‑1 (MQ=255)

                                                                                                                                                                                                                                                   tcccctccccTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCaa                                                                                                                                           >  8:279176/1‑149 (MQ=255)

                                                                                                                                                                                                                                                     ccctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAAt                                                                                                                                          >  5:602340/1‑148 (MQ=255)

                                                                                                                                                                                                                                                      cctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATa                                                                                                                                         >  5:413044/1‑148 (MQ=255)

                                                                                                                                                                                                                                                       ctcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCttt                                                                                                                                                               >  7:243461/1‑125 (MQ=255)

                                                                                                                                                                                                                                                        tcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATCGGGACGTACATCTTTATATCTTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAATGCTTCTGAATGAACGCGAGTAACTTACGGAGTGAGGAAAAATGGCACATAATTTCACTGACTGGCCTCCCGCTGCCGTAATTACGGGTTTtcg        >  4:375778/1‑275 (MQ=255)

                                                                                                                                                                                                                                                        tcccctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCa                                                                                                                                            >  2:385288/1‑143 (MQ=255)

                                                                                                                                                                                                                                                            ctcATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACGGGAACTTTTTAGAGCAATAg                                                                                                                                        <  4:385286/143‑1 (MQ=255)

                                                                                                                                                                                                                                                              cATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACttttt                                                                                                                                                  >  5:544009/1‑131 (MQ=255)

                                                                                                                                                                                                                                                               aTCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTAt                                                                                                                                                                                                                                 >  4:663192/1‑51 (MQ=255)

                                                                                                                                                                                                                                                                 cccGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGa                                                                                                                              >  5:76857/1‑148 (MQ=255)

                                                                                                                                                                                                                                                                    gcgcAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGa                                                                                                                                                                                                    >  2:626467/1‑75 (MQ=255)

                                                                                                                                                                                                                                                                    gcgcAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGAt                                                                                                                    >  2:787312/1‑155 (MQ=255)

                                                                                                                                                                                                                                                                     cgcAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTAt                                                                                                                                                                                                                                 >  2:452560/1‑45 (MQ=255)

                                                                                                                                                                                                                                                                        aGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGAtgtg                                                                                                                                                                                                <  4:626465/75‑1 (MQ=255)

                                                                                                                                                                                                                                                                        aGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTc                                                                                                                <  4:787310/155‑1 (MQ=255)

TGCATCATGGCACACCAGAACAGGCACAGATGATCCGTACCGCCATCGAACAGGGTAACGGTCGCCATCTTCTGGAACCGGTTCTGGAAGCAATGAACGCTTGTGGATCTCTTGAATGGACGCGTCAGCGTGCCGAGGAAGAAGCAGACAAAGCCATCGCAGCGTTACAGGTG‑CTCCCGGACACCCCTTGGCGAGAAGCACTCATCGGCCTCGCGCACATCGCTGTTCAACGCGATCGTTAATCCCCTCCCCTCATCCCGCGCAGC‑GCGCGGGATGGTTCCAGTAAGTTCCATAAAAACACTTATTCAGCTCTATAACACCTGAAAAATAGGGATGTGCATCTTTATATATTCTGAATATTCACACTCTTTACAGGAACTTTTTAGAGCAATAGGCCATCAGGAGTATAGTGATGCTCGACAGAAGAAGTGTTCTGAATGAAAGCGACTAACTTAAGGAGTGAGGAAAAATGGAAAGTAATTTCATTGACTGGCATCCCGCTGACATCATTGCGGGTTTGCGCAAGAA  >  NC_000913/3334440‑3334967

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