BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,447,964	T→A	6.1%	L24* (TTA→TAA)	ytfP →	GGCT‑like protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,447,964	0	T	A	6.1%	407.7 / 13.4	115	L24* (TTA→TAA)	ytfP	GGCT‑like protein
Reads supporting (aligned to +/- strand): ref base T (59/49); new base A (3/4); total (62/53)
Fisher's exact test for biased strand distribution p-value = 7.02e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.13e-01

Only 100 of 115 total aligned reads displayed.

TGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGG  >  NC_000913/4447818‑4448109

tGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTAc                                                                                                                                                 <  3:16838/149‑1 (MQ=255)

 gCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACt                                                                                                                                                <  1:493068/149‑1 (MQ=255)

   cTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTgg                                                                                                                                              <  3:499173/149‑1 (MQ=255)

    tAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTggg                                                                                                                                             >  1:26163/1‑149 (MQ=255)

       gCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGa                                                                                                                                          >  4:621056/1‑149 (MQ=255)

        cTATATCTGGAAGCCGTGTGTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACGGTCACTGGATGACCAATGCCCAGTTACTGGGCGAt                                                                                                                                         <  4:470949/149‑1 (MQ=255)

         tataTCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGTCCAGTAACTGGGcattg                                                                                                                                        <  2:119927/149‑5 (MQ=255)

         tataTCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTAACTGGGcattg                                                                                                                                        <  4:223762/149‑5 (MQ=255)

           tatCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTc                                                                                                                                      >  2:52520/1‑149 (MQ=255)

             tCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAg                                                                                                                                    >  2:364390/1‑149 (MQ=255)

                ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt                                                                                                                                 <  3:49643/149‑1 (MQ=255)

                ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt                                                                                                                                 >  3:201310/1‑149 (MQ=255)

                ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt                                                                                                                                 <  3:80123/149‑1 (MQ=255)

                ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTAACTGGGcattg                                                                                                                                        <  3:275690/142‑5 (MQ=255)

                ggAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTAACTGGGcattg                                                                                                                                        >  4:275690/1‑138 (MQ=255)

                 gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATc                                                                                                                                >  2:153900/1‑149 (MQ=255)

                 gAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATc                                                                                                                                >  4:630276/1‑149 (MQ=255)

                  aaGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCg                                                                                                                               >  3:15669/1‑149 (MQ=255)

                      cGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATaa                                                                                                                           <  1:160316/149‑1 (MQ=255)

                      cGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATaa                                                                                                                           <  2:326279/149‑1 (MQ=255)

                        tgtCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACt                                                                                                                         <  1:153900/149‑1 (MQ=255)

                         gtCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTa                                                                                                                        >  4:461304/1‑149 (MQ=255)

                         gccTGGTCTACACCAGACACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTa                                                                                                                        >  3:619588/3‑149 (MQ=255)

                                 tAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGt                                                                                                                <  4:612810/149‑1 (MQ=255)

                                    aCCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGtata                                                                                                             <  1:364390/149‑1 (MQ=255)

                                     ccAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAg                                                                                                            >  4:430223/1‑149 (MQ=255)

                                     ccAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAg                                                                                                            >  4:510668/1‑149 (MQ=255)

                                      cAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGc                                                                                                           <  3:412376/149‑1 (MQ=255)

                                      cAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGc                                                                                                           <  3:104129/149‑1 (MQ=255)

                                       aGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGcc                                                                                                          >  1:254746/1‑149 (MQ=255)

                                         gCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTg                                                                                                        <  2:585859/149‑1 (MQ=255)

                                         acACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTg                                                                                                        <  4:619588/148‑1 (MQ=255)

                                           aCTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTggg                                                                                                      >  1:595386/1‑149 (MQ=255)

                                            cTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTa                                                                                                                        <  3:20983/130‑1 (MQ=255)

                                            cTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTa                                                                                                                        >  4:20983/1‑130 (MQ=255)

                                             tGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGcc                                                                                                    >  3:194191/1‑149 (MQ=255)

                                              ggATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCa                                                                                                   <  3:409157/149‑1 (MQ=255)

                                                    gCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATcc                                                                                             <  3:645564/149‑1 (MQ=255)

                                                      tctATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAg                                                                                           >  1:356956/1‑149 (MQ=255)

                                                        tATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCTCAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt                                                                                                                                 <  2:565392/109‑1 (MQ=255)

                                                        tATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTAt                                                                                                                                 >  1:565392/1‑109 (MQ=255)

                                                        tATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGgc                                                                                         >  1:634809/1‑149 (MQ=255)

                                                        tATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGgc                                                                                         <  4:193176/149‑1 (MQ=255)

                                                         aTCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGgcg                                                                                        >  3:612046/1‑149 (MQ=255)

                                                          tCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGgcgc                                                                                       >  2:535810/1‑149 (MQ=255)

                                                          tCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGgcgc                                                                                       >  3:661015/1‑149 (MQ=255)

                                                             gTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGt                                                                                    >  3:126190/1‑149 (MQ=255)

                                                              ttCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGtt                                                                                   >  3:180744/1‑149 (MQ=255)

                                                                cGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTcc                                                                                 <  3:592279/149‑1 (MQ=255)

                                                                 gAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCg                                                                                <  3:118171/149‑1 (MQ=255)

                                                                 gAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCg                                                                                >  1:658032/1‑149 (MQ=255)

                                                                      ttAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGggaa                                                                           <  3:510668/149‑1 (MQ=255)

                                                                         gCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGggaacgg                                                                        >  2:419870/1‑149 (MQ=255)

                                                                           aaTGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGggaacggaa                                                                      >  3:153243/1‑149 (MQ=255)

                                                                           aaTGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGggaacggaa                                                                      <  3:430223/149‑1 (MQ=255)

                                                                                 aaTATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTac                                                                >  1:26316/1‑149 (MQ=255)

                                                                                        gTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGa                                                         >  4:239591/1‑149 (MQ=255)

                                                                                         tCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGaa                                                        >  4:266163/1‑149 (MQ=255)

                                                                                         tCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGaa                                                        <  2:356956/149‑1 (MQ=255)

                                                                                           tACGGCAGTTTACGCCACAAACAAGGCAACAGTCACAGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCGGGGGAACGGAACGGTACACGGTGAAGt                                                      <  3:461304/149‑1 (MQ=255)

                                                                                            aCGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGtt                                                     <  1:419870/149‑1 (MQ=255)

                                                                                             cGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTgg                                                                                                       <  2:672825/98‑1 (MQ=255)

                                                                                             cGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTgg                                                                                                       >  1:672825/1‑98 (MQ=255)

                                                                                                cAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATc                                                 <  2:634809/149‑1 (MQ=255)

                                                                                                cAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGAAACGGAACGGTACACGGTGAAGTTTATc                                                 <  4:612046/149‑1 (MQ=255)

                                                                                                   tttACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTa                                              <  3:621056/149‑1 (MQ=255)

                                                                                                   ttaaCGCCACAAACAAGGAAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATACAGGCGCAGCTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTa                                              >  4:614445/4‑149 (MQ=255)

                                                                                                     tACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTAtt                                            <  2:139187/149‑1 (MQ=255)

                                                                                                          cacaAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGtata                                                                                                             <  2:91702/79‑1 (MQ=255)

                                                                                                          cacaAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGtata                                                                                                             >  1:91702/1‑79 (MQ=255)

                                                                                                          cacaAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACaa                                       <  4:201310/149‑1 (MQ=255)

                                                                                                           acaaacaaGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATcc                                                                                             <  2:569698/94‑1 (MQ=255)

                                                                                                           acaaacaaGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATcc                                                                                             >  1:569698/1‑94 (MQ=255)

                                                                                                           acaaacaaGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACaa                                       <  4:349645/148‑1 (MQ=255)

                                                                                                           acaaacaaGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACaa                                       >  3:349645/1‑148 (MQ=255)

                                                                                                           acaaacaaGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAAc                                      >  1:570155/1‑149 (MQ=255)

                                                                                                           acaaacaaGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAAc                                      >  2:678902/1‑149 (MQ=255)

                                                                                                           acaaacaaGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAAc                                      >  3:269903/1‑149 (MQ=255)

                                                                                                               acaaGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTTCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCa                                  >  3:570708/1‑149 (MQ=255)

                                                                                                                   ggCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCt                              >  3:471730/1‑149 (MQ=255)

                                                                                                                      aaCAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGc                           >  3:568231/1‑149 (MQ=255)

                                                                                                                           tCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAAc                      <  4:471730/149‑1 (MQ=255)

                                                                                                                               tGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGa                  >  2:648072/1‑149 (MQ=255)

                                                                                                                                ggATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGAt                 <  4:126190/149‑1 (MQ=255)

                                                                                                                                  aTGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGc               >  2:203138/1‑149 (MQ=255)

                                                                                                                                    gACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCt             <  3:239591/149‑1 (MQ=255)

                                                                                                                                      ccAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTAt                                             >  1:498410/1‑115 (MQ=255)

                                                                                                                                      ccAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTAt                                             <  2:498410/115‑1 (MQ=255)

                                                                                                                                         tcgCCCAGTAACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACgg                                                            <  1:333321/95‑1 (MQ=255)

                                                                                                                                         tcgCCCAGTAACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACgg                                                            >  2:333321/3‑97 (MQ=255)

                                                                                                                                         tcgCCCAGTAACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTgcgc        >  1:317461/3‑149 (MQ=255)

                                                                                                                                           gCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCAc      >  4:591465/1‑149 (MQ=255)

                                                                                                                                             ccAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGAt                 <  1:632901/136‑1 (MQ=255)

                                                                                                                                             ccAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGAt                 >  2:632901/1‑136 (MQ=255)

                                                                                                                                             ccAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCa    >  1:109780/1‑149 (MQ=255)

                                                                                                                                             ccAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCa    >  3:359642/1‑149 (MQ=255)

                                                                                                                                               aGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCgg                                                                               >  2:281684/1‑72 (MQ=255)

                                                                                                                                               aGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCgg                                                                               <  1:281684/72‑1 (MQ=255)

                                                                                                                                               aGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAgg  <  2:658032/149‑1 (MQ=255)

                                                                                                                                               aGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAgg  <  2:26163/149‑1 (MQ=255)

TGCCTAAGCTATATCTGGAAGCCGTGTCTGGTGTAGACCAGGCACTGGATTTGCTCTATCAGTTCGAGTTTTAGCAATGCGAATATTTGTCTACGGCAGTTTACGCCACAAACAAGGCAACAGTCACTGGATGACCAATGCCCAGTTACTGGGCGATTTCAGTATCGATAACTACCAGTTGTATAGCCTGGGCCACTATCCAGGCGCAGTTCCGGGGAACGGAACGGTACACGGTGAAGTTTATCGTATTGACAACGCCACGCTGGCCGAACTTGATGCCTTGCGCACCAGG  >  NC_000913/4447818‑4448109

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