BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	499,459	A→C	5.4%	L172R (CTC→CGC) ‡	aes ←	acetyl esterase

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	499,459	0	A	C	5.4%	452.6 / 13.6	130	L172R (CTC→CGC) ‡	aes	acetyl esterase
Reads supporting (aligned to +/- strand): ref base A (59/63); new base C (3/4); total (62/68)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.15e-01

Only 100 of 130 total aligned reads displayed.

AAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTG  >  NC_000913/499312‑499605

aaaTCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAg                                                                                                                                                   >  1:131001/1‑149 (MQ=255)

    ccTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGAGAGCGcc                                                                                                                                               <  2:587820/149‑1 (MQ=255)

       gTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAg                                                                                                                                             >  3:889625/1‑148 (MQ=255)

         tGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGCTCGCCAgtg                                                                                                                                           <  4:39526/148‑3 (MQ=255)

          gCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATg                                                                                                                                         <  2:872046/149‑1 (MQ=255)

            gTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGc                                                                                                                                       <  2:426260/149‑1 (MQ=255)

              tAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCa                                                                                                                                      >  4:540691/1‑148 (MQ=255)

              tAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCa                                                                                                                                      >  3:994456/1‑148 (MQ=255)

              tAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCa                                                                                                                                      >  3:81769/1‑148 (MQ=255)

                aGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCt                                                                                                                                   >  2:604176/1‑149 (MQ=255)

                 gCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTg                                                                                                                                  >  1:35014/1‑149 (MQ=255)

                 gCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTg                                                                                                                                  >  2:449787/1‑149 (MQ=255)

                  ccATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTg                                                                                                                                  >  3:467233/1‑148 (MQ=255)

                        cAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAAt                                                                                                                            <  4:977308/148‑1 (MQ=255)

                           acacCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCAcc                                                                                                                        <  2:1025007/149‑1 (MQ=255)

                            cacCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATg                                                                                                                                         >  4:202944/1‑131 (MQ=255)

                            cacCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATg                                                                                                                                         <  3:202944/131‑1 (MQ=255)

                            cacCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCAcc                                                                                                                        <  3:317072/148‑1 (MQ=255)

                              ccGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGc                                                                                                                     >  1:457083/1‑149 (MQ=255)

                                gCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCaa                                                                                                                   >  1:434748/1‑149 (MQ=255)

                                 cccAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCaa                                                                                                                   >  4:1039231/1‑148 (MQ=255)

                                   cAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAACCGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGc                                                                                                                <  1:355497/149‑1 (MQ=255)

                                      cAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCTACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCaa                                                                                                             <  1:768970/149‑1 (MQ=255)

                                         acgacgAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGc                                                                                                          <  2:150498/149‑1 (MQ=255)

                                          cgacgaGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAAAGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGc                                                                                                          >  4:179473/1‑148 (MQ=255)

                                           gacgaGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCgg                                                                                                        <  1:762298/149‑1 (MQ=255)

                                           gacgaGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCgg                                                                                                        <  1:292260/149‑1 (MQ=255)

                                             cgaGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCggg                                                                                                       <  3:267805/148‑1 (MQ=255)

                                             cgaGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCggg                                                                                                       <  4:223745/148‑1 (MQ=255)

                                                gTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACa                                                                                                    <  4:331570/148‑1 (MQ=255)

                                                gTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACa                                                                                                    <  4:842107/148‑1 (MQ=255)

                                                   aCGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGGACTGGCGGGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACatat                                                                                                 <  4:194246/148‑1 (MQ=255)

                                                    cGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATAtt                                                                                                >  4:707706/1‑148 (MQ=255)

                                                     ggAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATAttg                                                                                               >  4:1002947/1‑148 (MQ=255)

                                                          cccGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATAttgatttg                                                                                          <  3:70377/148‑1 (MQ=255)

                                                          cccGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATAttgattag                                                                                          <  4:218140/148‑3 (MQ=255)

                                                          cccGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATAttgacttg                                                                                          <  3:70361/148‑1 (MQ=255)

                                                            cGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATa                                                                                       <  1:64182/149‑1 (MQ=255)

                                                            cGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATa                                                                                       <  2:469558/149‑1 (MQ=255)

                                                              tAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAAt                                                                                     >  1:903712/1‑149 (MQ=255)

                                                                 tCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAAtcct                                                                                  >  1:934455/1‑149 (MQ=255)

                                                                  ccGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAAtcctc                                                                                 <  1:663629/149‑1 (MQ=255)

                                                                  ccGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAAtcct                                                                                  <  4:973237/148‑1 (MQ=255)

                                                                    gTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCg                                                                               >  2:160688/1‑149 (MQ=255)

                                                                         gCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCctg                                                                           <  3:1039231/148‑1 (MQ=255)

                                                                             aTACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCtggt                                                                      >  1:269686/1‑149 (MQ=255)

                                                                                ccACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGa                                                                   >  1:508296/1‑149 (MQ=255)

                                                                                  acaGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCGCTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGaa                                                                  <  3:700734/148‑1 (MQ=255)

                                                                                        aaaCGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATaacaa                                                           >  2:652761/1‑149 (MQ=255)

                                                                                         aaCGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATaacaa                                                           >  3:897409/1‑148 (MQ=255)

                                                                                           cGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAacaaca                                                         <  3:772364/148‑1 (MQ=255)

                                                                                           cGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAacaaca                                                         <  3:776198/148‑1 (MQ=255)

                                                                                             cccGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAagc                                                      <  1:610682/149‑1 (MQ=255)

                                                                                              ccGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCAc                                                                                                                         <  1:542075/81‑1 (MQ=255)

                                                                                              ccGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCAc                                                                                                                         >  2:542075/1‑81 (MQ=255)

                                                                                                gCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAagcagc                                                   <  1:919356/149‑1 (MQ=255)

                                                                                                gCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAagcag                                                    >  4:629514/1‑148 (MQ=255)

                                                                                                 cAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAagcagc                                                   >  3:443817/1‑148 (MQ=255)

                                                                                                         ccGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTc                                          <  2:457083/149‑1 (MQ=255)

                                                                                                          cGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTcc                                         >  2:1072617/1‑149 (MQ=255)

                                                                                                          cGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTcc                                         <  2:903712/149‑1 (MQ=255)

                                                                                                          cGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTcc                                         <  2:825245/149‑1 (MQ=255)

                                                                                                           gCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCct                                        <  1:160688/149‑1 (MQ=255)

                                                                                                              aTCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTa                                     >  1:23146/1‑149 (MQ=255)

                                                                                                               tCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTAt                                    >  2:152537/1‑149 (MQ=255)

                                                                                                                  aTCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGc                                 >  1:392595/1‑149 (MQ=255)

                                                                                                                     tGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCtggt                                                                      <  2:430774/109‑1 (MQ=255)

                                                                                                                     cggcaaTCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCtt                               >  4:407746/6‑148 (MQ=255)

                                                                                                                      gTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGc                             <  1:449787/149‑1 (MQ=255)

                                                                                                                           tCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGc                             >  2:871392/1‑144 (MQ=255)

                                                                                                                           tCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGc                             <  1:871392/144‑1 (MQ=255)

                                                                                                                           tCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGaa                         >  3:447356/1‑148 (MQ=255)

                                                                                                                           tCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGaa                         >  4:458126/1‑148 (MQ=255)

                                                                                                                            cACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGaaaa                       <  2:818555/149‑1 (MQ=255)

                                                                                                                            cACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGaaa                        <  3:550536/148‑1 (MQ=255)

                                                                                                                            cACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGaaa                        >  3:882086/1‑148 (MQ=255)

                                                                                                                            cACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGaaa                        >  3:393097/1‑148 (MQ=255)

                                                                                                                             aCGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAAc                      >  2:29030/1‑149 (MQ=255)

                                                                                                                              cGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAAc                      >  4:1058818/1‑148 (MQ=255)

                                                                                                                              cGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAAc                      >  3:650850/1‑148 (MQ=255)

                                                                                                                               gCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAAcgc                    >  1:18848/1‑149 (MQ=255)

                                                                                                                                 aaCCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAagcag                                                    <  2:977525/115‑1 (MQ=255)

                                                                                                                                 aaCCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAagcag                                                    >  1:977525/1‑115 (MQ=255)

                                                                                                                                 aaCCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACgcgc                  <  2:396285/149‑1 (MQ=255)

                                                                                                                                    cacaACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCtt                <  3:313986/148‑1 (MQ=255)

                                                                                                                                    cacaACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTc               >  2:773910/1‑149 (MQ=255)

                                                                                                                                    cacaACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCAGCctgc                                                                          >  1:858436/1‑90 (MQ=255)

                                                                                                                                      caACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCGCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAg             <  1:837594/149‑1 (MQ=255)

                                                                                                                                       aaCGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGaaaa                       >  2:1057931/1‑138 (MQ=255)

                                                                                                                                       aaCGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGaaaa                       <  1:1057931/138‑1 (MQ=255)

                                                                                                                                        aCGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAAcgcg                   <  4:378351/141‑1 (MQ=255)

                                                                                                                                         cGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGt           >  3:1047935/1‑148 (MQ=255)

                                                                                                                                           tgcTGGCGCTCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTc                                          >  4:276954/3‑115 (MQ=255)

                                                                                                                                           tgcTGGCGCTCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTc                                          <  3:276954/113‑1 (MQ=255)

                                                                                                                                           tgcTGGCGCTCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAAcgcg                   >  4:31862/3‑138 (MQ=255)

                                                                                                                                           tgcTGGCGCTCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAAcgcg                   <  3:31862/136‑1 (MQ=255)

                                                                                                                                              tGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAg      <  3:943010/148‑1 (MQ=255)

                                                                                                                                              tGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAg      <  3:849646/148‑1 (MQ=255)

                                                                                                                                              tGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAg      <  3:622237/148‑1 (MQ=255)

                                                                                                                                              tGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAg      <  3:179473/148‑1 (MQ=255)

AAATCCTGTTGCGTTAAGCCATCCCAGACACCGCCCAACAGACGACGAGTCACGGAATCCCGTAATCCGTAAAGCCCATACCACAGCAAAACGCCCGCAACTTTACCGCAATCGATCTGTTTATCACGCAACCACAACGCACTGGCGAGCGCCAGCATGGCACCTGCGGAATCACCGGCAAAGCCAATGCGGGACATATTGATTTGATAATCCTCCGCCTGCTGGTGGAAATAACAACAAGCAGCCACAATTTCCTCTATCGCTTGCGGAAAACGCGCTTCAGGTGAAAGGGTG  >  NC_000913/499312‑499605

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