BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	3,689,667	T→A	5.7%	T199S (ACC→TCC)	bcsZ ←	endo‑1,4‑D‑glucanase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,689,667	0	T	A	5.7%	276.3 / 4.2	101	T199S (ACC→TCC)	bcsZ	endo‑1,4‑D‑glucanase
Reads supporting (aligned to +/- strand): ref base T (51/42); new base A (2/4); total (53/48)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.38e-01

Only 100 of 102 total aligned reads displayed.

GTAAACGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCC  >  NC_000913/3689518‑3689812

gTAAACGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGt                                                                                                                                                     <  2:1661314/150‑1 (MQ=255)

  aaaCGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTcc                                                                                                                                                   >  2:634191/1‑150 (MQ=255)

    aCGCGGATAGCGTCGTTGCTGCTGTTCAATGTTTTTTCTGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAACAACCTTTCGGGGCGGGTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCAc                                                                                                                                                 <  1:81147/150‑1 (MQ=255)

      gcgGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACg‑‑g                                                                                                                                             >  2:1360976/1‑150 (MQ=255)

          aTAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑gcgcg                                                                                                                                         <  2:1402731/150‑1 (MQ=255)

            ttcgtcgtAGCTTCTGTTGAATGTTTTTTCTGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCTCACCCATTCTGGCGAAAAACCTTTCGGGTCTGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCGCCGTGGACCACG‑‑atgcgcg                                                                                                                                       <  2:719858/148‑2 (MQ=255)

             gCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCa                                                                                                                                      >  1:1579404/1‑150 (MQ=255)

              cgtcgtAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGC‑‑CGTGGACCACG‑‑CTGCGCGaaac                                                                                                                                   <  1:1438908/150‑2 (MQ=255)

                  gtAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCg                                                                                                                                 <  1:1000089/150‑1 (MQ=255)

                    aGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGGTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCggg                                                                                                                               <  2:1229074/150‑1 (MQ=255)

                      ctgctgATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTg                                                                                                                             >  2:119930/1‑150 (MQ=255)

                      ctgcggCTCAATGTTTTTTCGGCTTTTAGCGGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGGCTGGCGAAAAACCTGTCGGGGCGGTTTCCAGCAATACACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTg                                                                                                                             <  2:1435627/150‑1 (MQ=255)

                      atgctgATCAATGTTTATTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCCACAAACCTTTCGGTGCGGTTTCCAGCAATAAACGATGATTGGATTCGCGCAGCGTGGTCAACG‑‑GCGCTCCAAAGCGCGTg                                                                                                                             >  2:125549/2‑150 (MQ=255)

                        gctgATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCGCCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGaa                                                                                                                           >  1:39943/1‑150 (MQ=255)

                           gATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCCCG‑‑GCCCGCCAAAGCGGGTGAAATa                                                                                                                        >  2:1355562/1‑150 (MQ=255)

                            aTCAATGTTTATTCGGCTTTTAGCTGCCAGACTTTGACTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTGGGGGCGGTTTCCAGCAATAAACGATGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATAc                                                                                                                       <  2:310933/150‑1 (MQ=255)

                                aTGTTTATTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCTCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCATAGCGGTTCAAAGACAgcg                                                                                                                   >  1:1075828/1‑150 (MQ=255)

                                    tttttCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACGCAGTCTGTCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCTCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAg                                                                                                               <  1:1242008/150‑1 (MQ=255)

                                      tttCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCg                                                                                                             <  1:1030657/150‑1 (MQ=255)

                                       ttCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGt                                                                                                            >  1:161637/1‑150 (MQ=255)

                                         cGGCTTTTAGCTGCCAGCCTTTGTCTTACTCATAGCGCACCCAGTCTGGCGAAAACCCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTcg                                                                                                          >  2:651152/1‑150 (MQ=255)

                                          ggCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCCCG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTcgg                                                                                                         >  1:1285643/1‑150 (MQ=255)

                                          ggCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTcgg                                                                                                         >  2:1558306/1‑150 (MQ=255)

                                          ggCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTcgg                                                                                                         >  2:1706587/1‑150 (MQ=255)

                                          ggCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTcgg                                                                                                         >  1:432843/1‑150 (MQ=255)

                                          ggCTTTTAGCTGCCAGACTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGACAGAGTcgg                                                                                                         >  2:1293979/1‑150 (MQ=255)

                                           gCTTGTAGCTGCCTGCCTTTGTCTTTCGCATAGCGCACCCAGTCGGGCGAAAAACCTTTCGGGGCGGGCTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTcggc                                                                                                        <  2:1056537/150‑1 (MQ=255)

                                               ttAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCa                                                                                                    >  2:119033/1‑150 (MQ=255)

                                                tAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAg                                                                                                   <  2:1189492/150‑1 (MQ=255)

                                                 aGCTGCCAGCATTTGTCATTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGGTTCCAGCAATAAACGTGGATTGGTTTCGCGCAGCGTGGTCCCCG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAgg                                                                                                  <  1:213185/150‑1 (MQ=255)

                                                   cTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTa                                                                                                >  1:1720640/1‑150 (MQ=255)

                                                        cgCCATTGTCTGTCTCATCGCGCCCCCCGTCTGGCGAACACCCTTTCGGGGCGGTTTCCAGCAATACACGATGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGGGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTgg                                                                                           <  2:1662008/149‑1 (MQ=255)

                                                        aGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTgg                                                                                           >  1:632530/1‑150 (MQ=255)

                                                          ccTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTgggg                                                                                         <  2:256165/150‑1 (MQ=255)

                                                          ccTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTgggg                                                                                         >  2:568092/1‑150 (MQ=255)

                                                          ccTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTgggg                                                                                         <  2:947461/150‑1 (MQ=255)

                                                            tttGTCTTTCTCATAGCGCGCCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGtt                                                                                       <  2:39943/150‑1 (MQ=255)

                                                            tttGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGtt                                                                                       <  2:942098/150‑1 (MQ=255)

                                                            tttGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGACCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGtt                                                                                       >  1:760003/1‑150 (MQ=255)

                                                            tttGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGGTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGtt                                                                                       <  2:904513/150‑1 (MQ=255)

                                                              tGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTaa                                                                                     >  1:472464/1‑150 (MQ=255)

                                                                 cTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAAc                                                                                  >  1:203775/1‑150 (MQ=255)

                                                                  tttCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACg                                                                                 <  1:1663965/150‑1 (MQ=255)

                                                                    tctcATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGACCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGcc                                                                               <  2:760003/150‑1 (MQ=255)

                                                                     ctcATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAACGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCa                                                                              >  1:1590021/1‑150 (MQ=255)

                                                                     ctcATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCa                                                                              <  1:544972/150‑1 (MQ=255)

                                                                     ctcATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGCGGTTAAAACGCCa                                                                              >  1:620180/1‑150 (MQ=255)

                                                                       cATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGCGTTTTAACCCCAGc                                                                            <  1:408176/150‑1 (MQ=255)

                                                                        aTAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCt                                                                           <  1:200593/150‑1 (MQ=255)

                                                                        aTAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCt                                                                           >  1:230674/1‑150 (MQ=255)

                                                                        aTAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTCGTTTCGCGCAGCGTGGTCCACG‑‑CCGCGCCAAAGCGGCTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGcg                                                                           >  2:546690/1‑149 (MQ=255)

                                                                          aGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGt                                                                         >  1:1161633/1‑150 (MQ=255)

                                                                           gcgcACCCAGTCTGGCGAAACACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGtt                                                                        >  1:75065/1‑150 (MQ=255)

                                                                               aCCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATcc                                                                    <  1:1356360/150‑1 (MQ=255)

                                                                               aCCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATcc                                                                    <  1:566469/150‑1 (MQ=255)

                                                                                cccAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCt                                                                   <  2:161637/150‑1 (MQ=255)

                                                                                cccAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGGTATCCt                                                                   >  1:1553262/1‑150 (MQ=255)

                                                                                 ccATTCTGGCGAAAAACCTTTCGGGGCTGTTTCCAGCAATAAACGTAGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGACAGCTGTTATCCTc                                                                  <  2:418687/150‑1 (MQ=255)

                                                                                 ccAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTTGGGTTAAAACGCCAGCTGTTATCCta                                                                  >  2:1266066/1‑149 (MQ=255)

                                                                                      cTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCaa                                                             >  2:1709025/1‑150 (MQ=255)

                                                                                        ggCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCaaac                                                           >  2:856001/1‑149 (MQ=255)

                                                                                         gCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGATATCCTCAGCAAAAc                                                          >  1:652122/1‑150 (MQ=255)

                                                                                          cGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAAcc                                                         >  1:309620/1‑150 (MQ=255)

                                                                                           gAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAAccc                                                        >  2:502409/1‑150 (MQ=255)

                                                                                            aaaaaCCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCa                                                       <  1:1709025/150‑1 (MQ=255)

                                                                                              aaaCCTTTCGGGGCGGCGTCCAGCAATAAACGTTGATTGGGTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGGTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACt                                                     <  1:327612/150‑1 (MQ=255)

                                                                                               aaCCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACtt                                                    >  2:1029314/1‑150 (MQ=255)

                                                                                                   tttCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGcc                                                <  2:849575/150‑1 (MQ=255)

                                                                                                   tttCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGcc                                                >  2:1389210/1‑150 (MQ=255)

                                                                                                      cGGGGCGGTTTCCAGCAATAAACGTCGAGTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCgg                                             <  1:1080362/150‑1 (MQ=255)

                                                                                                       agggCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGt                                            <  1:1135818/149‑1 (MQ=255)

                                                                                                        gggCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTTAAACGCCAGCTTTTATACTCAGCAACACCCCCTATGccagcta                                           >  2:1563602/1‑145 (MQ=255)

                                                                                                          gCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTaac                                         >  1:1500468/1‑150 (MQ=255)

                                                                                                           cGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTaaca                                        <  1:119930/150‑1 (MQ=255)

                                                                                                             gTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTaacaac                                      >  1:632445/1‑150 (MQ=255)

                                                                                                               ttCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACAt                                    <  2:152690/150‑1 (MQ=255)

                                                                                                                tCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATg                                   <  1:634191/150‑1 (MQ=255)

                                                                                                                 ccAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAACACGCCCGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATgg                                  >  2:1221569/1‑150 (MQ=255)

                                                                                                                           aCGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCATCGTCGGCTGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAAcccagccc                        <  2:312086/150‑1 (MQ=255)

                                                                                                                           aCGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAAcccagccc                        >  1:1256841/1‑150 (MQ=255)

                                                                                                                            cGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAAcccagccca                       >  1:1138336/1‑150 (MQ=255)

                                                                                                                              ttGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAgg                     <  1:394559/150‑1 (MQ=255)

                                                                                                                               tGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGc                    <  1:42101/150‑1 (MQ=255)

                                                                                                                                gATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGTTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCCGCAAACCCCACTTTGCCCGGTAACAACATGGAGCCCCGCCCCGGCa                   >  1:123930/1‑150 (MQ=255)

                                                                                                                                gATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGcc                   >  1:624722/1‑149 (MQ=255)

                                                                                                                                 aTTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCAc                  <  2:1218509/150‑1 (MQ=255)

                                                                                                                                 aTTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCAc                  <  2:22402/150‑1 (MQ=255)

                                                                                                                                 aTTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCAc                  <  1:659498/150‑1 (MQ=255)

                                                                                                                                   tGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCg                <  1:946304/150‑1 (MQ=255)

                                                                                                                                    ggTTTCGCGCAGCGTGGGCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGt               <  2:1721527/150‑1 (MQ=255)

                                                                                                                                     gTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCCCTTTGCCCGGTAACAACATGGCACCCAGCCCAGGCACCGTc              >  1:192502/1‑150 (MQ=255)

                                                                                                                                     gTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGGTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTc              >  1:1686962/1‑150 (MQ=255)

                                                                                                                                     ctttGGCGCGCCGTGGACCACGCTGCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCg                >  1:219067/2‑150 (MQ=255)

                                                                                                                                        tCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTcacc           <  2:1248954/150‑1 (MQ=255)

                                                                                                                                        tCGCGCAGCGTGGACCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGAAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTcacc           <  1:960369/150‑1 (MQ=255)

                                                                                                                                          gcgcAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACAGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTTAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTcccccc         >  1:461749/1‑148 (MQ=255)

                                                                                                                                               gCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTtcct    <  2:1579404/150‑1 (MQ=255)

                                                                                                                                                cGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCCGCCCAGGCACCGTCACCACTtcctc   >  2:969678/1‑150 (MQ=255)

                                                                                                                                                cGTGGGCCACG‑‑GAGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCCGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTtcctc   <  2:833437/150‑1 (MQ=255)

                                                                                                                                                 gTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCCGCCCAGGCACCGGCACCCCTtcctcc  >  2:1280045/1‑150 (MQ=255)

GTAAACGCGGATAGCGTCGTAGCTGCTGATCAATGTTTTTTCGGCTTTTAGCTGCCAGCCTTTGTCTTTCTCATAGCGCACCCAGTCTGGCGAAAAACCTTTCGGGGCGGTTTCCAGCAATAAACGTTGATTGGTTTCGCGCAGCGTGGTCCACG‑‑GCGCGCCAAAGCGGGTGAAATACTGCGCCAGCGTCGGCGGCAGGTAGCTGGGGTTAAAACGCCAGCTGTTATCCTCAGCAAAACCCACTTTGCCCGGTAACAACATGGAACCCAGCCCAGGCACCGTCACCACTTCCTCC  >  NC_000913/3689518‑3689812

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