BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	annotation	gene	description
RA	NC_000913	3,948,805	(G)_6→7	coding (1165/1521 nt)	yifB ←	putative magnesium chelatase YifB

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	3,948,799	1	.	G	98.8%	321.3 / ‑3.3	85	coding (1171/1521 nt)	yifB	putative magnesium chelatase YifB
Reads supporting (aligned to +/- strand): ref base . (1/0); new base G (34/50); total (35/50)
Fisher's exact test for biased strand distribution p-value = 4.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.44e-01

Only 100 of 111 total aligned reads displayed.

CAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCC‑‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGC  >  NC_000913/3948660‑3948930

cAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑Gcc                                                                                                                                        <  2:151059/140‑1 (MQ=255)

 aGAATTGGTGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCC‑‑g                                                                                                                                     <  2:790647/140‑1 (MQ=255)

  gAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCC‑‑gg                                                                                                                                    <  2:539965/140‑1 (MQ=255)

      tGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCC‑‑gggggg                                                                                                                                <  2:602621/140‑1 (MQ=255)

      tGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCC‑‑gggggg                                                                                                                                <  1:818090/140‑1 (MQ=255)

       ggCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACTGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCC‑‑ggggggg                                                                                                                               <  1:1280331/140‑2 (MQ=255)

       ggCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCC‑‑ggggggg                                                                                                                               <  1:269675/140‑2 (MQ=255)

       ggCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCC‑‑ggggggg                                                                                                                               <  1:331310/140‑2 (MQ=255)

       ggCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCC‑‑ggggggg                                                                                                                               <  2:235135/140‑2 (MQ=255)

       ggCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCC‑‑ggggggg                                                                                                                               <  2:1491334/140‑2 (MQ=255)

       ggCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCC‑‑ggggggg                                                                                                                               <  1:782698/140‑2 (MQ=255)

       ggCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCC‑‑ggggggg                                                                                                                               <  1:1274273/140‑2 (MQ=255)

       ggCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCC‑‑ggggggg                                                                                                                               <  1:914656/140‑2 (MQ=255)

       ggCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCC‑‑ggggggg                                                                                                                               <  2:744504/140‑2 (MQ=255)

       ggCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCC‑‑ggggggg                                                                                                                               <  2:1170150/140‑2 (MQ=255)

        gCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCC‑‑gggggggt                                                                                                                              <  2:1022832/140‑3 (MQ=255)

        gCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCC‑‑gggggggg                                                                                                                              >  1:935436/1‑138 (MQ=255)

         cGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCC‑‑gggggggtg                                                                                                                             <  1:405714/140‑4 (MQ=255)

                 cGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATggg                                                                                                                     <  1:579303/140‑1 (MQ=255)

                  ggACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGa                                                                                                                    <  2:1267685/140‑1 (MQ=255)

                  ggACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGa                                                                                                                    <  1:220963/140‑1 (MQ=255)

                  ggACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGa                                                                                                                    <  1:664061/140‑1 (MQ=255)

                  ggACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGa                                                                                                                    <  2:167730/140‑1 (MQ=255)

                   gACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGAt                                                                                                                   <  1:742067/140‑1 (MQ=255)

                      tATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCGCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGAtctc                                                                                                                <  2:545679/140‑1 (MQ=255)

                        tCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCa                                                                                                              >  2:1474616/1‑140 (MQ=255)

                        tCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCa                                                                                                              >  2:945705/1‑140 (MQ=255)

                        tCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCa                                                                                                              >  2:1208082/1‑140 (MQ=255)

                        tCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCa                                                                                                              >  2:1258468/1‑140 (MQ=255)

                        tCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCa                                                                                                              >  1:335535/1‑140 (MQ=255)

                        tCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCa                                                                                                              >  1:1521057/1‑140 (MQ=255)

                        tCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑GGGTAATGGGATCTCCa                                                                                                              >  1:406198/1‑140 (MQ=255)

                        tCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGATGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCa                                                                                                              >  1:799934/1‑140 (MQ=255)

                                  cATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATc                                                                                                    >  1:185020/1‑140 (MQ=255)

                                    ttCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAATGGCG‑‑GGGGGG‑AGGTAATCGGATCTCCAGTGAGAGGTCGa                                                                                                  >  1:186611/1‑140 (MQ=255)

                                      cAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCGGGGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGaa                                                                                                 >  2:863896/1‑140 (MQ=255)

                                       aGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGt                                                                                                            >  1:1002014/1‑127 (MQ=255)

                                       aGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGt                                                                                                            <  2:1002014/127‑1 (MQ=255)

                                       aGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGc                                                                                               >  1:759222/1‑140 (MQ=255)

                                                tGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGaa                                                                                      <  1:808191/140‑1 (MQ=255)

                                                 gCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGaaa                                                                                     >  1:1305082/1‑140 (MQ=255)

                                                 gCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGaaa                                                                                     >  1:986588/1‑140 (MQ=255)

                                                 gCCGCTTAAATTGGCGCTCTCTGGCGGCCATTAAGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAACT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGaaa                                                                                     >  2:432986/1‑140 (MQ=255)

                                                  ccGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAg                                                                                    >  1:1445939/1‑140 (MQ=255)

                                                   cGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAgg                                                                                   >  2:445667/1‑140 (MQ=255)

                                                   cGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAgg                                                                                   <  2:1474226/140‑1 (MQ=255)

                                                   cGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAgg                                                                                   <  2:1511730/140‑1 (MQ=255)

                                                     cTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGc                                                                                 >  2:681272/1‑140 (MQ=255)

                                                         aaTTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCg                                                                             <  1:587864/140‑1 (MQ=255)

                                                             ggCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGagag                                                                         <  1:659309/140‑1 (MQ=255)

                                                              gcgcTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCGGGGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGagag                                                                         <  1:863896/140‑1 (MQ=255)

                                                              gcgcTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTGTTACTCAAAAT‑GCC‑‑GGGGGGTTGGTAATGGGATCTCCAGTGAGAGGTCGAAGCGGTCGAGAAAGGGCCCCGAGAGc                                                                        <  2:668800/140‑1 (MQ=255)

                                                               cgcTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGcc                                                                       >  2:736386/1‑140 (MQ=255)

                                                                  tctctGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGt                                                                    <  1:1033479/140‑1 (MQ=255)

                                                                    tctGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTg                                                                  <  2:237891/140‑1 (MQ=255)

                                                                      tGGCGGTCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTgag                                                                >  2:149640/1‑140 (MQ=255)

                                                                      tGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTgag                                                                <  2:72770/140‑1 (MQ=255)

                                                                      tGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTgag                                                                <  2:1331314/140‑1 (MQ=255)

                                                                      tGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTgag                                                                <  2:437602/140‑1 (MQ=255)

                                                                       ggcggtCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTgaga                                                               <  1:954526/140‑1 (MQ=255)

                                                                       ggcggcCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTgaga                                                               <  2:92204/140‑1 (MQ=255)

                                                                       ggcggcCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTgaga                                                               <  2:243192/140‑1 (MQ=255)

                                                                        gcggcCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGAt                                                              <  2:935436/140‑1 (MQ=255)

                                                                        gcggcCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGAt                                                              <  1:592023/140‑1 (MQ=255)

                                                                        gcggcCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGAt                                                              <  1:488215/140‑1 (MQ=255)

                                                                            ccATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACg                                                          <  1:676709/140‑1 (MQ=255)

                                                                             cATTACACGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGt                                                         <  2:1137045/140‑1 (MQ=255)

                                                                              aTTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTa                                                        >  2:1036347/1‑140 (MQ=255)

                                                                                tACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAAt                                                      <  1:1474616/140‑1 (MQ=255)

                                                                                  cgcgTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGt                                                    <  1:1038633/140‑1 (MQ=255)

                                                                                  cgcgTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGt                                                    <  2:1521057/140‑1 (MQ=255)

                                                                                  cgcgATGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGt                                                    <  2:799934/140‑1 (MQ=255)

                                                                                      ttgttTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGt                                                >  1:213255/1‑140 (MQ=255)

                                                                                                 tGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCa                                     <  1:945705/140‑1 (MQ=255)

                                                                                                 tGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCa                                     <  1:1258468/140‑1 (MQ=255)

                                                                                                 tGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCa                                     <  1:33107/140‑1 (MQ=255)

                                                                                                 tGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑AGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCa                                     >  1:826620/1‑140 (MQ=255)

                                                                                                  ggCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTgaag                                    <  1:1114926/140‑4 (MQ=255)

                                                                                                   gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGc                                   >  2:906843/1‑140 (MQ=255)

                                                                                                   gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGc                                   >  1:469820/1‑140 (MQ=255)

                                                                                                   gcgcTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGc                                   >  2:218915/1‑140 (MQ=255)

                                                                                                        gctTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTa                              <  1:1208082/140‑1 (MQ=255)

                                                                                                         ctTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTAt                             >  2:518605/1‑140 (MQ=255)

                                                                                                         ctTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTAt                             <  2:131255/140‑1 (MQ=255)

                                                                                                          tttCTCCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATg                            <  2:607015/140‑1 (MQ=255)

                                                                                                              tcCCGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGttt                        <  1:681272/140‑1 (MQ=255)

                                                                                                               cccGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTc                       >  2:1480266/1‑140 (MQ=255)

                                                                                                                 cGGCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTccc                     <  2:366190/140‑1 (MQ=255)

                                                                                                                   gCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTg                   >  1:68425/1‑140 (MQ=255)

                                                                                                                   gCACTACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTg                   >  2:907663/1‑140 (MQ=255)

                                                                                                                       tACCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGGCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATaa               <  1:60502/140‑1 (MQ=255)

                                                                                                                        aCCGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAAt              >  2:546141/1‑140 (MQ=255)

                                                                                                                         ccGTTTTACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATg             <  1:144641/140‑1 (MQ=255)

                                                                                                                              ttACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCgg        >  1:1507709/1‑140 (MQ=255)

                                                                                                                               tACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGt       <  2:135071/140‑1 (MQ=255)

                                                                                                                               tACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGt       <  1:810203/140‑1 (MQ=255)

                                                                                                                               tACTCAAAAT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGt       <  2:12660/140‑1 (MQ=255)

                                                                                                                                 cTCAAAAT‑GCC‑‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAgg    >  1:672669/1‑140 (MQ=255)

                                                                                                                                    aaaaT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGc  <  1:1050994/140‑1 (MQ=255)

                                                                                                                                    aaaaT‑GCCG‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGc  <  1:937524/140‑1 (MQ=255)

CAGAATTGGCGTATTTCCGGACTATCCAGCCAGGCATTCAGTTTATTCTGCCGCTTAAATTGGCGCTCTCTGGCGGCCATTACGCGTTGTTTAACGGTGGCGCTGCTTTCTCCCGGCACTACCGTTTTACTCAAAAT‑GCC‑‑GGGGGG‑TGGTAATGGGATCTCCAGTGAGAGATCGAAGCGGTCGAGAAAGGGCCCCGAGAGCCGGTTGAGATAACGTAATGTCTGTTCTGGCGTGCAGCGGTTATGGTTTCCCTGATAATGTCCGGTAGGGC  >  NC_000913/3948660‑3948930

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