BRESEQ :: Evidence

Predicted mutation
evidence	seq id	position	mutation	freq	annotation	gene	description
RA	NC_000913	4,261,155	A→T	48.2%	D186V (GAT→GTT)	yjbM →	uncharacterized protein

Read alignment evidence...
	seq id	position		ref	new	freq	score (cons/poly)	reads	annotation	genes	product
*	NC_000913	4,261,155	0	A	T	48.2%	‑1.1 / 283.4	185	D186V (GAT→GTT)	yjbM	uncharacterized protein
Reads supporting (aligned to +/- strand): ref base A (43/53); new base T (39/50); total (82/103)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.33e-01

Only 100 of 187 total aligned reads displayed.

AGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACTTA  >  NC_000913/4261056‑4261253

agTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGa                                                                                                    >  1:517447/1‑100 (MQ=255)

 gTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGAt                                                                                                   <  1:4161175/100‑1 (MQ=255)

   aCCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTcc                                                                                                 <  1:173525/100‑1 (MQ=255)

   aCCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTcc                                                                                                 <  2:903179/100‑1 (MQ=255)

    ccTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCt                                                                                                <  2:3072284/100‑1 (MQ=255)

    ccTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCt                                                                                                >  1:663626/1‑100 (MQ=255)

     cTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTg                                                                                               <  1:4182507/100‑1 (MQ=255)

     cTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTg                                                                                               >  2:3015516/1‑100 (MQ=255)

      tGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGa                                                                                              >  2:1189780/1‑100 (MQ=255)

       gggAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGat                                                                                             <  2:3562074/100‑1 (MQ=255)

       gggAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGat                                                                                             >  1:1084111/1‑100 (MQ=255)

       gggAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGat                                                                                             <  1:3071482/100‑1 (MQ=255)

         gagagaGGAAGCGTCTGGCTGATGGGGGAGAAATAAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGatat                                                                                           <  1:2694918/100‑1 (MQ=255)

           gagaTGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATAttt                                                                                         <  2:3715992/100‑1 (MQ=255)

           gagaTGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATAttt                                                                                         <  2:2469936/100‑1 (MQ=255)

           gagaTGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATAttt                                                                                         <  1:2347669/100‑1 (MQ=255)

             gaTGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTa                                                                                       <  1:2621026/100‑1 (MQ=255)

               tGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAg                                                                                     <  2:1378700/100‑1 (MQ=255)

               tGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGATCGTTCCTGATATTTTAAg                                                                                     <  1:2488085/100‑1 (MQ=255)

                gAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGc                                                                                    >  1:1065861/1‑100 (MQ=255)

                 aaGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCa                                                                                   >  1:4220252/1‑100 (MQ=255)

                  ngCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCaa                                                                                  >  2:1597542/2‑100 (MQ=255)

                  aGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCaa                                                                                  >  2:1541371/1‑100 (MQ=255)

                  aGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCaa                                                                                  >  2:2298771/1‑100 (MQ=255)

                     gTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAAt                                                                               >  2:3117371/1‑100 (MQ=255)

                         ggCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCga                                                                           >  2:664722/1‑100 (MQ=255)

                          gCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCgag                                                                          >  1:797000/1‑100 (MQ=255)

                          gCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCgag                                                                          <  2:630500/100‑1 (MQ=255)

                           cTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCgaga                                                                         >  2:2232870/1‑100 (MQ=255)

                           cTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCgaga                                                                         <  1:533272/100‑1 (MQ=255)

                              aTGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGAtgc                                                                      <  2:4122544/100‑1 (MQ=255)

                                 ggggAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGAtgctgc                                                                   >  1:97857/1‑100 (MQ=255)

                                 ggggAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGAtgctgc                                                                   <  2:2911889/100‑1 (MQ=255)

                                   ggAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGcac                                                                 <  1:1075135/100‑1 (MQ=255)

                                   ggAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGcac                                                                 >  1:3829495/1‑100 (MQ=255)

                                   ggAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGcac                                                                 <  2:4218021/100‑1 (MQ=255)

                                         aTGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaat                                                           >  2:2015879/1‑100 (MQ=255)

                                           gAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaatta                                                         >  1:3157969/1‑100 (MQ=255)

                                           gAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaatta                                                         <  2:847205/100‑1 (MQ=255)

                                             ccgatcttGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaattaat                                                       <  1:1445341/93‑1 (MQ=255)

                                             aTTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaattaat                                                       >  1:2172279/1‑100 (MQ=255)

                                             aTTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaattaat                                                       >  1:967/1‑100 (MQ=255)

                                              ttttGATGGCTATGTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaattaat                                                       <  2:839273/99‑1 (MQ=255)

                                              ttttGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaattaatt                                                      >  2:1430037/1‑100 (MQ=255)

                                               tttGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAattaatta                                                     <  1:1792465/100‑1 (MQ=255)

                                                ttaaTGGCGATCTGGATAAGATGACACCAGTGTTGAATGAAGGAACGCACAATCCTGAAATTTAAAGCAAAATGCAAGATGCTGCACCTCaattaatt                                                      <  1:885481/95‑1 (MQ=255)

                                                ttGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAg                                                    <  2:1743272/100‑1 (MQ=255)

                                                ttGATGGCTATCTGGATAAGATGTAATATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAattaatta                                                     <  1:1584443/99‑1 (MQ=255)

                                                 tGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGc                                                   <  2:767462/100‑1 (MQ=255)

                                                 tGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGc                                                   <  2:2102267/100‑1 (MQ=255)

                                                 tGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGc                                                   <  2:2034495/100‑1 (MQ=255)

                                                   aTGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCtt                                                 <  1:204547/100‑1 (MQ=255)

                                                    tGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACaattaatagatcgg                                                >  2:1445345/1‑93 (MQ=255)

                                                      gCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTcgc                                              >  1:3113370/1‑100 (MQ=255)

                                                       cTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTcgcg                                             <  2:261434/100‑1 (MQ=255)

                                                       cTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTcgcg                                             >  2:3229357/1‑100 (MQ=255)

                                                        tATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGa                                            >  1:4199160/1‑100 (MQ=255)

                                                        tATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGa                                            >  1:2648390/1‑100 (MQ=255)

                                                         aTCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGaa                                           <  2:2436739/100‑1 (MQ=255)

                                                               aTAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGtt                                     <  1:2655822/100‑1 (MQ=255)

                                                                 aaGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAATGCGTTTg                                   <  1:351115/100‑1 (MQ=255)

                                                                  aGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGa                                  >  2:33846/1‑100 (MQ=255)

                                                                  aGATGTCACATATGGTGAATGAAGGAACGCTTGATCCTAAAAAAAAAAAAAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGa                                  <  1:1394118/100‑1 (MQ=255)

                                                                   gATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGttt                                    <  1:2255450/97‑1 (MQ=255)

                                                                    aTGTCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGAcc                                >  2:3312304/1‑100 (MQ=255)

                                                                    aTGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATTCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGAcc                                <  1:2371513/100‑1 (MQ=255)

                                                                    aTGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGAcc                                >  2:522653/1‑100 (MQ=255)

                                                                       tCACATCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCaa                             >  2:2735900/1‑100 (MQ=255)

                                                                        cacaTCTGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAAc                            >  1:3191235/1‑100 (MQ=255)

                                                                        cacaTCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAAc                            >  2:87169/1‑100 (MQ=255)

                                                                        cacaTCTGGTGAATGAAGGAACGCTCGACCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACACTTAATTAGCTTCGCGAAGTCGTTTGACCCAAc                            >  1:3971322/1‑100 (MQ=255)

                                                                         acaTCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAAcc                           <  1:1464626/100‑1 (MQ=255)

                                                                           aTCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTc                         >  2:1087634/1‑100 (MQ=255)

                                                                            tCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCa                        <  2:81352/100‑1 (MQ=255)

                                                                              tGGTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCaaa                      <  2:2909552/100‑1 (MQ=255)

                                                                              tGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAAGGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCaaa                      <  1:1724518/100‑1 (MQ=255)

                                                                               taaaatatcaGGAACGAGCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATGAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAg                     >  1:2887920/10‑100 (MQ=255)

                                                                               ggTGAATGAAGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAg                     <  2:2768827/100‑1 (MQ=255)

                                                                                gTGAATGAAGGAACGCTCGTTCCTGCTATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAgg                    <  1:1327711/100‑1 (MQ=255)

                                                                                  gaatgaaGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGaa                  >  2:1211492/1‑100 (MQ=255)

                                                                                  gaatgaaGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGaa                  >  1:3029906/1‑100 (MQ=255)

                                                                                   aatgaaGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAag                 <  2:3275661/100‑1 (MQ=255)

                                                                                    atgaaGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAaga                >  2:2230317/1‑100 (MQ=255)

                                                                                     tgaaGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAagag               >  2:3784464/1‑100 (MQ=255)

                                                                                     tgaaGGAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAagag               <  2:573752/100‑1 (MQ=255)

                                                                                     atcaGGAACGAGCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAagag               >  1:1309675/4‑100 (MQ=255)

                                                                                       aaGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGAt             >  2:1328759/1‑100 (MQ=255)

                                                                                       aaGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGAt             <  2:3168220/100‑1 (MQ=255)

                                                                                        aGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGAtt            >  2:396224/1‑100 (MQ=255)

                                                                                         ggAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTa           >  2:3605667/1‑100 (MQ=255)

                                                                                          gAACGCTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTaa          <  2:437013/100‑1 (MQ=255)

                                                                                          gAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTaa          <  1:1521872/100‑1 (MQ=255)

                                                                                            aCGCTCGATCCTGATATTTTAAGCAAGATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTaaa         <  1:1909582/99‑1 (MQ=255)

                                                                                               cTCGTTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATAc     <  2:3540919/100‑1 (MQ=255)

                                                                                               cTCGTTCCAGATATAATAAACAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATa      <  1:2195159/99‑1 (MQ=255)

                                                                                                tCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAgagaagatcggaag    >  1:3152553/1‑90 (MQ=255)

                                                                                                tCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACt    <  2:417407/100‑1 (MQ=255)

                                                                                                tCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACt    <  1:745187/100‑1 (MQ=255)

                                                                                                  naTCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTGGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACTTa  >  1:2327172/2‑100 (MQ=255)

                                                                                                  gATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACTTa  <  2:4040044/100‑1 (MQ=255)

AGTACCTGGGAGAGATGAAGCGTCTGGCTGATGGGGGAGAAATGAATTTTGATGGCTATCTGGATAAGATGTCACATCTGGTGAATGAAGGAACGCTCGATCCTGATATTTTAAGCAAAATGCGAGATGCTGCACCACAATTAATTAGCTTCGCGAAGTCGTTTGACCCAACCTCAAAGGAAGAGATTAAAATACTTA  >  NC_000913/4261056‑4261253

Alignment Legend
Aligned base mismatch/match (shaded by quality score): `ATCG/ATCG` `< 3 ≤` `ATCG/ATCG` `< 12 ≤` `ATCG/ATCG` `< 32 ≤` `ATCG/ATCG` `< 41 ≤` `ATCG/ATCG`
Unaligned base: `atcg` Masked matching base: `atcg` Alignment gap: `‑` Deleted base: `‑`