BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	1,585,156:1	+C	5.4%	coding (1331/2280 nt)	ydeP ←	putative oxidoreductase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	1,585,156	1	.	C	5.4%	492.7 / 21.7	112	coding (1331/2280 nt)	ydeP	putative oxidoreductase
Reads supporting (aligned to +/- strand): ref base . (52/54); new base C (2/4); total (54/58)
Fisher's exact test for biased strand distribution p-value = 6.80e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 114 total aligned reads displayed.

GCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCAG  >  NC_000913/1585010‑1585303

gCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg‑ca                                                                                                                                                   >  4:286911/1‑149 (MQ=255)

gCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTg‑ca                                                                                                                                                   <  2:849670/149‑1 (MQ=255)

    aGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACTGTACAGCGCTCGCTTCCCGATCTGGCATTGCCACCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTGGCATGCT‑GCAATTGCCAGCAt                                                                                                                                               <  2:162291/149‑1 (MQ=255)

    aGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCAt                                                                                                                                               <  1:138599/149‑1 (MQ=255)

      tGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCT‑GCAATTGCCAGCATgc                                                                                                                                             <  2:844415/149‑2 (MQ=255)

        gCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGt                                                                                                                                             <  3:164364/147‑1 (MQ=255)

        gCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTcc                                                                                                                                           <  3:58238/149‑1 (MQ=255)

         cTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCa                                                                                                                                          >  4:214807/1‑149 (MQ=255)

          tACTTGTACCTCCAAATGTAATTGCGTTAACGGTACAGCGCTTGTTTCCCGATCTGGCATTGCCAGATCAAAATTGCCCCCCATACAGATCAATGCTGTAGCCTGCCCTGTACATTTTGCTGGAATGCTGGCAATTG‑CAGCATGTCCAg                                                                                                                                         <  4:447498/149‑1 (MQ=255)

          tACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAg                                                                                                                                         >  3:540793/1‑149 (MQ=255)

          tACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAg                                                                                                                                         >  2:790782/1‑149 (MQ=255)

          tACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAg                                                                                                                                         >  4:511731/1‑149 (MQ=255)

          tACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAg                                                                                                                                         <  1:353943/149‑1 (MQ=255)

                tACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCAt                                                                                                                                   <  3:452318/149‑1 (MQ=255)

                  ccGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGa                                                                                                                                 <  2:540091/149‑1 (MQ=255)

                    gCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAgg                                                                                                                               >  3:82627/1‑149 (MQ=255)

                    gCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAgg                                                                                                                               >  1:890464/1‑149 (MQ=255)

                     ccAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGt                                                                                                                                       >  1:471071/1‑140 (MQ=255)

                     ccAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAgg                                                                                                                               >  1:370370/1‑148 (MQ=255)

                     ccAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAgg                                                                                                                               <  2:370370/148‑1 (MQ=255)

                     ccAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAgg                                                                                                                               >  3:613756/1‑148 (MQ=255)

                     ccAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAgg                                                                                                                               <  4:613756/148‑1 (MQ=255)

                     ccAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGt                                                                                                                              >  1:841770/1‑149 (MQ=255)

                     ccAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGt                                                                                                                              >  4:754589/1‑149 (MQ=255)

                     ccAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGt                                                                                                                              >  3:180124/1‑149 (MQ=255)

                       aaaTCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTgg                                                                                                                            <  3:286911/149‑1 (MQ=255)

                         aTCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTgggg                                                                                                                          <  4:340398/149‑1 (MQ=255)

                          tCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGt                                                                                                                         <  1:699561/149‑1 (MQ=255)

                             aaTTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGaa                                                                                                                      <  2:387327/149‑1 (MQ=255)

                              aTTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAg                                                                                                                     >  4:806305/1‑149 (MQ=255)

                              aTTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAg                                                                                                                     >  2:670744/1‑149 (MQ=255)

                               ttGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGc                                                                                                                    >  4:552413/1‑149 (MQ=255)

                               ttGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGc                                                                                                                    >  4:213112/1‑149 (MQ=255)

                                   gTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCGAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATa                                                                                                                <  3:754589/149‑1 (MQ=255)

                                   gTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATa                                                                                                                <  4:782505/149‑1 (MQ=255)

                                    ttAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAg                                                                                                               >  4:486470/1‑149 (MQ=255)

                                       aCGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAgcgc                                                                                                            <  4:818402/149‑1 (MQ=255)

                                           tACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCAc                                                                                                        >  4:716330/1‑149 (MQ=255)

                                              aGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCa                                                                                                     >  4:678184/1‑149 (MQ=255)

                                                       ttCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTAGATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCa                                                                                            >  4:62163/1‑149 (MQ=255)

                                                       ttCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGc                                                                                                                    >  3:561911/1‑125 (MQ=255)

                                                       ttCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCa                                                                                            <  2:513901/149‑1 (MQ=255)

                                                       ttCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCa                                                                                            >  2:849610/1‑149 (MQ=255)

                                                       ttCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCa                                                                                            >  2:892210/1‑149 (MQ=255)

                                                         cccGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATCG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGa                                                                                          >  4:607939/1‑149 (MQ=255)

                                                           tgaccTGGCATTGCTAGCGCAAAATTGCCCCTCATGCAGATCAATGCTCGAGCCTGCCCCGTATATATGGCTTGCATGCTGGCAATAG‑CAGCGTGTCCAGGTGCACGAGGTGCGGTGAAGCCATAGAGCTCACCCAGACGAGGCAGaaa                                                                                        <  3:545797/145‑1 (MQ=255)

                                                           cGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGaaa                                                                                        <  4:146690/149‑1 (MQ=255)

                                                           cGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGaaa                                                                                        <  1:546625/149‑1 (MQ=255)

                                                           cGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGaaa                                                                                        <  3:618404/149‑1 (MQ=255)

                                                             aTCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAAct                                                                                      >  2:456277/1‑149 (MQ=255)

                                                               cTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAActct                                                                                    <  2:583310/149‑1 (MQ=255)

                                                               cTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAActct                                                                                    <  1:488168/149‑1 (MQ=255)

                                                                 ggCATTGCCAGCGCAAAATTGCCCCCCTTGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGc                                                                                  <  1:679908/149‑1 (MQ=255)

                                                                  gCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCa                                                                                 <  1:5411/149‑1 (MQ=255)

                                                                    aTTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGa                                                                               >  3:536345/1‑149 (MQ=255)

                                                                            cgcAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTtct                                                                       <  4:513943/149‑1 (MQ=255)

                                                                              cAAAACTGCCACCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGATTCTCg                                                                     <  4:498566/149‑1 (MQ=255)

                                                                               aaaaTCGACCCCAAAGCAGAACAATTCGCGAGTGTGCCCCGGACATATTGGTTGCATCATGGCAATTG‑CAGCATGTCCAGGTCCGTGAGGTGGTGCCAAGCCTGAGGGTTAATCCAGATGAGCCGGAAACTCAGCAGACGGTTTCTCgg                                                                    <  2:445907/149‑1 (MQ=255)

                                                                                aaaTTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGt                                                                   <  1:558426/149‑1 (MQ=255)

                                                                                aaaTTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGt                                                                   <  1:131653/149‑1 (MQ=255)

                                                                                      ccccccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATAcc                                                             <  2:841770/149‑1 (MQ=255)

                                                                                         cccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGAc                                                          <  4:759518/149‑1 (MQ=255)

                                                                                         cccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGAc                                                          >  1:862081/1‑149 (MQ=255)

                                                                                          ccATGCAGATCTATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACg                                                         >  3:514412/1‑149 (MQ=255)

                                                                                          ccATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACg                                                         >  1:116290/1‑149 (MQ=255)

                                                                                           cATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACgg                                                        <  1:629986/149‑1 (MQ=255)

                                                                                                 gATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTCCGGt                                                  >  4:869895/1‑149 (MQ=255)

                                                                                                     aaTGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCAcc                                                                                                       <  4:64096/92‑1 (MQ=255)

                                                                                                     aaTGCTCGAGCCTGACCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCAcc                                                                                                       >  3:64096/1‑92 (MQ=255)

                                                                                                      aTGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGccc                                             >  4:269884/1‑149 (MQ=255)

                                                                                                         cTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGt                                          >  2:376/1‑149 (MQ=255)

                                                                                                            gAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACa                                       <  2:115861/149‑1 (MQ=255)

                                                                                                             aGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTAACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACAt                                      <  2:814846/149‑1 (MQ=255)

                                                                                                              gCCTGCCCCGTACATATTGGTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACAtt                                     >  3:696459/1‑149 (MQ=255)

                                                                                                                 tGCCCCGTACATATTGCTTGGATGCGGGCAATTG‑CAGCATGTCCAGGTGCATGAGGGGGGGTGAAGCCATAGCGCTCACCCATACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTaga                                  <  4:173247/149‑1 (MQ=255)

                                                                                                                 tGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATTTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTaga                                  <  3:759870/149‑1 (MQ=255)

                                                                                                                   ccccGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAgt                                <  3:213112/149‑1 (MQ=255)

                                                                                                                    cccGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAgtg                               >  4:579932/1‑149 (MQ=255)

                                                                                                                     ccGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTg                                           <  2:818958/136‑1 (MQ=255)

                                                                                                                       gTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTcc                            <  2:516956/149‑1 (MQ=255)

                                                                                                                          cATATTGCTTGCATGCTGGCAATTG‑CAGTATGTCCAGGTGTATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGCAACACTGCCGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCCCg                         <  4:238140/149‑1 (MQ=255)

                                                                                                                           atatTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGt                        >  4:156958/1‑149 (MQ=255)

                                                                                                                              ttgcttgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGt                     >  2:469229/1‑149 (MQ=255)

                                                                                                                                  ttgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGa                                                                 >  2:485912/1‑101 (MQ=255)

                                                                                                                                  ttgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGa                                                                 <  1:485912/101‑1 (MQ=255)

                                                                                                                                  ttgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGc                 >  2:143487/1‑149 (MQ=255)

                                                                                                                                  ggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTaga                                  >  3:443452/4‑132 (MQ=255)

                                                                                                                                  ggacATGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTaga                                  <  4:443452/129‑1 (MQ=255)

                                                                                                                                   tgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTaga                                  >  4:636075/1‑131 (MQ=255)

                                                                                                                                   tgcATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTaga                                  <  3:636075/131‑1 (MQ=255)

                                                                                                                                     cATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGa              <  3:678184/149‑1 (MQ=255)

                                                                                                                                      aTGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGAt             <  2:886069/149‑1 (MQ=255)

                                                                                                                                      aTGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGAt             >  2:188243/1‑149 (MQ=255)

                                                                                                                                      aTGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGt                                                                   <  3:107012/95‑1 (MQ=255)

                                                                                                                                      aTGCT‑GCAATTGCCAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGt                                                                   >  4:107012/1‑95 (MQ=255)

                                                                                                                                       tGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATa            >  2:660243/1‑149 (MQ=255)

                                                                                                                                           ggCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCg        >  3:785616/1‑149 (MQ=255)

                                                                                                                                            gCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGc       >  2:233685/1‑149 (MQ=255)

                                                                                                                                                ttG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCa   >  1:555241/1‑149 (MQ=255)

                                                                                                                                                 tG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCAg  <  3:607939/149‑1 (MQ=255)

GCTTAGTGGCTACGTGTACCGCCAAATCTAATTGCGTTAACGGTACAGCGCTCGCTTCCCGATCTGGCATTGCCAGCGCAAAATTGCCCCCCATGCAGATCAATGCTCGAGCCTGCCCCGTACATATTGCTTGCATGCTGGCAATTG‑CAGCATGTCCAGGTGCATGAGGTGGGGTGAAGCCATAGCGCTCACCCAGACGAGCCAGAAACTCTGCAGACGGTTTCTCGGTGATACCGACGGTTCGGTCGCCCTGTACATTAGAGTGTCCACGTAGTGGGCAGATACCCGCACCAG  >  NC_000913/1585010‑1585303

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