BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,217,943	T→C	5.7%	Y279H (TAT→CAT)	aceA →	isocitrate lyase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,217,943	0	T	C	5.7%	290.9 / 6.2	87	Y279H (TAT→CAT)	aceA	isocitrate lyase
Reads supporting (aligned to +/- strand): ref base T (47/35); new base C (3/2); total (50/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.69e-01

CCCGTACCGATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCAGAAAAACCT  >  NC_000913/4217800‑4218091

cccGTACCGATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATgcg                                                                                                                                                 >  2:246116/1‑149 (MQ=255)

 ccGTACCGATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATgcgc                                                                                                                                                >  1:622917/1‑149 (MQ=255)

      ccGATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCatat                                                                                                                                           <  1:105059/149‑1 (MQ=255)

      ccGATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCatat                                                                                                                                           <  2:725139/149‑1 (MQ=255)

           gCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGa                                                                                                                                      <  4:171848/149‑1 (MQ=255)

            cTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGAc                                                                                                                                     <  4:434294/149‑1 (MQ=255)

             tGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATgc                                                                                                                                                  <  2:134166/135‑1 (MQ=255)

             tGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATgc                                                                                                                                                  >  1:134166/1‑135 (MQ=255)

                      ggATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACctggtctggt                                                                                                                           >  4:697246/1‑149 (MQ=255)

                      ggATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACctggtctggt                                                                                                                           <  4:258668/149‑1 (MQ=255)

                       gATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATTCTGACCTGGTCTGgtg                                                                                                                          <  3:484656/149‑1 (MQ=255)

                        aTCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGgtgt                                                                                                                         >  3:211390/1‑149 (MQ=255)

                         tCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGtgtg                                                                                                                        >  2:783682/1‑149 (MQ=255)

                         tCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGCATacgc                                                                                                                                                >  2:134583/1‑125 (MQ=255)

                         tCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGCATacgc                                                                                                                                                <  1:134583/125‑1 (MQ=255)

                            gATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGaaa                                                                                                                     >  2:345641/1‑149 (MQ=255)

                            gATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGaaa                                                                                                                     >  1:359790/1‑149 (MQ=255)

                             aTCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAAc                                                                                                                    <  4:243114/149‑1 (MQ=255)

                                  cTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTc                                                                                                                 <  3:222904/147‑1 (MQ=255)

                                  cTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTc                                                                                                                 >  4:222904/1‑147 (MQ=255)

                                   tCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCat                                                                                                                                             >  2:311485/1‑118 (MQ=255)

                                   tCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCat                                                                                                                                             <  1:311485/118‑1 (MQ=255)

                                   tCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCAc                                                                                                              >  1:63055/1‑149 (MQ=255)

                                    ccGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACg                                                                                                             >  2:742470/1‑149 (MQ=255)

                                      gATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGcc                                                                                                           >  4:106745/1‑149 (MQ=255)

                                      gATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGcc                                                                                                           >  3:649759/1‑149 (MQ=255)

                                      gATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGcc                                                                                                           >  2:226906/1‑149 (MQ=255)

                                        ttGCGACCCGTATGACAGCGAATTTCTTACCGGCGAGCGTACCAGTGAAGGCTTCGTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCgg                                                                                                         >  4:673118/1‑149 (MQ=255)

                                            gACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCt                                                                                                     <  1:783682/149‑1 (MQ=255)

                                            gACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCt                                                                                                     >  1:784381/1‑149 (MQ=255)

                                              cccGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTgg                                                                                                   <  4:649759/149‑1 (MQ=255)

                                               ccGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGa                                                                                                  >  4:787541/1‑149 (MQ=255)

                                                cGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGaa                                                                                                 <  1:598527/149‑1 (MQ=255)

                                                       cAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGgcgc                                                                                          <  3:787541/149‑1 (MQ=255)

                                                       cAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGgcgc                                                                                          <  1:246590/149‑1 (MQ=255)

                                                        aGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGcgcg                                                                                         >  2:744556/1‑149 (MQ=255)

                                                         gCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCt                                                                                                     >  1:760178/1‑136 (MQ=255)

                                                         gCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCt                                                                                                     <  2:760178/136‑1 (MQ=255)

                                                             aTTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCt                                                                                    <  2:63055/149‑1 (MQ=255)

                                                                 attaCCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGc                                                                                <  3:572942/149‑1 (MQ=255)

                                                                   taCCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGcac                                                                              >  1:509562/1‑149 (MQ=255)

                                                                   taCCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGcac                                                                              >  1:147837/1‑149 (MQ=255)

                                                                     ccGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACaa                                                                            <  3:270903/149‑1 (MQ=255)

                                                                      cGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAg                                                                           <  2:622917/149‑1 (MQ=255)

                                                                       ggCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGc                                                                          >  1:369498/1‑149 (MQ=255)

                                                                          gAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTgg                                                                                             <  1:356697/127‑1 (MQ=255)

                                                                          gAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTgg                                                                                             >  2:356697/1‑127 (MQ=255)

                                                                                aCCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCAcgc                                                                 >  2:688119/1‑149 (MQ=255)

                                                                                  cAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGa                                                               >  1:667159/1‑149 (MQ=255)

                                                                                  cAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGa                                                               >  1:265137/1‑149 (MQ=255)

                                                                                  cAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGa                                                               >  1:172534/1‑149 (MQ=255)

                                                                                       aaGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAatat                                                          >  3:337919/1‑149 (MQ=255)

                                                                                       aaGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAatat                                                          <  1:246116/149‑1 (MQ=255)

                                                                                       aaGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAatat                                                          <  1:815359/149‑1 (MQ=255)

                                                                                           cttcttCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCgg                                                      <  4:324103/149‑1 (MQ=255)

                                                                                           cttcttCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCgg                                                      <  4:365572/149‑1 (MQ=255)

                                                                                               ttcCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCaa                                                  >  4:159723/1‑149 (MQ=255)

                                                                                                tcCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCaaa                                                 <  2:784381/149‑1 (MQ=255)

                                                                                                  cGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAAct                                               <  1:688119/149‑1 (MQ=255)

                                                                                                      cTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAActgctg                                           >  4:706166/1‑149 (MQ=255)

                                                                                                       tCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTgg                                          >  1:380200/1‑149 (MQ=255)

                                                                                                         aTGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGc                                                                                <  2:416913/109‑1 (MQ=255)

                                                                                                         aTGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGc                                                                                >  1:416913/1‑109 (MQ=255)

                                                                                                           gCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTta                                      <  3:706166/149‑1 (MQ=255)

                                                                                                             gggCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGAt                                                                                                       >  4:716141/1‑82 (MQ=255)

                                                                                                             gggCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGAt                                                                                                       <  3:716141/82‑1 (MQ=255)

                                                                                                             gggCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCt                                                                                    >  2:111480/1‑101 (MQ=255)

                                                                                                             gggCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCt                                                                                    <  1:111480/101‑1 (MQ=255)

                                                                                                                 aTTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTg                                >  2:586121/1‑149 (MQ=255)

                                                                                                                 aTTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTg                                >  2:328651/1‑149 (MQ=255)

                                                                                                                    gAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATaa                                   <  1:105410/143‑1 (MQ=255)

                                                                                                                    gAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATaa                                   >  2:105410/1‑143 (MQ=255)

                                                                                                                      gCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGc                           <  3:159723/149‑1 (MQ=255)

                                                                                                                       cAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGcc                          >  1:319810/1‑149 (MQ=255)

                                                                                                                       cAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATACGGGCAAACTGCTGGCTTATAACTGCTCGcc                          >  2:769285/1‑149 (MQ=255)

                                                                                                                        aaGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCcg                         >  1:815010/1‑149 (MQ=255)

                                                                                                                          gCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCcgtc                       >  1:284550/1‑149 (MQ=255)

                                                                                                                            gATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCcgtcgt                     >  2:307003/1‑149 (MQ=255)

                                                                                                                            gATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCcgtcgt                     >  4:638776/1‑149 (MQ=255)

                                                                                                                             aTCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGtt                    >  2:766629/1‑149 (MQ=255)

                                                                                                                             aTCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGtt                    >  3:18532/1‑149 (MQ=255)

                                                                                                                                  ccGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACt               >  4:64213/1‑149 (MQ=255)

                                                                                                                                   cGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTg              <  2:369498/149‑1 (MQ=255)

                                                                                                                                     gcatatGGCGCATACGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCgg                                                      >  3:429849/6‑107 (MQ=255)

                                                                                                                                     gcatatGGCGCATACGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCgg                                                      <  4:429849/102‑1 (MQ=255)

                                                                                                                                       atatGGCGCATACGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCAg          >  3:774776/4‑149 (MQ=255)

                                                                                                                                             cGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCAGAAAAAc    <  2:380200/149‑1 (MQ=255)

                                                                                                                                               tATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCAGAAAAACCt  >  2:446214/1‑149 (MQ=255)

                                                                                                                                               tATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCAGAAAAACCt  <  1:226906/149‑1 (MQ=255)

                                                                                                                                               tATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGACGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCAGAAAAACCt  >  3:116708/1‑149 (MQ=255)

CCCGTACCGATGCTGATGCGGCGGATCTGATCACCTCCGATTGCGACCCGTATGACAGCGAATTTATTACCGGCGAGCGTACCAGTGAAGGCTTCTTCCGTACTCATGCGGGCATTGAGCAAGCGATCAGCCGTGGCCTGGCGTATGCGCCATATGCTGACCTGGTCTGGTGTGAAACCTCCACGCCGGATCTGGAACTGGCGCGTCGCTTTGCACAAGCTATCCACGCGAAATATCCGGGCAAACTGCTGGCTTATAACTGCTCGCCGTCGTTCAACTGGCAGAAAAACCT  >  NC_000913/4217800‑4218091

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